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news and views 


ADMINISTRATIVE COMMITTEE 


Henry Booker and John Smyth have been re-elected 
to the Administrative Committee of the PGAP, and 
J. I. Bohnert and D. D. King have been elected. These 
are for three year terms beginning June 1, 1955. Henry 
Booker has also been elected Vice-Chairman for the 
coming year. 


PGAP SPONSORSHIP AND PARTICIPATION IN 
NATIONAL TECHNICAL SYMPOSIA 


The Administrative Committee, under the effective 
jurisdiction of Mr. Bolljahn, is evolving a policy for 
PGAP’s role in symposia. It has been agreed that we 
plan to take a very positive active part in the symposia 
of which we consider ourselves co-sponsor, and that in 
other meetings where we take a subordinate role and 
assist only by such things as contacting our membership 
‘and other minor promotional activities, we should con- 
sider ourselves as only participants, and not be classified 
as a co-sponsor. 

As a tentative program for the coming year, it has 
been agreed that we should take an active part in the 
IRE National Convention in New York. We would also 
plan to sponsor sessions at the National Electronics 
Convention and the WESCON. Plans for this are al- 
ready well under way. It has also been planned that, 
subject to acceptance by the appropriate URSI Com- 
missions, we should take part as co-sponsor in the 
Spring and Fall National URSI Meetings. In addition 
to these, we should sponsor other national meetings 
where our group would have a definite, positive con- 
‘tribution to make. 

We have received an invitation from W. W. Mum- 
ford, Chairman of PGMTT, inviting us to participate 
and co-sponsor a technical symposium in Philadelphia 
on February 2 and 3, 1956. The theme of this session is 
to be Microwave Techniques. It has been agreed that we 
will accept this invitation and some PGAP representa- 
tives will be named to help prepare the arrangements 
and program. 


Mr. Bolljahn has requested suggestions for interesting 
session titles. Roy Spencer has suggested “Fourier and 
Information Theory Techniques Applied to Antenna 
Scanning.” Please forward your suggestions to Mr. 
Bolljahn. 

For the National Convention, it has been suggested 
that the general format for the contributed paper ses- 
sions be changed to accommodate more papers. The 
general feeling is that, when necessary, authors could 
reduce the length of their papers without sacrificing 
pertinent details, and that discussion could be reduced 
to a minimum. The general plan is to allow 15 minutes 
for each paper with no discussion following the paper. 
However, the authors would be requested to remain and 
written questions could be submitted to the chairman 
who would address them to the appropriate author dur- 
ing the general discussion period at the end of the ses- 
sion. The papers of general interest, such as new applica- 
tions, would be emphasized and the authors of more 
scientific papers of significance would be encouraged to 
submit their papers for the URSI Sessions. D. C. Ports 
is advising the Convention Committee of our plans in 
this regard. 


PGAP NEws 


An annual index to the 1955 volume of IRE 
TRANSACTIONS on Antennas and Propagation will 
appear in the January, 1956 issue in order that it 
may include the Proceedings of the URSI-Michi- 
gan Symposium, which will be published shortly. 

—The Editor 


The July issue of TRANSACTIONS contained abstracts 
of the papers presented at the IRE-URSI Symposium 
held May 2-5, 1955, in Washington D. C. The PGAP, 
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as co-sponsor of the USA URSI Meetings, will hence- 
forth publish in the July issue of the TRANSACTIONS ab- 
stracts of the papers presented, thus providing a means 
of distributing this information to those who did not 
attend the sessions. This will serve as a permanent ref- 
erence. 

The Washington Chapter of PGAP is quite active and 
has formulated an ambitious plan for a technical sym- 
posium on “Communication by Scatted Techniques,” 
November 14 and 15, 1955. It is being co-sponsored by 
the Professional Group on Communication Systems and 
George Washington University, who will provide the 
auditorium and other meeting facilities. The general 
theme is on application rather than theoretical aspects. 


CHAPTER NEWS 


Washington, D. C.—At its May 23, 1955 meeting, the 
following officers were elected for the coming year: 


Chairman Coleman Goatley 
Vice-Chairman Clarence H. Stewart 
Secretary William W. Balwanz 


A paper entitled “Behavior of Ferrites in Microwave 
Components” was presented by John C. Cacheris. 


Los Angeles—The summer meeting of the Los Angeles 
Chapter of the PGAP, held on August 11 at the I.A.S. 
Building, 7660 Beverly Boulevard, featured two out- 
standing speakers in the antenna field: V. H. Rumsey of 
the University of Illinois and L. C. Van Atta of Hughes 
Aircraft Company. The dinner before the meeting was 
held as usual at the Encore Restaurant, 806 N. La 
Cienega, with 54 members attending. A total of 106 
members and guests were present at the general meet- 
ing. 

Professor Rumsey, who is spending the summer at the 
Lockheed Missile Systems Division in Van Nuys, gave 
a highly stimulating talk on “Some Recent Develop- 
ments in Antennas” which included some very useful 
ideas in the design of antennas. Configurations specified 
only by angles were shown to have infinite bandwidth 
and some practical approximations were discussed. One 
interesting topic was the 60 Pi antenna resulting from a 
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configuration having a Babinet complement identical to 
the antenna itself and having always an input imped- 
ance of 60 Pi regardless of shape. 

Dr. Van Atta presented a summary of developments 
in surface wave antennas including two dimensional 
arrays of slots, printed dipoles, leaky waveguides, and 
wave trapping surfaces. Some experimental data were 
presented which showed the excellent beam shaping that 
can be obtained with corrugated surface wave antennas. 


CALENDAR OF COMING EVENTS 


The National Electronics Conference, Hotel Sher- 
man, Chicago, IIl., Oct. 3-5. 


AIEE Fall General Meeting, Morrison Hotel, 
Chicago, Ill., Oct. 3—7. 


AIEE Aircraft Electronic Equipment Conference, 
Los Angeles, Calif., Oct. 11-13. 


IRE East Coast Conference on Aeronautical & 
Navigational Electronics, Lord Baltimore Hotel, 
Baltimore, Md., Oct. 31—Nov. 1. 


World Symposium on Applied Solar Energy, Con- 
ducted under leadership of Standford Research 
Institute, Phoenix, Arizona, Oct. 31—Nov. 4. 


Nuclear Congress and Atomic Exposition, spon- 
sored by the Engineers Joint Council, Cleveland, 
Ohio, Dec. 12-17. 


AIEE Winter General Meeting, Statler Hotel, 
New York, N. Y., Jan. 30-Feb. 3, 1956. 


Symposium on Microwave Theory and Tech- 
niques, University of Pennsylvania, Phila- 
delphia, Penna, Feb. 2-3, 1956. 


The 34th Annual Convention of the National 
Association of Radio and Television Broad- 
casters, Conrad Hilton Hotel, Chicago, IIL, 
April 15-19, 1956. 


National Aeronautical and Navigational Electron- 
ics Conference, Dayton, Ohio, May 14-16, 1956. 


COS So 


1955 
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contributions 


Prediction of Oceanic Duct Propagation from 
Climatological Data* 
L. J. ANDERSON} anv E. E. GOSSARD} 


Summary—Techniques are described for predicting the effect of 
the “oceanic duct” in increasing the coverage of low-sited microwave 
radars. Prediction is made on a monthly basis using statistical data 
, on sea temperature, air temperature, humidity, and wind speed. A 
large area in the northwest Pacific is used to illustrate the methods 
used. The basic prediction approach is to deduce refractive index 
profiles from available data parameters and to use Schelkunoff’s 
ducting criteria to predict radar propagation effects. Monthly con- 
tours are then drawn over the area for each radar band showing the 
probability of extended coverage. 


INTRODUCTION 


HE EVENTUAL goal of all tropospheric propa- 
“ee research is the prediction of field strengths, 
at frequencies above 50 mc, from routine meteoro- 
logical data. For most situations where such predictions 
are desired the problem is sufficiently complex that pre- 
dictions require careful analysis by experienced person- 
nel. It should be stated parenthetically that such pre- 
dictions can be made and in fact have been made and 
verified on several occasions with surprising accuracy. 
In such cases, however, semi-empirical methods must be 
employed and cookbook routines cannot be specified. 
There is one important case, however, where a routine 
method has been developed which has greatly simplified 
the mechanics of preparing the predictions. This is the 
case of predicting the coverage of low-sited microwave 
radars as affected by the oceanic duct. The meteorologi- 
cal parameters used are routine shipboard observations 
of sea temperature, air temperature, and humidity at 
bridge height, and wind speed. The purpose of this paper 


* Original manuscript received by the PGAP, February 18, 1955; 


revised manuscript received, May 6, 1955. ; 
+ U.S. Navy Electronics Lab., San Diego 52, Calif. 


is to outline the method of reducing these data to sta- 
tistical radar coverage prediction charts for a large area 
of the northwest Pacific. 


BACKGROUND 


In order to assess the effect of the meteorological con- 
ditions upon propagation, one must work with vertical 
profiles of refractive index. Therefore, the first step in 
the prediction process is to deduce these profiles from 
the routine meteorological parameters enumerated 
above. Having the shape of the profile, one then must 
determine its effect upon the propagation of radar 
energy of a given wavelength. 

Details of these two steps are discussed in a previous 
paper! and will be outlined only briefly here. The tem- 
perature and humidity data allow one to compute the 
refractive index difference AM, between the sea surface 
and bridge height. Using the work of E. L. Deacon,? one 
can then deduce the shape of the profile using a profile 
index 6 which can be obtained from the air-sea tempera- 
ture difference and wind data. The wind data also en- 
ables one to pick a value for Zo, a length parameter char- 
acteristic of the roughness of the sea surface. Thus from 
the four simple observables listed in the introduction, 
the profile is determined. Detailed profiles* have been 
found to agree with those computed by the foregoing 
method. To assess the effect of a given profile upon 

1L. J. Anderson and E. E. Gossard, “The effect of the oceanic 
duct on microwave propagation,” Trans. Amer. Geophys. Union, 
vol. 34, p. 695; October, 1953. 

2 E. L. Deacon, “Vertical diffusion in the lowest layers of the at- 
mosphere,” Quart. Jour. R. Met. Soc., vol. 75, p. 89; January, 1949. 

3, J. Anderson and E. E. Gossard, “The Oceanic Duct and its 


Effect on Microwave Propagation,” U. S. Navy Electronics Lab., 
Rep. No. 305; August, 1952. 
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radar energy we use the approximate method of 
S. A. Schelkunoff.* ; 

Practical experience has shown that as trapping con- 
ditions are approached for a given wavelength, radar 
coverage is extended two to five times normal, whereas 
normal coverage is obtained from standard to almost 
trapping. The best one can predict seems to be whether 
radar coverage will be normal or extended two to five 
times. The gamut between these extremes is run with 
very slight changes in meteorological conditions. Inas- 
much as measured M profiles exhibit wide changes in 
both time and space one must work with time or space 
averages of the profiles to arrive at reliable results. 
Hence, one is hardly justified in attempting to refine the 
method to predict anything beyond the simple yes-no 
answer. When conditions are borderline between trap- 
ping and no trapping, they are so variable with time 
that further refinement would be unusable even if one 
could achieve it. 

With this philosophy in mind, one can solve for a 
minimum AM, required to just trap a given wavelength.! 
This yields: 


.060[3770\2(3 — B)/B |6/8 
3-8 ort 
2—-— B b1-8 = Zoe 


AM; (min) = 0.8, 


where 0 is the bridge height in feet and X is the wave- 
length in cm. This expression is the criterion used for 
predicting whether or not extended coverage should be 
experienced. It should also be stated that inherent in the 
method is the assumption that both the radar and the 
target are located within the oceanic duct, usually 
within 50 feet of the sea surface. This condition is met in 
the case of surface ship targets, but most shipboard 
radars are mounted as high as possible and only in the 
case of surfaced submarines is this below 50 feet. Fur- 
thermore, oceanic duct intensities normally encountered 
limit the wavelengths which are affected to 10 cm and 
shorter. On rare occasions 25 cm radars should also be 
affected. 


PREDICTION PROCEDURE 


Routine surface meteorological observations over the 
sea have been taken for many years and have been com- 
piled as monthly averages for 5 degree latitude and 
longitude squares in the “Atlas of Climatic Charts of the 
Oceans.”® Three such parameters that we need; namely, 
sea surface temperature, air temperature, and wet bulb 
depression have been tabulated and plotted in this way. 
In the case of wind speed, we have the percentage of 


4S. A. Schelkunoff, “Microwave Transmission in a Non-Homo- 
Seren Atmosphere,” Bell Tel. Lab. Rep. MM-44-110-53; July, 
5 “Atlas of Climatic Charts of the Oceans,” U. S. th 
W. B. No. 1247, U.S. Gov't. Printing Office: hee etn 
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observations where the wind is Beaufort Force 3 (9 
knots) or less. This is sufficient for determining the 
value of Zo to use. 

These average values are not sufficient to determine 
the percentage of time extended coverage should be en- 
countered since they tell us nothing about the frequency 
distribution over a given month of 6 and AM; all we get 
are monthly average values. In order to determine such 
distributions we must go to the individual ship observa- 
tions and determine 8 and AM, for each set and plot a 
histogram of a hundred or more such values. Figs. 1 
and 2 are histograms of 8 and AM, respectively for the 
Far East area in December. Plotting similar data for 
July and for the north and south extremes of the area 
shows that the shapes of the histograms are similar and 
are merely shifted along the abscissa. In Fig. 2 the histo- 
grams are shown for both the stable (6 <1) (solid curve) 
and the unstable (8>1) (dashed curve) cases and are 
seen to have virtually the same shapes. 

Using this fact we assume that since we know the 
average values of 6 and AM, for each month and each 
5 degrees square, we can simply shift the histogram 
curves along the abscissae until the average coincides 
with the monthly average and we will then have the re- 
quired histograms for the particular month and square. 

From the criterion equation we can draw curves, for 
each wavelength of interest, of AM; (min) vs 8. We must 
actually draw two such curves; one for the smooth sea 
case (Z)=0.00005 feet) and one for the rough sea case 
(Z) =0.02 feet). Furthermore, if the correlation between 
AM, and # is not significant we can obtain the prob- 
ability of AM, min (as given by the 8 vs AM, min 
curves) being exceeded for a given wavelength. There is 
actually some correlation between 6 and AM, but this 
has been partially taken into account by computing the 
probabilities separately for stable and unstable condi- 
tions. 

Fig. 3 illustrates the procedure. The numbers across 
the top of the figure are the percentages of observations 
of AM, in unstable air that fall in the class intervals 
shown along the bottom; the row of numbers across the 
middle of the figure (along the 8 =1.00 line) are similar 
percentages for stable air. The numbers just outside the 
left border are percentages of 8 observations in the class 
intervals indicated at the far left. One now multiplies 
each 8 percentage with each AM, percentage along its 
horizontal row and obtains a cross probability product. 
Starting at the right end of each row and working to- 
ward the left, adding the incremental cross probabilities 
as one proceeds, yields the numbers typed in the body 
of the grid. For the unstable cases (top half of the 
figure) one uses the AM, percentages at the top of the 
figure; for the lower half one uses those across the middle 
of the figure. Shifting histograms in this way gives a 
slight discontinuity in the vertical columns which is 
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the average AM, for the whole set of AM, values (both 
smoothed by computing the first row above and below _ stable and unstable) used in plotting the histograms of 
the discontinuity both ways and taking an average’ Figs. 1 and 2. 
value. In this case the differences were at most a few The criterion curves are now plotted as an underlay 
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Fig. 3—Cumulative probability product and X-band criteria, Far East, December. 
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Fig. 4--Extended coverage expectancy East Asia Coast, December X-band. 


for \=3 cm; the solid curve being for a rough sea sur- 
face and the dashed curve for a smooth sea surface. In 
order to determine the probability that 3 cm radars will 
encounter extended coverage for given average 8 and 
AM, values, one shifts the grid with its arrows along the 
AM axis until the AM arrow lies on the given average 
AM, and vertically until the 8 arrow falls on the given 
average 8 value. One then reads the number just to the 
right of the criterion curve in each row. The sum of these 
nine values gives the required probability in per cent. 
In practice these probabilities are read for both the 


smooth and rough sea cases and the final answer is the 
mean of the two weighted according to the per cent of 
time the wind is less than Beaufort Force 3. 

This process is repeated for each month and for each 
5 degree square in the area and the numbers plotted on 
a series of monthly charts. Contours are then drawn and 
yield charts similar to Figs. 4 and 5. Such charts have 
been prepared for each month for wavelengths of 3 cm, 
10 cm, and 25 cm.§ 


6“Far East Radio-Radar Propagation Conditionsfor Fleet 
Units,” U. S. Navy Electronics Rep. No. 319; July, 1952. 
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The Curved Passive Reflector* 
E. BEDROSIAN}{ 


Summary—The curved passive reflector is analyzed theoretically 
by the aperture-field method and the results are presented in the 
form of gain characteristics and radiation patterns. It is seen that 
curving the reflector produces a focusing effect in addition to the 
basic diffraction phenomenon, thereby providing increased gain as 
compared with the flat reflector. The limiting value for the flat re- 
flector gain is 6 db while the gain of the curved reflector increases 
continually with reflector size when properly illuminated. Radiation 
patterns are computed and it is shown that the effect of reducing 
tower height is to widen the main lobe and to move the side lobes 
outward and depress them. It is concluded that, in general, the 
curved reflector causes lower side lobes than the flat reflector. The 
dish illumination is shown to be relatively unimportant with respect 
to the total antenna gain but the side lobe levels are reduced when 
the dish illumination is made more uniform. 


INTRODUCTION 


ARLY microwave radio-relay systems employed 
iD antenna pairs mounted atop suitable towers 

which provided the terrain clearance between 
relay points. The relatively large tower heights required 
to allow reasonable path lengths made the placing of 
receiving and transmitting equipment at the top of a 
tower difficult and expensive while the base location led 
to inefficiencies because of transmission-line losses. The 
solution employed in modern microwave systems is the 
use of a vertically directed antenna at the base illumi- 
nating a passive reflector mounted at the top of the 
tower. 

The reflector is usually a flat elliptical sheet of con- 
ducting material mounted at a 45-degree angle and 
illuminated by a parabolic dish. Such a reflector has 
been the subject of a number of recent articles. Unger! 
assumed a typical illumination of the dish aperture and 
performed the diffraction analysis graphically. Jakes? 
performed essentially the same analysis using analog 
computation techniques assuming a 10 db-tapered 
parabolic dish illumination. His work was verified ex- 
perimentally by Drexler.’ Greenquist and Orlando! re- 
ported a theoretical and experimental analysis of a 
rectangular reflector based on a method of applying 
successive corrections to the far-field approximation 
assuming a dish illumination of the form (1—p?)!?. 
Finally, Crosby® analyzed the flat reflector using an ap- 
proximate graphical procedure with a 10 db-tapered 


* Original manuscript received by the PGAP, February 23, 1955. 

+ Motorola Inc., Riverside, Calif. 

1H. G. Unger, “Ebene spiegel zur strahlumlenkung bei richtan- 
tenen,” Frequenz, vol. 6, pp. 272-278; September, 1952. 

*W. C. Jakes, Jr., “A theoretical study of an antenna-reflector 
problem,” Proc. IRE, vol. 41, pp. 272-274; February, 1953. 

J. Drexler, “An experimental study of a microwave periscope,” 
Proc. IRE, vol. 42, p. 1022; June, 1954. 

4R. E. Greenquist and A. J. Orlando, “An analysis of passive re- 
rors antenna systems,” Proc. IRE, vol. 42, pp. 1173-1178; July, 

5 D. R. Crosby, “Theoretical Gain of Flat Microwave Reflector,” 
1954 IRE Convention ReEcorp, Part 1, “Antennas and Propaga- 
tion,” pp. 71-75. 


parabolic illumination. These works are in excellent 
agreement and show that the flat reflector can yield an 
increase in the far-field signal strength over that due to 
the parent antenna alone. 

The curved reflector is a natural outgrowth of its flat 
predecessor. The flat reflector yields an increased far- 
field signal strength because of the Fresnel-ring type of 
interference obtained by illuminating properly shaped 
apertures. The preceding analyses show the limiting 
increase for such a system to be 6 db. The curved re- 
flector utilizes not only the diffraction effect but also its 
natural focussing ability to yield an additional increase. 


THE DIFFRACTION PROBLEM 


A comprehensive treatment of the methods available 
for computing the diffracted field is given by Silver.® The 
current-distribution method used by Crosby assumes a 
current distribution on the reflector surface such that 
the incident field which would appear at the reflector 
surface in its absence is properly terminated. The cur- 
rent distribution is then used to compute the reradiated 
field. The aperture-field method used by the others 
assumes that the reflector may be replaced by the 
aperture formed by its projection in the principal direc- 
tion of reflection. The aperture is then considered to be 
illuminated by the field obtained by simple geometrical 
reflection from the reflector. 

The concept of gain is introduced as a means of evalu- 
ating the reflector performance. Unger and Crosby de- 
fined the gain as the ratio of the distant field intensities 
obtained with and without the reflector in the system 
when serving as a transmitting antenna. Jakes con- 
sidered a receiving system and defined the gain (in- 
directly) as the ratio of the received signal powers 
obtained in the two cases. The reciprocal properties of 
antenna systems make these definitions equivalent. The 
analysis to follow uses the aperture-field method and 
assumes a transmitting system. 


MATHEMATICAL ANALYSIS 


The curved reflector is assumed to be a segment of a 
very large paraboloid of revolution whose focus lies at 
the center of the illuminating dish throat and whose 
axis is directed toward the distant receiving point. A 
typical system is illustrated in Fig. 1 from which it is 
seen that the practical reflector is only slightly curved 
and resembles a flat ellipse mounted directly above the 
dish at a 45-degree angle. 

The aperture system is shown in Fig. 2(a). The feed 
aperture is taken as the dish throat which is circular of 

6S, Silver, “Microwave Antenna Theory and Design,” Rad. Lab. 


Series, vol. 12, McGraw-Hill Book Co., Inc., New Yorks N\ieve, 
chaps. 3-6; 1949. 
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radius P. It is separated by a distance h, corresponding 
roughly to the tower height, from the reflector aperture 
which is taken as the projection (assumed circular of 
radius R) in the direction of the path. While the re- 


REFLECTOR 


/ h 
gees 


/ PARABOLA 


h/2 


DISH 
Fig. 1—Typical passive-reflector antenna system. 


flector aperture itself is flat it appears curved to the sig- 
nal rays because of the reflector curvature. Since the 
reflector is a segment of a paraboloid it is seen that its 
aperture is a segment of a sphere of radius h. The far- 
field geometry is shown in Fig. 2(b), where the true re- 
flector aperture is shown in relationship to an arbitrary 
» point in space. 

By the Kirchoff-Huygens approximation, the diffrac- 
tion field E(x, y, z) at an arbitrary point in space due to 
radiation from a surface S on which the field is de- 
scribed by E’(é, n, ©) where (£, n, ¢) are coordinates of 
the surface is given by 

1 gt ar 
R(x,92)=j—f BGnt) as, e=>- 
Ads l r 

where / is the distance from the point (&, 7, ¢) to the 
point (x, y, z), \ is the wavelength and @ is the phase 
constant. In general, both Z and E’ are complex quan- 
tities. The approximations involved are that the dimen- 
sions and the radii of curvature of the reflector are large 
| compared to a wavelength and that the point (x, y, 2) 
does not depart greatly from the axis. 

The dish illumination is assumed herein to be given 
by the parabolic function 


Ep(a) = L— kx?, (2) 


where # is an adjustable parameter and x is the normal- 
ized radius variable. When k=O, the illumination is 
uniform, when k=1.0, it is completely tapered to zero 
at the edge. The case k = 0.684 gives the common 10 db 
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taper (i.e., 10 db down at the edge). The phase at the 
dish aperture is assumed constant. 

The distance / between points on the apertures is well 
approximated by 


p” — 2rp cos 
2h 


l a h + ) (3) 
where / is the average dish-to-reflector spacing, w is the 
angular spacing between the arbitrary points and p and 
r are the dish and reflector variables, respectively. Tak- 
ing dish aperture as surface and (3) for / leads, upon sub- 
stitution into (1), to reflector-aperture illumination 


Bal AP D8 nay : 
Erx(r) = j ee it Ep(p/p)Jo (=) ep" l2hodp, (A) 


0 


REFLECTOR APERTURE 


DISH APERTURE 


(a) DISH REFLECTOR GEOMETRY 


4 
(4,0, &) 
EER in FIELD 
& POINT 
ae 


REFLECTOR APERTURE 


(b) REFLECTOR- FAR-FIELD GEOMETRY 


Fig. 2—Aperture configuration for curved-reflector system. 


where Fp is the dish illumination. Making the changes 
of variables 
p r BP? R 


0 Sy ean ) Of) caer 


5 
P R 2 h 1 (8) 


and using the dish illumination of (2) the reflector il- 
lumination becomes, upon integration, 


En(y) = je rei? f(y), (6) 
where 
k 
f(y) = [U1 + jU2] — a8 [Us + jUs] 
+ kq?y?[Us + jU4] (7) 
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Fig. 3—Gain characteristics of curved reflector. 


and U,=U,(p,pqy) are Lommel’s functions’ of two 
variables defined by 


co w n+2m™ 
U,,(w, z) a 3 (ee 1) i (=) Jerome): (8) 


m=0 3 


The radiated field at a distant point is found by ap- 
plying the procedure again using the reflector aperture 
as the radiating surface and (6) for the illumination. 
Referring to Fig. 2(b), the distance / from a point (7, 8) 
on the reflector aperture to the point (d, ¢, a) is seen to 
be approximated by 


l=d—r cosy sin a, (9) 


where d is the distance from the center of the reflector 
aperture, y is the angular spacing between the points 
and a is the angle off the axis. Performing these substi- 
tutions and integrating leads to the distant field 


BR? 1 
B= ——einrroe ion fH) (by)ydy, (10) 
0 
where 
b = BRsina. (11) 


The distant on-axis field which would appear if the 
dish alone were used is obtained simply as 


EI ey or ( _ *) (12) 
E—) —€ 4 ; 
‘ed 4 


while the distant on-axis field from the reflector is ob- 
tained from (10) by setting a to zero. The power gain is 
then the square of the magnitude of the ratio of these 
fields and is given by 

2 


(13) 


Rd 4q° 1 
e an 3 “ip f(y) ydy 


7G. N. Watson, “A Treatise on the Theory of Bessel Functions,” 
2nd ed., Cambridge Univ. Press, London, England, pp. 537-550; 
1982. 
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Fig. 4—Total gain characteristics of curved and 
flat reflector systems. 


The radiation pattern of the reflector is the square of 
the magnitude of the ratio of the far-field E as given by 
(10) to its on-axis value. Thus 


9 


f aM coe 


f tow | 


The solutions to the various integrals appearing in (13) 
and (14) are summarized in the Appendix. 


ve (14) 


GAIN CHARACTERISTICS 


The gain characteristics of the curved reflector are 
presented in Fig. 3 as a function of 4/BRP for various 
values of R/P, the reflector-dish radii ratio. The dish 
illumination has been taken as 10 db-tapered to con- 
form with illuminations actually observed. Each curve 
has a maximum corresponding to the point representing 
the best compromise between the amount of power in- 
tercepted by the reflector and the uniformity of the 
intercepted illumination. Also, the curves rise without 
limit as R/P is increased reflecting the focusing ability 
of the curved reflector. 

Of more practical interest to the systems engineer 
than the gain of the reflector over that of the dish is the 
total gain of the tower, i.e., the gain of the reflector-dish 
combination. The dish gain over an isotropic radiator 
can be expressed as 


Gp = G(BP)? = GBh(PVB/h)’, (15) 


where G is the dish efficiency. If the reflector gain is ex- 
pressed as a function of P./8/h then the total gain in 
decibels can be written 


Gr — 10 log GBh = Gr(PV/B/h) + 20 log P»/B/h. (16) 


Fig. 4 is a plot of total gain with R\/B/h as a parameter. 
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Since the operating frequency is normally specified and 
the tower height is determined on the basis of providing 
proper terrain clearance, it is apparent that the curves 
amount to a plot of total gain per tower versus dish size 
for various reflector sizes. The optimum reflector-dish 
combination may be readily selected from Fig. 4. The 
smaller reflectors are seen to be relatively insensitive to 
dish size while the larger reflectors yield significant in- 
creases in total gain and require but small dishes. 

The total gain characteristics of the flat reflector were 
computed from Jakes’ curves and are shown on Fig. 4 
for comparison. Though not shown, the curve for 
Ry/8/h=1.0 probably falls only slightly below the cor- 
responding curve for the curved reflector while curves 
for larger values of Ry/8/h fall below that for 2.5. The 
curved-reflector system is clearly at an advantage when 
large reflectors are used with small dishes. 


1 ce 
Uso oe 
mee NG 
UNG a 
Ne 
Gis alien 


ate 

= +4 
< 
oO 

fe) 

-1 

=e. 

=o 

fe) 0.2 0.4 0.6 0.8 1.0 ik 
h/BRP 


Fig. 5—Effect of dish illumination on reflector gain. 


The effect of dish illumination is illustrated in Fig. 5 
where the reflector gain has been plotted for a fixed 
value of R/P and the values of k corresponding to com- 
plete taper, 10 db taper and uniform dish illuminations. 
These curves show optimum dish and reflector gains to 
be incompatible since optimizing dish gain by making 
its illumination uniform results in a sharper feed beam 
which causes the reflector illumination and gain to de- 
teriorate. In the vicinity of the optimum, the variations 
in dish gain are comparable to the compensating varia- 
tions in reflector gain making the total gain relatively 
insensitive to dish illumination. 


RADIATION PATTERNS 


Radiation patterns were computed for the curved re- 
flector for various values of BP?/h and are presented in 
Fig. 6 as a function of b, the angle parameter. The re- 
flector-dish radii ratio was chosen as 2 and the dish 
illumination assumed uniform for simplicity. It is seen 
that as the tower height is reduced the main lobe widens 
and the side lobes become depressed and move out- 
ward. At very low tower heights the main lobe begins 
to split and display a shallow minimum in the forward 
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direction. The effect of dish illumination is illustrated in 
Fig. 7 which is the radiation pattern at a fixed value of 
BP?/h and R/P. The principal effect of making the dish 
illumination more uniform is to reduce the side lobe 
levels. 

While radiation patterns for the flat reflector are not 
available in order to permit a direct comparison with the 
curved reflector it is possible to draw general conclusions 
from the fact that the principal difference between the 
two lies in the phase distribution of their aperture il- 
luminations. The curved reflector obtains its increased 
gain by presented equi-phase wavefronts to its aperture 
thus minimizing phase error. The flat reflector, on the 
other hand, produces a phase error at its aperture as it 
reflects the curved wavefront it encounters. Silver® 
points out that ihe effect of a quadratic phase error such 
as that displayed by the flat reflector is to increase the 
side-lobe level. 


db 


RELATIVE POWER 


BRSING 


Fig. 6—Typical curved-reflector radiation patterns. 
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Fig. 7—Effect of dish illumination on radiation pattern. 


At very large tower heights both reflectors may be 
considered uniformly illuminated thus giving identical 
patterns similar to that for p=0 in Fig. 6. As tower 
height is initially reduced neither reflector experiences 
much illumination taper but the flat reflector has an 


8 Silver, op. cit., p. 189. 


172 


appreciable phase error even at large tower heights so 
that its side lobes increase in amplitude. At lower tower 
heights the illumination taper probably predominates 
over the phase error thereby permitting the side lobes to 
recede. The radiation patterns of Fig. 6 show that the 
side lobe levels for the curved reflector decrease con- 
tinually as tower height is reduced, the decrease being 
due principally to the tapering of its illumination since 
phase errors are corrected except at very low tower 
heights. It may be assumed from the foregoing that the 
curved reflector will generally have a lower side-lobe 
structure than the flat reflector in the range of practical 
tower heights. 


CONCLUSION 


The curved reflector offers increased gain as com- 
pared with the flat reflector under all practical condi- 
tions; values of from 2 to 4 db more gain per tower 
being typical. Furthermore, for a given reflector the 
optimum dish size is smaller enabling a reduction in 
cost. Dish illumination is relatively unimportant with 
respect to the total antenna gain but the side-lobe levels 
are reduced when the dish illumination is made more 
uniform. In general, the curved reflector has lower side- 
lobe levels than the flat reflector thereby reducing stray 
radiation and providing greater immunity to interfer- 
ence. 

Practical curved reflectors are formed by indenting 
flat reflectors at several points; the most modern types 
being continually adjustable. The value of curving has 
been demonstrated repeatedly in the field by erecting 
such reflectors when flat, aligning the system for maxi- 
mum gain and then adjusting the curvature for a new 
maximum. The system gain invariably increases by this 
procedure. Rough gain measurements on systems which 
employ curved reflectors are in good agreement with the 
theoretical results thus tending to confirm the approxi- 
mations made in the analysis. However, actual reflectors 
exhibit greater beam splitting and side lobe formation at 
very lower tower heights than predicted herein. It is felt 
this may be attributed to the nonperfect curvature. 


APPENDIX 


The integrals in (13) and (14) are all represented by 
the following general forms: 


1 
(a) f 0.06, pay) yay v= 1, 2,3 
0 
1 
(b) f Up, pay) y'dy y= 3,4 
ze (17) 
() if Pe PhguaTulenlyoweay anit ied 
0 
1 
(a) f U,(p, pay)Folby)y%dy » = 3, 4 
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Type (a) is evaluated as 


t 1 
f U,(p, pay)ydy = Fu [U,-1(p, 0) — Uy-1(2, 99) |, (18) 


where 
pete?) 
=="Si =) 
2 


Uo(?, 0) = cos, Ui(p, 0) 


Ux(p, 0) = 1— cos : (19) 
Type (b) is given by 


1 1 
f U,(p, pay) y"dy = — — {Olt 9) 
0 PT 


2 
ae me [U,2(p, pg) — Us-2(f, 0) i - (20) 


Types (c) and (d) are cases of the general form 


1 
f U.(p, pay) Inlby) y™*1dy (21) 
0 
since 
1 
i Use-pay tuo w ae 
0 
De 1 
-— f Ul plpoy\Tibylydy 
0 
zs f Eee ear eee (22) 
0 


Integrating by parts gives 


il 
if Ualb, pavFalby lye tdy 
0 


ee) Pg? ™ ; 
> S) Urem(Ps baInzzzm(b) (23) 


which, for convenience of computation, may also be ex- 
pressed as 


1 
f U.(d, pay)Ir(by) y"tdy 
0 


bn Pe m 
Ga iwauce lp Ace 
m=0 enon 2 ) Uvim(P, PQ) Ant14-m(b) (24) 
where 
2°! 
A, (x) = . Seley (25) 


For the case of uniform illumination (k=0) these 
integrals become 


1 1 b 
f soysay = <4 | cos = SRT po) 
0 2g 2 
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At an infinite tower height (p=0) the limiting form of 
the radiation pattern is simply 


pas (=). 
b 
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A Multiple Telemetering Antenna System for 
Supersonic Aircraft* 


R. E. ANDERSON,} C. J. DORRENBACHER,{ R. KRAUSZ,t 
AND D. L. MARGERUM}+ 


Summary—tThe existing angle-of-attack indicator spike and ram- 
pressure tubes on a high-speed aircraft are modified to permit their 
use as telemetering antenna elements. The antenna system utilizing 
these elements provides for operation of three airborne telemetering 
transmitters. Operation of each transmitter and its antenna is inde- 
pendent of damage to or malfunction of the other components in the 
system, and no multiplexing loss is suffered in the antenna system. 
Since current flow on the fuselage is limited to the vicinity of the 
antenna elements, antenna performance would be little affected by 
variations in configuration of aircraft with similar basic design, and 
direct adaptation of this system to other aircraft of this type, there- 

- fore, is relatively straightforward. 


INTRODUCTION 
Gece OBJECTIVE of the work described in this 


paper was the development of an antenna system 

suitable for the transmission of telemetered in- 
telligence from supersonic aircraft. The antennas de- 
veloped provided for the radiation of energy from three 
independent airborne transmitters operating in the 216 
to 235 megacycle telemetering band. 

The antenna system discussed here utilizes existing 
ram-pressure tubes and the Angle-of-Attack Indicator, 
which were modified to permit use for the desired an- 
tenna functions. A pair of dummy ram-pressure tubes 
were added to serve as the antenna for the third trans- 
mitter. 

ANGLE-OF-ATTACK INDICATOR ANTENNA 


The Angle-of-Attack Indicator, or spike, normally 
protruding from the nose of a supersonic test vehicle 


* Original manuscript received by the PGAP, November 10, 1954; 
revised manuscript received, April 1, 1955. ; 

+ Stanford Res. Inst., Mount Lee Lab., Los Angeles, Calit. 

t Douglas Aircraft Co., Inc., Santa Monica, Calif. 


offers obvious possibilities for antenna purposes. Since 
this protrusion is about \/4 in length the impedance 
match is easily obtained over the telemetering band. 
The problems in adapting such an installation are: 
(1) excitation of the spike as an antenna with provisions 
for carrying the wiring from the Angle-of-Attack Indi- 
cator across the feed gap, and (2) realization of a satis- 
factory radiation pattern in the forward direction. 
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Fig. 1—Angle-of-Attack Indicator antenna installation. 


Feedpoints for the Angle-of-Attack Indicator antenna 
were provided by replacing a short section of the nose 
casting with a Fiberglas structure located about a 
quarter-wavelength from the end of the spike, as shown 
in Fig. 1. The spike is connected to the center conductor 
of the coaxial feed cable through a shunt matching ca- 
pacitor, the outer conductor being bonded to the fuse- 
lage. The Angle-of-Attack Indicator wiring is brought 
across the isolated section through the center conductor 


174 


which is isolated by a coaxial quarter-wave choke and, 
therefore, does not disturb the feedpoint impedance. 
Leakage of the signal into the fuselage interior is also 
prevented. 

The radiation pattern of a stub in the nose of a thin 
aircraft fuselage longer than a quarter-wavelength be- 
comes directive in the aft direction. Since this type of 
pattern has insufficient coverage forward for most tele- 
metering applications, steps were taken to prevent cur- 
rents from flowing on the fuselage aft of a quarter- 
wavelength from the feed point. This was accomplished 
by incorporating a coaxial quarter-wave choke into the 
nose section of the aircraft as shown in Fig. 1. Attempts 
were made to use a choke whose inner conductor was 
only slightly smaller than the skin of the nose section; 
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FUSELAGE LENGTH=6A 
Fig. 2—Radiation patterns of Angle-of-Attack Indicator antenna 
for circularly polarized receiving antenna. 


however, it was found that such configurations were in- 
effective in preventing current flow on the remainder of 
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Fig. 3—Impedance plot of tuned AAI antenna. 
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the fuselage. The lack of effectiveness was attributed to 
the large fringe capacity across the input to the choke 
and the low characteristic impedance of the coaxial sec- 
tion, which, in the presence of even small losses, will 
_result in a lowered impedance at the input to the choke. 
By reducing the center conductor of the choke section 
so that the characteristic impedance was about 100 
ohms, the desired performance was obtained. The re- 
sults of this choking action can be seen in the radiation 
patterns shown in Fig. 2, on the previous page. 

The tuning capacitor consisted of a brass plate (con- 
necting the center conductor of the feed line to the 
spike) spaced from the fuselage bulkhead by a Teflon 
gasket (see Fig. 1). The excursion of impedance with 
frequency about the design frequency after tuning is 
shown in Fig. 3, on the previous page. 


LONGITUDINAL 


Fig. 4—Radiation patterns of two opposite out-of-phase ram- 
pressure-tube antenna elements for circularly polarized receiving 
antenna. 


RAM-PRESSURE- TUBE ANTENNA 


Other protruding elements offering possibilities for 
adaptation to antenna functions on supersonic test 
vehicles are the ram-pressure pickup tubes. The ram- 
pressure-tube installation on the airframe considered 
consisted of two L-shaped stainless steel pickups located 
180 degrees apart on the nose section of the fuselage. 
These elements were increased in size to permit imped- 
ance matching, and additional dummy tubes were added 
between the existing two to make a total of four ele- 
ments spaced 90 degrees apart about the periphery of 
the nose. These four elements were to act as antennas 
for two independent telemetering transmitters. 

Exciting opposite antenna elements out of phase can- 
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cels the longitudinal currents on the fuselage; conse- 
quently, each element operates as a bent stub in a 
ground plane rather than as a feed for exciting long- 
wire currents on the fuselage. Choosing the relative 
lengths of the transverse and longitudinal portions of 
these stubs allowed amplitude control of the relative 
radiation in the Es and £, polarizations. Fig. 4 shows a 
longitudinal pattern which is entirely satisfactory; the 
transverse pattern, however, does have two deep nulls 
with this antenna arrangement. It should be pointed out 
that both £, and Ey components of the field are present 
in these directions, and to ground stations antennas 
having polarizations other than circular, the maximum 
reduction in signal may occur at slightly different as- 
pects of the aircraft with the depths of the nulls also 
modified. 


GENERAL DISCUSSION 


Inherent in the system is the isolation between pairs 
of antennas, as can be seen from consideration of the 
diagram in Fig. 4. The radial currents on one pair are 
perpendicular to the elements in the opposite pair and 
thus induce no currents, and the longitudinal currents 
in one pair are phased so that no net current is induced 
in the longitudinal portions of the elements in the oppo- 
site pair. Measured decoupling between any two of the 
three telemetering antenna systems was in excess of 
33 db. ; 

A drawing of the ram-pressure-tube antenna installa- 
tion is shown in Fig. 5. Impedance matching of these 
antennas was accomplished by choosing the distance 
from the external antenna feedpoint to the feedline in- 
put and controlling the length of a coaxial shunt stub 
built into the antenna mount. A Smith Chart plot of the 
impedance of a pair of ram-pressure-tube antennas is 
shown in Fig. 6 (next page). 
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Fig. 5—Ram-pressure-tube antenna installation. 
CONCLUSION 


The antenna system described in this paper provides 
for transmission from three airborne transmitters in the 
217 to 235 megacycle frequency range. All of the an- 
tennas are located in the nose section of the aircraft. The 
electrical performance of the antenna system, although 
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Fig. 6—Impedance plot for ram-pressure-tube antennas. 


somewhat dependent on the diameter of the fuselage, is 
essentially unaffected by its length and shape. Since 
each transmitter operates with separate antennas there 
are no losses incurred (such as are commonly encoun- 
tered in duplexing networks), and damage or malfunc- 
tion of one transmitter or its antenna does not affect the 


operation of the remainder of the system. 

The mechanical design of the antennas permits their 
installation on supersonic vehicles without appreciable 
increase in drag or weight, and since the tuning devices 


are built into the assembled antenna units, installation 
problems are minimized. 
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Measurement of Electric Field Distributions* 
R. JUSTICE} anp V. H. RUMSEY? 


; Summary—The usual method of measuring electric field dis- 
tributions comprises a receiving antenna at the point of observation 
connected by means of a transmission line to receiving apparatus a 
some distant point. The principal advantage of the method described 
here is that it eliminates the transmission line and thus makes possi- 
ble many measurements which cannot be done by the usual method. 
It is based on the fact that the echo from a thin straight conducting 
wire is proportional to the square of the component of electric field 
parallel to the wire. The practical application at centimeter wave- 
lengths is described. The method has proven most valuable in a 
variety of problems involving flush mounted directional antennas, 
radomes, artificial transmission lines, etc. 


INTRODUCTION 
IF on prot Principle applied to antenna radia- 


tion problems states that the radiation pattern 

of any source is determined by the distribution 
of the tangential components of electric or magnetic 
field intensity on any surface which encloses the 
source.!:? This well-known result is particularly useful in 
problems concerning radiating apertures in perfectly 
conducting closed surfaces such as a slot antenna 
mounted in an infinite cylinder or in an infinite ground 
plane, or more generally, a mathematically closed sur- 
face. In such instances the radiation patterns are de- 
termined by tangential E& in the aperture. The consider- 
able interest in such antennas shows the importance of a 
system for measuring accurately electric field distribu- 
tions in the vicinity of such sources of radiation. Con- 
ventional probing techniques are frequently impractical 


_ since they introduce more disturbing elements into the 


field than can be tolerated. In particular the transmis- 
sion line from the probing antenna to the receiving ap- 
paratus may so disturb the field as to invalidate the 
measurement. Furthermore, in making phase measure- 
ments, it is necessary to transmit the rf signal picked up 
by the moving probe back to the fixed receiving ap- 
paratus, and thus involve passing the rf through a 
flexible transmission line (or waveguide). The distortion 
of the rf signal caused by flexing the transmission line is 
a further objection to the usual probing method which 
is serious at short wavelengths (e.g., 3 cm). The purpose 
of this paper is to describe a reflection method of meas- 
uring electric field distributions in the vicinity of the 
source: its principal advantage is the elimination of the 


* Original manuscript received by the PGAP, December 20, 1954; 
revised manuscript received, April 1, 1955. 

+ Ohio State Univ., Columbus, Ohio. : a, 

t University of Illinois, Urbana, III.; formerly at Ohio State Uni- 
versity. 

u eisaen London Math. Soc. Proceedings, vol. 1, p. 1; 1904. 

2S. A. Schelkunoff, Bell Sys. Tech. Jour., vol. 15, p. 92; 1936. 


transmission line from the receiving antenna to the re- 
ceiving apparatus. 


PRINCIPLE OF THE METHOD 


Consider the measuring system sketched in Fig. 1. 
An antenna, represented by a horn, is connected to a 
transmitter through a hybrid junction. The junction is 
balanced so that when the antenna is radiating in its 
normal environment no energy is detected in the receiv- 
ing arm of the junction. If now a conducting scatterer is 
introduced into the field of the antenna, a signal will be 
observed at the receiving arm of the hybrid junction. 


Test 
Antenna 


Matching 
Stubs 


Transmitter 


Detector 


Logarithmic 
Recorder : 


Fig. 1—Schematic of the measuring system. 


Dummy Load 


The principle of the method is to use this signal, denoted 
by V, asa measure of the electric field, Z, at the position 
of the scatterer. It follows from the reciprocity theo- 
rem? that V and E are related by an equation of the 


form 
V=C f f E-JdS, 
scatterer 


where C is a constant, E=field when the scatterer is 
absent, J=current distribution on the scatterer. If the 
scatterer is a very short length of thin straight wire, 
(1) reduces to 


(1) 


V = CEllil, (2) 


where E;=component of £ in the direction of the wire, 
and JJ/=current moment of the current on the wire. The 
value of Il is determined from the condition that the 
total tangential electric field at the scatterer must 
vanish, which implies that 


(3) 


3S. A. Schelkunoff, “Electromagnetic Waves,” D. Van Nostrand 
Co., Inc., New York, N. Y., p. 477; 1943. 

4V. H. Rumsey, “Reaction concept in electromagnetic theory,” 
Phys. Rev., vol. 94, p. 2; September, 1954. 


E\(J) = — (tangential E; at the scatterer), 
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where E(J) =field that would be generated by the cur- 
rent distribution J. If the scatterer is very small com- 
pared with the wavelength, E(J) at the scatterer wilt be 
very nearly the same as if the scatterer were in “free- 
space.” In this case we can write 


E,J) = K'tl 


and therefore, from (3), 
Vi IG Dp (4) 


where K is a constant. Substitution in (2) for Z/ from (4) 
gives | 
oS ABE. (5) 


where A is a constant, provided the scatterer stays in 
the same medium (same yp, € and a). Eq. (5) shows that 
the echo voltage V is proportional to the square of the 
component £; of electric field and the phase angle of V 
is twice the relative phase angle of E;. We note that E 
represents the field when the scatterer is absent; i.e., E is 
the thing we are trying to measure. This result is the 
basis of the method of measurement. If the scatterer is 
moved through the field to be measured and the echo 
signal is recorded as a function of position of the scat- 
terer, the result is a continuous measurement of the dis- 
tribution of the component of electric field parallel to 
the scatterer along the path followed by the scatterer 
(again note that the electric field referred to is the field 
that exists when the scatterer is absent). 

We note that the validity of the method depends on 
the scatterer being sufficiently small and thin to make 
(4) a reliable approximation. The approximation will 
begin to fail when the scatterer comes close to some re- 
flecting surface, for then we cannot assume that Hi( J), at 
points on the scatterer, can be calculated as if the scat- 
terer were in “free space.” To estimate the magnitude of 
this effect, calculations were made for a scatterer in the 
presence of a plane perfectly conducting reflector. (The 
calculations are based on an assumed distribution of 
current on the scatterer, and are therefore approximate, 
but they should be quite satisfactory for estimating the 
order of magnitude of the effect). The results are plotted 
in terms of the error in E due to assuming that (5) is 
correct (see Fig. 2). It will be seen that the error is less 
than one per cent provided the distance from the re- 
flecting surface is not less than the length of the scat- 
terer. It should be noted that Fig. 2 applies to a plane 
perfect reflector; presumably the effect will be more or 
less depending on whether the reflector is concave or 
convex. In an enclosed region such as a cavity or wave- 
guide, the effect will presumably be significantly greater 
or, to put it another way, the scatterer will have to be 
significantly smaller, as one would expect. 

We should note that the error due to the size of the 
scatterer is characteristic of any method of measure- 
ment. For example, the magnitude of the error in the 
standard field probing method is of the same order, if 


IRE: TRANSACTIONS—ANTENNAS AND PROPAGATION 


db 


October 


the probe consists of a dipole of total length which is 
connected to a detector of extremely high impedance 
(the error being due simply to the finite size of the 
antenna). 

2% 


Ground 
Plane 


1% 


sine 
ie Scatterer % Error 
it 
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Fig. 2—Scatterer near ground plane, coordinate system and 
error vs distance from ground plane. 


An experimental verification of (5) was obtained by 
measuring the transverse variation of the electric field 
in a rectangular waveguide operating in the lowest order 
TE mode. The excellent agreement between the cosine 
distribution of the true field and the measured result, as 
shown in Fig. 3, confirms the prediction that conducting 
surfaces in the vicinity of the scatterer do not appreci- 
ably alter the square law dependence of echo voltage on 
incident electric field. 
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Fig. 3—Measured and calculated fields ina T Eon 
rectangular waveguide. 


INSTRUMENTATION 


A measuring system of the type described above has 
been in use at the Ohio State University Antenna Labo- 
ratory for some time. Some of the problems encoun- 
tered in developing the system and their solutions are 
now discussed. An isolation between transmitter and 
receiver of about 90 db is usually required. The need for 
such a fine balance of the hybrid junction requires that 
the transmitter be as nearly monochromatic as possible. 
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Just how near to monochromatic it has to be depends 
on the antenna under test. If it is well-matched and 
broadband, a good balance can be obtained even if the 
signal has a certain frequency spread; if it is highly se- 
lective, an extremely pure signal is required. It has been 
found that a reflex klystron immersed in a large undis- 
turbed oil bath generates a signal of sufficient purity and 
stability® for the requirements of the measuring system. 
The success of this device depends on the thermal ca- 
pacity of the oil being large enough to hold the trans- 
mitter temperature constant for the duration of a meas- 
urement. The necessary tuning to achieve the required 
degree of isolation between the transmitting and receiv- 
ing arms of the junction can be obtained from carefully 
designed and manufactured triple stub tuners, provided 
the test antenna is well matched into a waveguide. 

In both the original 3,000 mc system, and the 10,000 
mc system installed later, it was found that crystal de- 
tection was satisfactory. The energy was transmitted 
unmodulated and then modulated at the receiver by 
means of a crystal. This scheme was used because it 
furnished a more nearly monochromatic signal than that 
obtained by modulating the beam voltage of the kly- 
stron. The output of the detector was then fed into a 
selective audio amplifier whose output in turn was used 

‘to drive a logarithmic recorder. 

A problem of some importance in the actual experi- 
ments is that of supporting the scatterer and moving it 
through the field. It was found that the scattered field 
due to a fine nylon thread was completely negligible 
compared to that due to a conducting scatterer. Hence 
it was possible to attach the scatterer, usually a fine 
wire, to the nylon thread which was then pulled through 
the field. Synchronization of the motion of the scatterer 

‘and of the recorder then made it a simple matter to in- 
terpret the recorded echo voltage in terms of the electric 
field distribution. 

Another problem of considerable importance involves 
the ability of the scatterer to discriminate between 
mutually perpendicular field components; i.e., how thin 
must the scatterer be to select only the component £7? 
Measurements indicate that the length to diameter 
ratio must be about 30 to reduce the error to less than 
10 per cent of the cross polarized component. 
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APPLICATIONS 


The 3,000 me system was assembled in order to facili- 
tate the study of traveling-wave slot antennas. In this 
application many measurements were made of the dis- 
tribution of tangential electric field in the aperture of 
such antennas. A typical antenna of the type investi- 


gated is illustrated in Fig. 4 together with its coordinate 
system. 


Ground Plane 


Fig. 4—Slot antenna in ground plane with 
coordinate system. 


Since the phase of the signal V varies as twice the 
phase of the component £, it is possible to use this sys- 
tem to make phase measurements as well as amplitude 
measurements. The simplest way to obtain phase in- 
formation is to unbalance the hybrid junction thus pro- 
viding a constant reference signal R and then to record 
the sum of the constant reference signal and the echo 
signal produced by the scatterer. Thus for each position 
of the scatterer three voltages are known; namely V, 
V+R and R where R is the reference signal. Since R is 
constant for a given measurement the calculation of the 
relative phase distribution for that measurement follows 
directly from the application of the law of cosines to a 
series of distinct positions of the scatterer. Alternatively, 
it is possible to measure the relative phase distribution 
along the path of the scatterer by directly measuring the 
phase variation of the signal V. This latter method has 
been used successfully at 10,000 mc.* The first method is 
well suited for the measurement of the phase velocity of 
a guided traveling wave field. A typical result showing 
the sum of the echo and a fixed signal is shown in Fig. 5. 


aaed 


16 we — TaN? ‘ 3 
- yy \ rN - x zi 
AINA fy a l ++ 4 4 a tox x ee 
_ 14h Ni [po T a TaN ia Sn Shera 
Pikes QA v1 \ Ly Mert oylkl Vie ay Hi ee eee WO BP dA Le 
i sm ea | Ly 1 mall Ty) it Wa i 
LEVEL 4 ; - \ a 1! \! tt 4 ' | ‘ee Y vo = = 4 
ao) aah 1 “tt nts ti il WT f q a a 
i yo wt) eaee Waa 

Wee ao id : : 
6 + 4 ¢ zl z 
e j 
| —— |E, |b) ] [i cl = 2 
<| fe — —— |e, + REFERENCE | it cs re 
2 1 

oe jis -= th =, cf aa 5. 6a cat Bi EN 10x m0 ma re 13. 

Fig. 5—Amplitude and phase measurement of normal component of electric field in slot antenna. 
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The minima in the dashed line occur where the echo is 
180 degrees out-of-phase with the constant reference 
signal and therefore the phase difference between echoes 
at adjacent minima is 360 degrees. Thus the phase dif- 
ference between the values of £; at adjacent minima is 
180 degrees. The measurements plotted in Fig. 5 (on the 
previous page) show that the electric field distribution 
under observation can be approximated by a traveling 
wave of slightly variable wavelength (or phase velocity), 
the wavelength being equal to twice the distance be- 
tween consecutive minima. 
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Other applications of the method include the study of 
the effect of radomes on the performance of antennas 
and antenna systems, and the study of the propagation 
along artificial transmission lines. 
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Determining the Reflector Surface of a Radar 
Antenna with Point Source Feed* 
PENTTI LAASONEN}{ 


Summary—In the construction of a radar antenna the basic prob- 
lem is to give the antenna reflector such a design that the scattering 
pattern will conform to the required power distribution. An approxi- 
mative method is given for the determining of the central section of a 
reflector with point source feed and enveloping the osculating pa- 
raboloids, as well as a formula for the computation of the beam re- 
flected by this surface. The central section obtained may be improved 
successively until adequate accuracy is attained. The method of 
computation is illustrated by an example. 


mathematical treatment of a radar antenna are 

to determine the beams of an antenna of given 
shape; i.e., the power distribution in the scattering pat- 
tern, and to find such a design of the reflector that it will 
radiate a pattern of the desired power distribution in the 
various directions. 

Inasmuch as either problem is the inverse of the other 
as far as its starting point and its aim are concerned, and 
as problems of the first kind are more easily solved than 
those of the second kind, the normal procedure will be to 
find the solution of a problem of the second kind by 
successive approximations involving the recurrent solv- 
ing of problems of the first kind. This is why the solution 
formulas of the first problem will be derived as a first 
step in the following considerations, the objective of 
which is to describe the solution of the second problem 
in the case of a radar antenna with point-feed source. 

For the reflected pattern, the following expression will 
apply for the gain function of the pattern formed upon 
reflection, in the direction of the vector R®: 


ogee TWO basic problems encountered in the 


e Original manuscript received, November 12, 1954; revised man- 
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where X is the wavelength, a and i are two unit vectors 
orthogonal to the direction in question and to each 
other, and J is a complex vector defined through the sur- 
face integral 


am bebe: 
di f Garcia hear [(2 x 7) x da] 
a xs 


extended over the entire reflector area a. 


G;=gain function of the feed 
k=27r/X=propagation constant in the nonlossy 
media 
7=vector from feed point to reflector surface 
é° =unit vector defining the polarization of the radia- 
tion incident upon the reflector 
dai=n da=surface element of the reflector surface 


f= 1. 


We shall limit the following discussion to a reflector of 
such symmetry that the reflector surface and the pri- 
mary pattern G; of the antenna have a common plane of 
symmetry and that the direction of polarization of the 
fictive dipole idealizing the radiating waveguide is 
orthogonal to this plane. We now take the reference 
system of co-ordinates with the y-axis orthogonal to the 
said plane and with the unit vectors 7, 7 and k. We shall 
restrict our problem to determination of gain function G 
in directions lying in the xz-plane of symmetry only. 
Consequently, then, R° lies in this plane of symmetry 
and J is orthogonal to it, whence G obtains value 
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The polarization vector é° pertaining to the ray 7 in- 
cident upon the reflector lies in the plane determined by 
the direction of polarization j of the transmitting dipole 
and by the ray 7, and it is orthogonal to the latter; i.e., 
it is obviously the unit vector pertaining to the vector 

eee Te, 


For the y-component of the vector (¢X#) Xdd@ in the 
integrand in the expression of J there is now found the 
expression 


j \(@X 7) X da] = (j-7)(@-da) — (7-8) (#-da) 
= (j-A)(j-da) — #-da 
— (r.da, + r.da,). 


Since further the square of the length of @ is (r,2-+1r,2)/r, 
we thus have 


ll 


r,da, + r,da,|* 
MAT ante hee 


When it is attempted to concentrate the power in the 
direction of the xz central plane, it is expedient to shape 
the reflector to conform with the enveloping surface of 
such confocal revolution paraboloids pertaining to the 
directions in this plane that are of interest that the 
osculating paraboloid corresponding to any direction 0 
in this plane will be tangent to the envelope along a 
plane curve; i.e., a parabola, with this particular eleva- 
tion 6 (Fig. 1). If the distance of the intersection of this 
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a 


central section 


generating 
parabola 


contour 
line 


Fig. 1 


parabola and the central plane, Q, from the focal point 0 
is p and the elevation of the radius 0Q is ¢, so that 
26 =$+6 is twice the angle of reflection, a simple calcu- 
lation shows the focal distance of the parabola to be 
p costs. Let now S(é, 0, ¢) denote the focus of the 
parabola, P(x, y, z) an arbitrary point on this parabola, 
and at the same time also a point on the reflector sur- 
face, and let the azimuth angle of this point be PSIO=¢. 
The said points have the co-ordinates: 
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&= p sin @ + p cos? a sin 8 
¢ = — pcos ¢+ p cos? o cos A 
x= psing+ pcos’o tg? $9 sin 0 
ys 2p cos? ¢ tg 49 
2 = — pcos¢-+ p cos’ a tg? S¢@ cos 0. 


Using the parabolic co-ordinate a defined by the equa- 
tion 

tg a = cosa tg dy, 
instead of the azimuth angle ¢, the following expressions 


are found for the co-ordinates of a point P on the re- 
flector surface: 


x= psing+ p tg? ga sin 8, 
y= 2p tg ga cosa, 
\2 = — pcos¢+ p tg? ga cos 6. 


Differentiation with regard to the variable ¢ shall be 
indicated by a prime. On the basis of the shape of the 
central section (a=0) of the reflector the following dif- 
ferential equation between the variables p, ¢ and @ is 
found: 


dp 
SR! IE MEI Cele SE (1) 
do 
and furthermore the functional determinants 
0(4, Z 
oe bated bile (1 + 6 tg? 4a) cos $a sin $(¢ — 9), 
0(¢, a) cos® 4a 
O(Z, x 2 
oP y= sanbiltr Toh omg alton iy 
0(¢, a) cos? 4a 
CG, 2 
(sy) = + ‘ (1 + 6’ tg? 4a) cos $a cos 3(¢ — 8). 
0(¢, a) cos? 4a 
In accordance with the theory of surfaces there is 
O(y, 2 JGae (4%, vy). 
ia = [rt Miah, 2 Eada 
d(¢, a)  9(g, a) A(, &) 


The differential expression appearing in the expression 
of G is then 
r,da, + 7,da;z 
Viet 
1+ @ tg? 4a) ae dod 
= SS ed SSS Q. 
al he /1 + tg? a sin? o 
Let 00 be the angle of elevation of the particular direc- 
tion R°® in the xz-plane for which the value of the gain 
function shall be determined. The exponent in the inte- 
grand in the expression of the gain function G(@) is 


— kj(r — 7-R) 


4 
= — — feos? 4(6 + %) + ta? dor sin® 400 — 00] 
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whence the expression of the gain function is ultimately 
found to be 


iL f (Cyl, 3? (L48"tg Ba) 


COS a 


/1+tg? asin? o 
2 
» @— Amp |») [oos” (@+09) /2+tg” @/2-sin® (6-80) /2)idgda : 


G(6o) ce 


(2) 


When the central section of the reflector has a given 
shape as determined in polar co-ordinates by a function 
p=p(¢), 6 (and consequently also a) will be determined 
through (1) in terms of ¢. The gain function Gy; ex- 
pressing the power distribution of the radiation from 
the waveguide is originally determined through two 
directional co-ordinates; e.g., the spherical angular co- 
ordinates 


‘ = arctg (— Va? + 97/2), 
y = arctg («/¥). 

Since, however, these are related with the variables ¢ 
and a through the equations 


| cos 8B = cos ¢ cos? 4a — cos @ sin? 3a, 
sin 8 cos WY = sin a cos a, 


the function G;(¢, a) will thus be determined. All in all, 
the evaluation of (2) then involves the integration of a 
definite function of the variables ¢ and a; ¢ designating 
a certain parabola section of the reflector and a@ repre- 
senting the parameter which defines the location of a 
point upon this parabola. The domain of integration is 
limited by the contour line of the reflector represented 
by the boundary function a=a.(¢), whence —a,(¢) Sa 
<+a.(¢). 

The double integral in the expression of G(@o) has to 
be computed for each requisite value of 6) by means of 
some method of numerical integration. The dependence 
of the integrand of variable a is regular enough, both as 
regards the exponential function and the factors pre- 
ceding it; that, e.g., Gauss’ method of mean values can 
be employed with advantage in the integration with re- 
gard to a. The integrand being even with regard to a, it 
is adequate if, for instance, six intermediate points are 
used, to compute the integrand for the different values 
of @ at three a values chosen in accordance with the 
mean value method. The dependence of the integrand 
of variable ¢ is mostly irregular to such a degree that the 
integration with regard to this variable should be 
carried out in accordance with some more versatile 
method; e.g., Simpson’s rule, but increasing the number 
of subdivisions as occasion demands. 

The factor in the integrand most strongly dependent 
on the shape of the central section of the reflector is the 
exponential function; at various points on the reflector 
it obtains complex values of absolute value unity, but 
with phase angles which may vary within a range of 
several full turns. In order to examine more closely the 
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nature of this exponential function, we shall expand its 
phase angle to a Taylor series in the vicinity of the 
value #9 corresponding, in accordance with (1), to the 


| 
: 
: 
) 


direction 6) under consideration. It will be seen, then, | 


that the first-order term always disappears; the zero- 
and second-order terms are 


-2 it] = {- 
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Arpo 


COS” do 


= <* oo'(1 + 69’ tg? ga)(o — bo)? + °° ae 


It is evident from this that the said phase angle has 
a stationary value along the osculating parabola 
section @o belonging to the direction 6, its steepness, 
however, increasing with the increase of the factor 
(apo/d)O0'(1+60'tg?3a). Due to the integral formula 


“upe 1 
ie 
f e(atet’)ig¢ = ccereinig /=, 
ps b 


we thus arrive at the approximative formula 


Po 
Cpe ee ee 
(Ao) x] ay] Cos? ao 


1+ 0 tg? 3 
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This formula also constitutes the starting point for 
the finding of the first approximate solution to the in- 
verse problem; i.e., that of determining such a central 
section p=p(¢) of the reflector that the gain function of 
the reflected pattern will be a given function of the di- 
rectional angle 6. If, namely, this formula is approxi- 
mated further by substituting for its integrand the 
value it obtains at the center (w=0), there results, with 
an integration interval of the length 2a.(¢0), the ap- 
proximate formula 


G0) = K~ cost H+ NGG, Ya(4))%, (3) 
where K is a constant. 

In addition to this equation relating the variables p, 
¢ and 6 we have the differential equation (1) containing 
these variables. Eliminating p and the indefinite con- 
stant K from these equations, we have the differential 
equation of the second order: 


ae Ge Liproees ae db 

d¢? G;(¢, 0) a(p) / do 
re (= + tg3(¢ + ®) (=) - 0 
G() do ‘ 


which has a solution satisfying the given boundary con- 
ditions 0(¢1) =6:, 0(g2) =62; to this solution corresponds 
a certain shape of the central section: 
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p = poelte 1/2(4+0(4)) a6. 


This first approximate solution is, however, most 
easily obtained directly from (3) by integration of its 
both members with regard to 6: 


ie G(0)dé6 


$ 
=K J Gx(6, O)(ac(9))*0(@) cost Ho + 0C6))d6. 4) 
1 
It is true that the right member of this equation is not 
exclusively a function of ¢, since the integrand contains 
the unknown function 6(¢) immediately as well as 
through p, but a process of iteration starting with an 
assumed function is generally very rapidly convergent. 
The constant K is naturally determined in such a way 
that even the second boundary condition 6(¢:) =62 is 
satisfied. 

When now the power distribution of the pattern 
actually reflected by this first reflector shape deter- 
mined in the way described above is calculated from (2) 
it will not generally coincide with the required G(@), but 
it has a certain imperfect configuration G(@)-+AG(@). A 
new, improved shape of the reflector is now found by re- 
peating the calculations outlined above, inserting, how- 
ever, for the required pattern the function G(@) modified 
in the opposite direction; i.e., G(@) —AG(@). The pro- 
cedure usually leads quite rapidly to a reflector shape of 
such perfection that it will reproduce the required pat- 
tern with adequate accuracy. This convergence has been 
verified in practice. Moreover, a reflector based on these 
calculations has been constructed and measured. 


EXAMPLE 


As an illustration of the application of this method we 
shall consider a problem of schematic simplicity: to de- 
sign a reflector meeting the following requirements: 


1. The power distribution in the primary pattern 
radiated by the waveguide conforms to the func- 
tion G;=cos? 6. 

2. The reflector area is confined to the region with 
G;20.25. 

3. The dimensions of the reflector are determined by 
the condition that the upper edge of the central 
section should have the distance 10 from the feed 
point. 

4. The central section of the beam of this antenna is 
required to conform to Fig. 2; i.e., there should be 


G(6) sin? 6=constant, when arctg 0,150530°, 
G(@) =constant, when 004 <arctg 0,1, 
G(6) =0, when 6<0 or 0>30°. 


It was seen above that the solving of this problem is 
essentially connected with the determination of the 
function 6(¢) corresponding to the shape of the central 
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section. The boundary values of 6 being 0 and 30° and 
those of ¢, with g=0 and consequently a=0, being ob- 
viously —60° and +60°, the following boundary con- 
ditions may be set up for this monotonously changing 
function: 


when ¢= —60°, we must have 6 =30°, 
when ¢= + 60°, we must have 6=0. 


8-arctg 0,1 


Fig. 2 


As a first trial function we shall employ the linear 
A(p) = 15° — ad. 


According to (1), the shape of the central section corre- 
sponding to this function is 


p = polcos ($ + 7°5) |-8/8, 


In Table I, these 0 and p functions of the first approxi- 

mation are listed under A. The a, values pertaining to 

the reflector edge have been obtained on the basis of the - 
equation 


cos B = cos ¢ cos? 4a — cos 8 sin? 4a (5) 


substituting 8 = 60°, which corresponds to the edge. 


TABLE I 

A B G 

co 0 p/r Qe 0 p/x Oe 6 
—60° aN? fee 0® | 30200, 81455 0700 30°00 
—40° D5 eer Or SOeraicanie2 0) ol 7.833 46°56 20°12 
—20° ZOE tO, S20 Se 9°68 7.469 57°08 8°79 
0° {52 ee O99. 81008 Soli Pete COROU 4°57 
+20° 10° A Dae Deehily Tp Ss a elasters 2°02 
+40° 5° 8.44) 462 0°53 8.482 45°67 0°44 
60° Om 0 O0REN Oy 0°00 10.000 0°00 0°00 


The second approximation of the functions 6, p and a, 
will be found by means of (4). Since the integrand in the 
left member of this equation is known in terms of @, and 
the integrand in the right member can be computed by 
substituting 6, p and a, in terms of and in accordance 
with the first approximation, these integrals can be 
evaluated by means of numerical or graphical integra- 
tion in terms of 6 and @, respectively. Plotting the re- 
sults in a rectangular (¢, 0) system of co-ordinates with 
the values of the left-hand integral measured along the @ 
axis and those of the right-hand integral along the 
@ axis, the relation between ¢ and 6 can be found by 
means of a simple nomographic construction shown in 
Fig. 3. The columns joined under B in the table give the 
function @ obtained in this way as well as the functions 
p and a, calculated once more according to equations 


184 IRE TRANSACTIONS—ANTENNAS AND PROPAGATION October 


9 
{a(e)de 


207 


Tra 


-60° -50° -yo® -30° -20° -10° O° 10° 20° 30° oO? 50° 60° 


d 


Fig. 3 


(1) and (5). For purposes of control, a third approxima- 
tion of 6(@) has been calculated from these functions, 
using the nomograph in Fig. 4. The beam corresponding 
to this reflector shape has then been computed from (2) 
for a few 4 values. The double integral appearing here 
has been evaluated numerically, integrating first with 
regard to a, applying the mean value formula of Gauss 
on the interval (—a, +a.) at two abscissa values 
[+(1/4/3)a.|, then with regard to ¢, employing Simp- 
son’s rule on the interval (— 60°, + 60°). The result is 
shown graphically in Fig. 5 which gives the shape of the 
central section of the beam; i.e., the gain function G(@). 
The desired shape is also indicated in the figure to make 
a comparison possible. 


-60? -50° -4O° -30° -20° -10° ©° 10° 20° 30° 40° 50° 60° 


Fig. 4 


pane 


Fig. 5 


The next step would be to repeat the calculation, only 
with the difference of substituting for the desired func- 
tion G(@) such a fictive G which tends to correct for the 
errors observed ; thus, the value corresponding to 6 = 20°; 
e.g., should be smaller in an appropriate degree than the 
actual value. As this calculation would be quite similar 
to that carried out in the first step (except for being 
shorter, since there will not be any great changes in the 
reflector shape), our example has not been carried 
further. 


COS) 


| 
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Multipath Phase Errors in CW-FM Tracking Systems* 
T. E. SOLLENBERGERt 


; Summary—High-precision, cw, missile-tracking systems are sub- 
ject to multipath errors due to ground reflections and reflections 
from neighboring objects. To reduce the effects of ground reflection, 
advantage may be taken of the characteristic of frequency modula- 
tion to discriminate against a weaker interfering signal. This paper 
discusses the magnitude of the phase error reduction to be obtained 
by use of frequency modulation. 


It is demonstrated that the phase error can be reduced to any 
desired value if sufficient band width is utilized in the frequency 
modulation signal and significant improvement over an amplitude 
modulation system can be obtained with a wide-band FM system, 


INTRODUCTION 


[: RECENT years a number of phase-comparison 
tracking or position location systems have been de- 
veloped for guided missile testing. In such systems, 
ranging between the missile and a ground location is ac- 
complished by measuring the phase difference between a 
reference cw transmission and a replica of the reference 
transmission which has traveled to the missile and back 
to the ground location. 

Fig. 1 illustrates the principle of operation of a repre- 
sentative system.' A cw transmission consisting of a 
modulated carrier is transmitted to the missile where it 
is retransmitted, while maintaining coherence of the 
modulation. This retransmitted signal is received by the 
three ground receivers. The signal from the ground 
transmitter is transmitted by ground links to the three 
ground receivers and serves as a reference from which 
the phase delay in the modulation is measured. This 
phase delay is a measure of the sum of the ranges from 
ground to missile and back to ground. With the measure 
of three range sums and knowledge of the receiver 
positions relative to the ground transmitter, the position 
of the missile in space can be determined. 

Because the measurement of phase is ambiguous for 
every half-wavelength displacement of the remote trans- 
ponder on the missile, a series of sine wave modulations 
are placed on the carrier. These modulations usually 
have a decade relationship and provide coarse and fine 
range scales. The low-frequency modulation serves to 
resolve the ambiguity of the next higher frequency 
modulation. 

All such systems employing a phase comparison tech- 
nique for distance measurement are to some degree 
susceptible to errors caused by multipath reflections 


* Original manuscript received, October 21, 1954; revised manu- 
script received, April 14, 1955. This work was supported by Air 
Force Contract No. AF-08(169)-94. . Rae 

+ The Ralph M. Parsons Co., Electronics Division, Pasadena 2, 
ak. W. Johnson, E. I. Bosch, and P. H. Reedy, “Precision Auto- 
matic Ranging System—PARAN,” Rep. 955-RI, Phase I Tech. Rep., 
Part I, AF Armament Center, Eglin AFB, Fla. Contract AF 18(600)- 
202; November, 1953. (Unclassified.) 
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Fig. 1—Typical baseline trajectory-measuring system. 


from neighboring objects in the field of the antennas and 
particularly from the ground plane of the antennas. 
While these multipath errors may be considered syste- 
matic, they are sometimes difficult to recognize and iso- 
late since they are changing with time due to motion of 
the missile. Methods of reducing errors due to multipath 
consist of using highly directive antennas and, in the 
case of ground reflection, lowering the height of the 
antenna so that the path difference between the direct 
and reflected rays is reduced. The first method compli- 
cates the system since, with directive antennas, it is 
necessary to rotate the antenna to follow the missile. 
Reduction of antenna height often means loss of range 
coverage. Even if ground reflection could be minimized, 
there would remain problems caused by reflections from 
neighboring objects such as poles, trees, and buildings. 

Fig. 2 illustrates the order of magnitude of phase 
errors that may be encountered with a double-side-band 
amplitude-modulation system. These are computed 
curves for a reflection coefficient of 0.9 and a carrier-to- 
modulation frequency ratio of 8. The sharp phase-error 
peaks are due to interference at the carrier frequency. 
The envelope is a result of interference at the modula- 
tion frequency. 


2R. A. Epstein and P. H. Reedy, “Ground-Reflection Phase 
Errors of CW-AM Tracking Systems,” AF Tech. Rep. No. 7, 
AFMTC, Patrick AFB, Cocoa, Fla., Parsons Aerojet Co.; September, 
1951. (Unclassified.) 
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The thought that frequency modulation may have 
some immunity to multipath error is supported by the 
fact that frequency modulation tends to suppress the 
weaker of two signals in the same channel. It has been 
previously demonstrated that multipath interference on 
an FM channel generates a multitude of harmonic dis- 
tortion components. However, when the modulation 
consists of one or more discrete frequencies, as is the 
case in tracking systems, there is a possibility of filtering 
out most of the distortion components after detection. 


Fig. 2—Phase error in double-side-band AM for carrier-to-modulation 
frequency ratio of 8. Ordinate is modulation frequency phase error 
in radians. Abscissa is phase delay modulation of reflected signal 
in radians. 


CALCULATION OF PHASE ERROR 


. The problem considered here is that of determining 
the magnitude of phase error on an FM transmission 
after filtering of distortion. Actual computation of 
phase error has been carried out for a simple system 
which has been idealized in several particulars: 

1. To start with it is assumed that the transmission 
consists of a carrier which is frequency-modulated by a 
single sinusoidal wave form. Later, multitone modula- 
tion is also considered. 

2. A single interfering multipath reflection, such as 
from an infinite ground plane of the antenna, is assumed. 
The case of more than one interfering multipath reflec- 
tion is not analyzed because of the extremely intractable 
character of the mathematical manipulations. The case 
chosen for study is believed to represent a practical 
situation, since considerable care should normally be 
exercised to choose an antenna site having no nearby 
reflecting objects other than the ground plane. 

3. Both the amplitude and phase shift of the coeffi- 
cient of reflection are assumed to be independent of fre- 
quency over the bandwidth of the transmission. 

4, The receiver contains a perfect limiter preceding 
the discriminator. 

5. Perfect filtering of the modulations takes place 
after detection so that only the component of distortion 
at the modulation frequency is retained in the output. 

This assumption of a perfect filter is made to facilitate 
the analysis. In any real system the post-detection band- 
width must be consistent with the information con- 


IRE TRANSACTIONS—ANTENNAS AND PROPAGATION 


October 


tained in the measurement of phase. This bandwidth is 
the same as the bandwidth of the doppler frequency 
shifts on the modulation due to motion of the remote 
transmitter or transponder and is usually quite small 
in comparison to the modulation frequency, so that no 
serious filtering problem arises. 

6. Delay associated with the reflection is assumed to 
be constant and independent of time. This means that 
the remote transmitter or transponder is not moving. If 
the distant transmitter were moving, it would be neces- 
sary to consider the doppler effect and consequently the 
frequency difference between the doppler-shifted direct 
and reflected signals. In the case of reflection from the 
ground plane, the paths of the direct and reflected sig- 
nals are nearly parallel when the distance from the 
transmitter to the receiving antenna is large compared 
to distance from ground plane to antenna. Both direct 
and reflected signals are then shifted in frequency almost 
equally so the frequency difference is quite small. 

In the computation of signal distortion, it is entirely 
feasible to include transmitter motion; however, as dis- 
cussed later, for the highest velocities encountered, the 
contributions to phase error by doppler effects are in- 
deed quite small and, therefore, ignored. 

The signal at the receiving antenna is the result of 
direct and reflected signals. Let the field strength of the 
direct signal at the receiving antenna have an amplitude 
Eo, and the field strength of the reflected signal at the 
receiving antenna an amplitude £,. The amplitudes Eo 
and £, contain the effects of the antenna radiation pat- 
tern and the separate space attenuations of direct and 
reflected signals. The amplitude £, contains the reflec- 
tion coefficient of the ground plane as a factor. So far as 
this discussion is concerned, these effects can be ex- 
pressed by writing the ratio of reflected-to-direct signals 
as E/E) =R exp (jg), where Rand ¢ may be thought of 
respectively as the amplitude and phase shift of an effec- 
tive reflection coefficient. As stated earlier, both Rand ¢@ 
are assumed independent of frequency. It is also as- 
sumed that R and ¢ are independent of time. Travel 
time of the direct signal will be designated by #4 and 
time of the reflected signal by ¢.. The multipath deiay is 
then, fg -—f=to. The frequency-modulated signal for the 
direct path may be expressed as 


—t 
éo = Ey exp lif ads] = Eo exp (jb1), (1) 
0 


where the instantaneous frequency, Q, is defined by 


Q = w, + Aw cos wil, 


with Aw = 27x frequency deviation and w; = modulation- 
angular frequency. For reflected signal, field strength is 


tte 
é, = EyR exp] f Qdt + 6| = E,)Rexp (jbe). (2) 
0 


It is assumed that eo and e, may be combined at the 
antenna by phasor addition to produce a resultant sig- 
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nal at the receiver input. Fig. 3 shows the phasor addi- an antenna pattern null. Unfortunately, this appears to 


tion to the direct and reflected signals. 
The instantaneous phase b of the resultant is 


R sin (bs 1 bi) 
1+ Rcos (b — 0;) , 


b = 6; + arctan 


(3) 


Fig. 3—Phasor addition of direct and reflected signals. 


Assuming the use of a perfect limiter, the output of the 
discriminator is the instantaneous frequency given by 


Resin (65 — 6 
— = — + a arctan ee 2 (4) 


di dt d 1+ Ros (b. — 6,) 
or 

db 

Fas + Aw cos wi(t — 41) 


R sin (by — by) 
1+ Rceos (b2 — b;) 


(S) 


+ — arctan 
dt 


The dc term w, need not be considered. The desired out- 
put is the second term, Aw cos w;(t—f), of which the 
phase lag wit; is to be measured by comparison with a 
reference signal. The last term represents nonlinear dis- 
tortion. We wish, then, to find the magnitude of the 
component of distortion at frequency a. All other com- 
ponents are assumed to be removed by a perfect filter. 
It is now convenient to make the following substitu- 
tions: 


69 = weto + d = phase delay of the carrier, (6) 
B = Aw/w, = frequency deviation ratio, (7) 
B = w(t — th — t/2), (8) 
z = 26 sin wylo/2. (9) 


With these substitutions, it can be shown that 


bo = by = A —"Z% COS 164, (10) 

Corrington* has developed the expansion of the de- 
rivative of the arctangent function of (5) which con- 
verges for R<1. It may seem unrealistic to restrict R 
to values less than unity since it is quite possible for R 
to be greater than unity as, for example, in the situation 
in which the direction of the transmitter coincides with 


3M. S. Corrington, “Frequency modulation distortion caused by 
multipath transmission,” Proc. IRE, vol. 33, pp. 878-891; December, 
1945. 


be a particularly bad situation for which there is no 
mitigation other than to avoid it altogether. Because the 
FM receiver tends to suppress the weaker of two signals, 
the measured phase closely follows the phase of the re- 
flected signal when R>1, obviously an undesirable cir- 
cumstance. Therefore, in order to prevent the measured 
phase from oscillating grossly between the phase of the 
direct signal and the phase of the reflected signal, care- 
ful design of the antenna and siting must be used to pre- 
vent & from becoming greater than unity. Starting with 
(3) and using the identity tan y/2=sin y/(1+cos ¥), it 
is not difficult to show that when R=1, the phase error 
is Wito/2, or one-half the phase delay on the modulation 
frequency. 

Following Corrington, the expansion of the arc- 
tangent in (5) with substitutions of (6), (7), (8), and 
(9) is 


R sin (@ — Z cos B) = 
1+ Rcos (9 —Zcos B) — a n 


arctan sin 769 


| Jains) +2 ss (—1)"Jom(mz) Cos 2mB | 
io] = IRE 
[pee 


n=0 112 


COS 20 


[2 SS (-1) J omsi(nz) cos (2m + a], (11) 
m=0 
where Jon(nz) and Jom41(z) are Bessel coefficients. 

In the process of differentiating (11) to find the signal 
distortion components, the quantities 60, z, and B may 
be considered functions of time by virtue of motion of 
the transmitter which causes the quantity ft) to vary 
with time. The phase B also contains time explicitly. 
Considering f) to be a function of time, the differentia- 
tion yields terms proportional to w,(dto/dt), Bun(dto/dt), 
w;(dto/dt), and w,. In the case of ground reflection on the 
assumption of parallel rays, the time delay fo is nearly 
2h/c sin 6 where h is the antenna height, c the velocity of 
propagation, and @ the elevation angle. If the elevation 
angle is changing at the excessive rate of one radian per 
second, dfy/dt is less than 2h/c. Suppose, for a possible 
tracking system, that 


h = 1 meter 
Ou 
@ y= 108 
6 = 10. 


Then, 
wedto/dt = (2/3)(10) 
Boordto/dt = (2/3)(10)(10°/108) = (2/3)(10~) 
wdly/dt = (2/3)(108/108) = (2/3)(10-), 


so that all terms proportional to dto/dt are negligible in 
comparison to the term proportional to w;(=10°). It 
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should be noted also that from (4) the dc term , would 
properly appear as w,(1 —dt,/dt) asa result of the motion 
of the transponder. This term, having only very low- 
frequency components corresponding to the motion of 
the transponder, is filtered from the output. The result, 
after retaining only the term proportional to 1, is 


d R sin (@) — 2 sin B) | 
= t 
dt E F 1 -+ Rcos (6 — 2 sin B) 
= 01 >, Smsin 2mB + 01 >) Cm sin (2m + 1)B, (12) 
m=1 m=1 
where 
Sin(Oo, R, 2) = (—1)™4m 
co — Ix n 
>> oy sin (160)J om(1z) (13) 
n=1 n 
Cm(90, k, Z) ae (= 1) ™2(2m = 1) 


Cake 


n 


cos (100)J om-+1) (nz). (14) 


foo) 
n=1 


The term in the expansion at frequency a is given by 
m=O, so that the output of the discriminator and filter 
is 
output = Aw cos wi(¢ — t1) + wiCo sin B 

= Aw|cos w(t — #:) + Co/B sin ai(t — ty — to/2) | 


= AwA cos [w(t — ty) — al, (15) 


where 
Co/B cos wito/2 


fot paral aie (16) 
1 — Co/B sin wyto/2 


a = arctan 


is the desired expression for the phase error due to multi- 
path interference in this FM system. 
The coefficient Cy is given by 


n=1 nN 


Oy x% 
“Ji E sin (= =), 
oy 2 


where a new parameter has been introduced as x =o, 
which is the difference in phase of the carrier for the two 
paths. 


(17) 


INTERPRETATION 


The following observations can be made about C). 
First, Co is periodic at both the carrier and modulation 
periods. Second, Cy has the form of a rapid oscillation 
modulated by a slowly varying envelope represented 
by the Bessel function (the ratio w/w, is normally quite 
small in a tracking system). Third, Cy is bounded inde- 
pendently of 8, so that the ratio Cy/8 becomes smaller as 
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»O, Dds) 


6B becomes larger. This follows from the fact that series 
for Cy is absolutely convergent, and its sum less than 


(oe) R* 
a>) 
n=l 1 


which is convergent for R less than unity and inde- 
pendent of 6. 

Fourth, when the ratio of modulation frequency to 
carrier frequency is small, and the argument of the 
Bessel function may therefore be considered as slowly 
varying and, for practical purposes, constant, the func- 
tion Cy has its maximum and minimum values at the 
points where (x—@) is zero or a multiple of 7. This can 
be seen by differentiating Co with respect to x, holding g 
constant in the Bessel function. It may thus be inferred 
that the envelope functions of Co are given very nearly 


by 


Co = — 23 (= 1)*(Re/n) Tile) (18) 
and 
Co=—2 S (R"/n)J (nz). (19) 


n=1 


Graphs of these functions are shown in Fig. 4 for several 
values of R. Using these envelope functions for Co, the 
envelope curves of phase error as a function of tp can be 
computed and are shown in Fig. 5 (following page) for 
R=0.9 and several values of B. 

A qualitative picture of the reduction of phase error 
as the deviation ratio is increased can be obtained by 
plotting the maximum values of the curves in Fig. 5 
versus 8. This plot is shown in Fig. 6, next page. The 
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Fig. 5S—Envelope of modulation-frequency phase error for FM as a 
function of phase delay, ue =mto of indirect signal modulation. 
Tae aay B=5, oa Sane aia B=1. 
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Fig. 6—Phase error at the first maximum of the phase error envelope 
as.a function of 8. — — — — lower envelope, upper en- 
velope. 


curve is roughly of the form 1/8, and that about 90 per 
cent of the error reduction takes place for f less than 10. 


MULTIPLE MODULATION 


When more than one modulation frequency is present, 
the phase-error formula on each frequency still has the 
same form as (16), but the parameter corresponding to 
Co is then defined by a different series which contains 
“cross-talk” terms between the modulation frequencies. 
The distortion term for the case of two modulation fre- 
quencies having the relationship w,= 10a is derived in 
the Appendix. The result after the harmonics have been 
filtered out is 


distortion = w3C1 sin By, + woe sin Bo, (20) 


where 


By, — w(t = ty = to/2) (21) 
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* EXPERIMENTAL 
o COMPUTED 
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Fig. 7—Phase error vs multipath delay time 8=10, R=0.5. 
Delay in degrees of modulation frequency. 


Bo == wo(t a ty as to/2), (22) 


and C, and C; are defined by the two series (40) and (41) 
of the Appendix. 
The receiver output is proportional to 


Aa, cos w(t ty) be w1Cy sin By 


+ Aw» cos welt ree ty) + w2C 2 sin By. (23) 


Following the procedure used in the previous section for 
a single modulation frequency, the phase error for each 
modulation frequency can be written immediately. For 
the w; component, the phase error is 


C/B1 cos (wto/2) (24) 


a, = arctan = ’ 
1 — Ci/Bx sin (w1 to/2) 
and for the wz component, the phase error is 
Cc COS (Wolo /2 
a, = arctan 2/8 C8 (25) 


{Cyr 


It is shown in the Appendix that both C, and C, are 
bounded independently of 8: and 6.. That the coeffi- 
cients C, and C, are also periodic in the carrier and 
modulation periods can be seen readily by inserting (28) 
in (40) and (41). Thus, it is observed from the form of 
C, and C, that the behavior of the phase errors is similar 
to the single modulation-frequency case. Since envelope 
curves for C; and C; are not easily derived, computation 
of phase error envelopes was not undertaken. It is 
sufficient, however, to note that, since C,; and C, are 
bounded, the phase errors a; and a, can be made as 
small as desired by sufficiently increasing the deviation 
ratios. 


EXPERIMENTAL DATA 


Fig. 7 shows a plot of experimental data for a single 
modulation frequency obtained with a mercury delay- 
line multipath simulator. The experiment was per- 
formed by Mr. Leonard Hawn of The Ralph M. Parsons 
Company Laboratory. Curves in Fig. 7 are the phase 
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error envelopes; therefore, periodicity at the carrier fre- 
quency is not shown. Computed curves of the envelopes 
obtained from (18) and (19) are also represented in the 
figure by the large circles. Although the agreement be- 
tween the experimental data and the computed data 
points is not precise, it is quite satisfying. 

The multipath simulator consisted of a transmitter, 
receiver, phase meter modulation source, and two paths 
between the transmitter and receiver. A calibrated 
mercury delay line provided a means of delaying the 
signal in one path by a known amount. An attenuator 
was inserted in the other path to enable the ratio of sig- 
nal strength of the two paths to be adjusted to a desired 
value. 

The transmitter consisted of a Colpitts oscillator 
operating at a center frequency of 1.0 mc, a reactance 
tube, a limiter, and impedance matching amplifiers to 
the two paths. Double-shielding and low-pass filters in 
the power-supply leads and modulation-frequency con- 
nections were used to prevent leakage to the receiver. 

The receiver consisted of a preamplifier, adder, S- 
meter, tandem limiter, wide-band discriminator, rf filter, 
audio amplifier, cathode follower, and output low-pass 
filter. The preamplifier compensated for insertion loss of 
the delay line and partially corrected for frequency and 
phase response of the delay line. The discriminator de- 
sign was based on relations given by Granlund,* who has 
demonstrated that, in order for the FM receiver to be 
effective in suppressing the weaker of two signals, the 
discriminator bandwidth must exceed the deviation by 
the factor, (1+ R)/(1—R). 

The modulation frequency used was 1000 cps. Effec- 
tive filtering of distortion components was obtained 
with an M-derived 1-kc filter having rejection fre- 
quencies of 2.0 and 3.0 ke. 


CONCLUSIONS 


Expressions for the phase error introduced by a single 
multipath reflection have been derived for the cases of 
one and two modulation frequencies. Experimental data 
obtained with a multipath simulator show good agree- 
ment with computed values of phase errors. Study of the 
behavior of the phase errors for these rather idealized 
cases substantiates the prediction that the use of wide- 
band frequency modulation can result in a significant 
reduction of phase error due to multipath signals. The 
improvement over the use of amplitude modulation is 
illustrated by comparison of Figs. 2 and 5. For 8=1, the 
phase error in the FM case has about the same order of 
magnitude as in the AM case. As @ is increased to a 
value of 10, the phase error is decreased about one-fifth. 
For a cw tracking system having a highest modulation 
frequency of 1 megacycle, the frequency deviation cor- 
responding to B=10 is 10 megacycles and requires a 


‘J. Granlund, “Interference in Frequency Modulation Recep- 
tion,” MIT Res. Lab. of Electronics, Tech. Rep. No. 42, p. 8; Janu- 
ary, 1949, 
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bandwidth at least this amount. Although a further in- 
crease in frequency-deviation ratio can result in further 
reduction of phase error, the necessary bandwidth may 
become excessive when considered from the point of 
view of broadband design and spectrum availability. 


APPENDIX 
MUvuLTIPLE MODULATION DERIVATIONS 


If there are N modulation components, the instan- 
taneous angular frequency of the frequency-modulation 
signal is 

Q =o. + ys Ayw COS wyt, (26) 
N 
where Ayw is the angular deviation associated with the 
Nth-frequency component. Following the same pro- 
cedure used in developing the resultant phase for one 
modulation frequency, 


bo — bi = wey + & — >, Zw COS By, (27) 
N 
where 
t 
7 ey yee ae (28) 
and 
to 

By = wn i) (29) 


Using the series expansion of the arctangent® and in- 
serting (27) and (28), 


R sin (be —< b1) 
1+ Rceos (b2 — 63) 


cc) R* 
= — >) (=1)* sin AG — }) zy cos Br). (30) 


n N 

The right side of (30) can be expanded further, using 
the appropriate trigonometric identities until the indi- 
vidual harmonics of By appear; however, the result is 
very complicated, and the task of sorting out the terms 
of a single frequency is extremely tedious if the number 
of modulation frequencies exceeds two. To illustrate the 
result, only two modulation frequencies are considered. 
Eq. (30) is first differentiated to find the instantaneous 
frequency. The result of the differentiation is 


arctan 


d | R sin (b2 — 8,) 
—| arctan ————_—__ 
dt 1+ Rcos (be as ra 


= ))(—1)"R" cos n(4) — 21 cos B, — Z2 CoS By) 


: Eon sin By + ZoWs sin By). 


(31) 


Using the identity for the cosine of the difference of two 
angles, (31) becomes 
SE. T. Whittacker and G. N. Watson, “A Course of Modern 


Analysis,” 4th ed., Cambridge University P : 
land, p. 190; 1952. ge University Press, Cambridge, Eng- 
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“A 
i) aretan 
dt 


R sin (bo a by) 
1+ Ros (by — =| 


a Dy fet )eRy [cos n(@) — 21 cos By) cos (nz2 cos By) 


n=1 


+ sin 7(8) — 21 cos By) sin (mz cos Be) | 
- [2101 sin By + zgw2 sin Bo] = — D> (—1)"R» 
n=1 


- [211 sin By cos n(8) — 21 cos By) cos (nzy cos By 
++ 1, sin By sin 2(4) — 2; cos Bj) sin (nze cos Be) 
+ Z2W2 Sin Be cos (8 — Z2 cos By) cos (nz, cos Bi) 
++ Zw2 sin By sin n(8) — 22 cos By) sin (nz; cos By) J. (32) 


Four identities are useful at this point: 


cos (nz cos B) = Jo(nz) + 2 D5 (— 1)*J2x(nz) cos 20B (33) 


v=1 


—2 ss (—1)*Joo41(nz) cos (20 + 1)B (34) 


v=0 


Z COs os 


sin (zz cos B) = 


sin B cos n(@ — 


=2(-p-= 


Jom4i(nz) sin [(2m + 1)B] cos (260) 


apy re 1)" — 


m=1 


— Fam(12) sin 2mB sin (6p) (35) 


sin B sin n(@) — z cos B) 


=2>(-1)"— 


m=0 


Samal sin [(2m + 1)B] sin (16) 


2m 
—2 Ss (- Lite — 5 J omknz) sin 2mB cos (4). 


m=1 


(36) 


The third identity, (35) is proved by Corrington,* and 
proof of (36) can be established by following his method 
as presented. Introduction of these four identities in (32) 
leads directly to 


d | R sin (be — bi) | 
2) GE RAR, Se ee ee 
dt 1 + Rcos (be = bi) 


= & 2 il 
7m > (—R)"Z101 {2 SD C1) sida J (Z2) J 2m41(n21) 
x a [(2m + 1)B,] cos ye 
1 
4 3 3 — 1) (mete) J 2u(22) J 2m41(221) 


m=0 v=1 


X cos 2vB, sin [(2m + Pee cos (169) 


= 2 . . 
+2 se (—1)” cst To(2)Jom(nz1) sin 2mB, sin (n6o) 


m=1 N21 
+4000 (- spite = = Fae(ts)Tan(ne) 
m=0 v=1 


X cos 20Be sin 2mB, sin (740) 
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m=0 v=0 NZ, 


X cos [(20 + 1) Be] sin be + 1) By] sin (16) 


440 (tyne TE 


m=1 v=0 BY 


J o041(032) I am41(n21) 


J ov41(122) IJ om(21) 


X cos [(20 + 1)B2] sin 2mB, cos my 


Same as above except with 2; 


— ee R)"Zewe}and z. interchanged and B, 37) 
= and Bz interchanged 


Inspection of (37) reveals that there is an infinite num- 
ber of terms at the fundamental frequencies, w; and w.. 
These terms are produced when the beat frequencies, 
which may be designated by kwi+qwe, equal either a; 
or we. The numbers k and q stand for the integers 2v+1, 
2v, and 2m, whenever applicable. No attempt will be 
made here to follow through the long process of sorting 
out the fundamental components, and only the first few 
terms for the condition that w.=10e,;, will be written 
down to illustrate the result. Thus, the following is ob- 
tained from the distortion term, (37). 

Fundamental Component of Distortion = 


= 204 si St 


n=1 


J o(n2)J1(nz1) cos (Ao) 


— 2J,(z2)J10' (21) sin (760) 

+ 2J2(nz2)J20'(nz1) cos (no) 

a 2I (12) J 30’ (2:1) sin (16) 

— 2F4(nz2)J 40'(nz1) cos (74) 

— 2F50(nz1)J5/(nz2) sin (26) 

+ an infinite series of higher order terms 
— 2wW» 5: Sak 


n=1 


Jo(nz1)J1(nz2) cos (280) 


— 2F10(mz1)J1'(nz2) sin (780) 

+ 2Fo0(nz1)J 2! (nz2) cos (760) 

+ 2 30(nz1)J3/(nze) sin (16) 

— 2F 49(nz1)J 4/(nz2) cos (nO) 

— 2F 50(nz1)J50'(nze) sin (76) 

+ an infinite series of higher order terms, (38) 
where J’ denotes the first derivative of the Bessel func- 
tion with respect to mz, or mzZ2 as required. 


The fundamental component of distortion, therefore, 
has the form 


w1Cy sin By, + woC 2 sin By. (39) 


With the introduction of Jo’(z) = — Ji(z), the coefficients . 
C, and C, can be written compactly as double summa- 
tions thus: 
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Cis =e Ds 
n=1 


Ear 


nN 


[-> (—1) *exJox(nz2)J’20%(021) COS (Ao) 


k=0 


—2 sy (- 1) *Fons1(122) I’ 10 (2K41) (M21) sin (nt) | (40) 


k=0 
and 
oe) — jf n 
eS: (>) 
n=1 n 


[- »s (—1) ®exJ 20%(Z1) J 2%(mZ2) COS (Oo) 


k=0 
— 2 5 (-1) "Siro eeg1y (021) 2n41(mz2) sin (nt), (41) 
k=0 


where e; is equal to two when & is not zero and e; is equal 
to one when & is zero. It is desirable to show that both 
C, and C, are bounded independently of 2: and zg,. Con- 
sider first the series for C;. Each term of the absolute 
series for C; is less than the corresponding term of 


) JRE ie) 
Ss —| >> ex| J ox(nZ2) | | J" o0x(n21) | 
n 


k=0 
+2 MS J ox41(n22) | | J" 10 2b41) (21) | ; (42) 
k=0 


Since J_n=(—1)"Jn, this series can be written con- 
veniently as 


n=1 
Using Cauchy’s inequality® 


Peer) 


k=—0 


R” e.) 
ye, | Ji(nz2) | | J" s0x(n21) | : 


N k=—co 


(43) 


6G. H. Hardy, “A Course of Pure Mathematics,” 9th Ed., The 
Macmillan Co., New York, N. Y., p. 490; 1946. 
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y | Je(n0z2) |?- Ss | (J10%(m21) |? 


Z 
= Y Vilnes) Eo Vaote (44) 
Now? - fe 
Ss Gian (45) 
and ee 
x [7sou(s) 2 < = [W@P=1. (46) 


This last inequality is obvious since the sum of the series 
formed by using only every 10th term is less than the 
sum of the series over all values of k. It is not difficult to 
show that 


dy Ive)? =1 
k=—2o 
by use of the recurrence formula 
2F«/(z) = Ju-1(2) — Juts(2), 
and the relations 


ss Ji(z)Ju42(z) = 0 (Ref. 7), and S3 J2@) Sai 


k=—o k=—oo 


Thus, it is clear that 


foo) 


DA, Ti(nz)s | | J" ox(mzi1) | <1 (47) 
k=—w20 
and, therefore, 
co Wee 
| C,| <2 S. (48) 
n=1 nN 


which is convergent for R<1. 
An identical result applies also to Cy. 


7G, N. Watson, “Theory of Bessel Functions,” The University 
Press, Cambridge, Mass., p. 31; 1952. 
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Application of the Reaction Concept to Scattering Problems* 
M. H. COHEN+ 


Summary—Rumsey’s “reaction concept” is applied to scattering 
problems by considering volume polarization currents which generate 
the scattered field. The scattering cross section is discussed. Ap- 
proximation formulas are developed and illustrated by calculations 
for the echo width of a circular cylinder as a function of dielectric 
constant. Calculations are made for four plane-wave type trial fields 
in the cylinder; these are compared with the exact solution and with 
integral equation solutions through the second order. In this case the 
reaction method gives better results than the integral equation 
method, for the same amount of labor. 


INTRODUCTION 


N A RECENT paper Rumsey! defined an electro- 
| magnetic reaction and showed why it is to be re- 

garded as a basic physical observable. His argument, 
briefly, is as follows. Any microwave measurement con- 
sists of measuring, not the electric field (force on a unit 
charge) at a point; but rather, the signal at the termi- 
nals of an antenna. Thus two antennas, for transmission 
and reception, are inherent in a measurement. (In the 
radar case these may degenerate into the same antenna.) 
The antennas may be considered as sources of two dis- 
tinct fields, and the reaction between these sources is 
the observable quantity. 

Let a and b represent the two sources,? consisting of 
distributions of electric currents J(a) and J(b), and 
magnetic currents K(a) and K(d). The fields radiated by 
these sources are E(a), H(a), and E(b), H(b). Then the 
reaction between the sources is denoted by the symbol 
(a, b) and it is defined by 


(a,b) = f f if [5@ JO) — H@)-K@Ieo. (1) 


The integration volume contains J(b) and K(é). Under 
suitable conditions the reciprocity theorem holds: 


(a, 6) = (6, a). (2) 


It can be shown that, as a phasor, the reaction (a, b) 
represents half the time-dependent portion of the in- 
stantaneous power expended by source a on source b. If 
source } is a unit current source, then (a, 0) represents 
the open circuit voltage received at b when a alone is 


* Original manuscript received, October 12, 1954. Portions of this 
paper were presented at the Washington Meeting of the American 
Physical Society, April 29, 30, and May 1, 1954; and at the joint 
URSI-IRE-AGU Meeting in Washington on May 3-6, 1954. It is 
partially contained in Report 475-14, Antenna Laboratory, Ohio 
State University, June 30, 1954. ; j 

+ School of Electrical Engineering, Cornell University, Ithaca, 
N. Y.; formerly Antenna Laboratory, Ohio State University, Co- 
lumbus, Ohio. : j ann b 

1V. H. Rumsey, “Reaction concept in electromagnetic theory, 
Phys. Rev., vol. 94, pp. 1483-1491; June, 1954. The notation used in 
this paper follows that used by Rumsey. 

2 It is implied that both sources oscillate at the same frequency. 
The time convention exp (—iwt) is used. 


radiating. Further significance of the reaction is dis- 
cussed by Rumsey.! 

An important application of the reaction concept con- 
sists in finding suitable approximation formulas in prob- 
lems where the exact solution is not known (or is useless 
because of computational difficulties). If one has an 
assumed current which generates an approximation to 
the desired field then parameters in the assumption can 
be “fixed” by requiring, as far as possible, that it have 
the same reactions with other sources that the correct 
current does. This generally gives the customary sta- 
tionary, or variational, formulas. 

In this paper the application of the reaction concept 
to scattering problems will be discussed by using volume 
currents which generate the scattered field. This follows 
directly from Rumsey’s work,! although he confined 
himself to surface currents. Approximation formulas 
will be developed and illustrated by calculations for an 
infinite circular cylinder. The value in these calculations 
is that the exact solution for this problem is known, and 
thus some measure of the worth of the different approxi- 
mations and current assumptions is obtained. It is 
hoped that this will furnish insight into suitable ap- 
proximation procedures to be used for other problems. 


Fig. 1—Scattering problem. 


SCATTERING BY A FINITE IsoTROPIC BopDy 


Let a source of monochromatic electromagnetic 
energy at g (Fig. 1) illuminate a finite isotropic scatter- 
ing body with the bulk constants yp, e, and a. Let E, E’, 
and E* represent the total, incident, and scattered elec- 
tric fields, respectively. The incident field is that ob- 
tained in the absence of the scatterer, and the scattered 
field is defined by 


Et =E-— E‘ (3) 
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It can be shown by a direct manipulation of Maxwell’s 
equations that the source of the scattered field may be 
taken to be the volume polarization and conduction cur- 
rents radiating in free space. In general there are both 
electric and magnetic currents, given by 


J = — tole — @)E + cE, (4) 
and 
K = — iw(u — bo). (5) 
For convenience, let 
K == — tle — €9) 7:0 (6) 
and 
y = — tw(u — bo), (7) 
so that 
J = KE (8) 
and 
K = oH. (9) 


This representation for the source of the scattered field 
is different from the surface currents used by Rumsey.! 
Some consequences of the difference are discussed. 

We wish to investigate the field at # (Fig. 1) and so 
we put an antenna there, according to the above dis- 
cussion. Let J(h), K(h) represent a distribution of 
sources which generates the field obtained when the an- 
tenna at h is radiating. Also, let c and d represent the 
sources of the scattered fields obtained when the sources 
g and h, respectively, are radiating. That is, 


J() = xE(g) J(d) = xE(h) 
K(c) = yH(g) K(d) = yH(h). 


Then the measure of the total field at 4, due to the 
source g, is the reaction (g, 1, h), where the internal 
numeral 1 means that the sources radiate in the presence 
of the scatterer. Similarly, the measure of the scattered 
field at hf is the reaction (c, 2, h), where the internal 
numeral 2 means that the sources radiate in the ab- 
sence of the scatterer. The formula for (c, 2, h) is 


(10) 


(6, 2,h) = f f f [B(c)-J(h) — H(c)-K(h) do, (11) 


where the integration volume contains J(h) and K(h). 
By the reciprocity theorem, 


(c, 2, hy = {h; 2, ¢) 


= [J [ 2o-s20 - wuy-Kole 
= [ff ee c-2@ - a U)-H@) I, (02) 
where the integration is over the volume of the scatterer. 


The reaction (c, 2, #) is proportional to the scattered 
signal received at h. When E(g) and H(g) are known (or 
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can be approximated), the reaction can be computed by 
(12), since E*(h) and H‘(h) are characteristic of the re- 
ceiving antenna. 

Consider (h, 2, c). It may be regarded as the signal 
induced at the “terminals” of the source ¢, by the inci- 
dent field due to 4. From (3) then, we have 


CAO RF Me PP) ES KY 10 WOES (13) 
and similarly, 

(g, 2, d) + (c, 2, d) = ¢g, 1, d) (14) 

(h, 2,.g) + (d, 2, 8) = (h, 1, g) (1S) 

(g, 2, h) + (c, 2, hk) = (g, 1, h). (16) 


From elementary combinations of these formulas (and 
the reciprocity theorem), the following reciprocal rela- 
tions are derived: 


(17) 
(18) 


(Gu2, #) = (a, 2, 2) 
{c, 1, h) = (d, 1, g). 


SCATTERING Cross SECTION 


The scattering cross section is ordinarily defined as 
the cross section of an isotropic scatterer which gives 
the same scattered intensity as the target does, for an 
incident plane wave. For a complete description, it is 
considered to be a two-by-two matrix with components 
A,x, where 


A,x = Lim 4nR,2 


Rr 


(19) 


| Ex! 


The subscript k refers to orthogonal linear components 
in the plane normal to the incident propagation vector; 
the subscript p refers to orthogonal linear components 
in the plane normal to the scattered propagation vector 
(along the line from the target to point h, Fig. 1). R, is 
the distance from the target to the point A. 

Let the source of the incident wave be the source g, a 
distant dipole with components J;. The polarization 
currents c; are those induced by E,. Let the source h bea 
distant dipole with components J,. Then the reaction 
(¢, 2, hk), proportional to the p-component of scattered 
field for a k-polarized incident wave, is given by 


Sermak wie f f ff E*(I;)-J(h,)dv = E,*(In)I». (20) 


As R,-> ©, | E,'| >wpol,/(42R;), at the target. Hence 


= (ono)? | (chy 2, hp) |? ; 


Aye = 
* 4e| Ex P?] E,' see 


z 


By (17), the scattering cross section satisfies the recipro- 
cal relation: 


A pk i. Aires (22) 


APPROXIMATIONS 


The “reaction method” of approximation consists of 
assuming functions which approximate the polarization 
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currents ¢ and d and forcing them, as far as possible, to 
have the same reaction relations that the correct sources 
do. The adjustable parameters contained in the assump- 
tion are thereupon fixed. 

The reaction relations satisfied by c and d are con- 
‘tained in (13) through (16). They may be generalized to 


ney oy Reed, y(t be) (23) 
and 


(h, 2, y) + (d, 2, y) = (h, 1, y), (24) 


where x and y are any sources. The approximations to 
the sources ¢ and d, denoted by sources a and b, are to 
be substituted into these equations. Adjustable con- 
stants in the approximations are fixed when sources are 
substituted for x and y. 

Now consider what the possible substitutions for x 
and y are. There are six sources in the problem: a, ), c, d, 
g, and h. Two of these, c and d, are unknown; and g and 
h cannot be used because (g, 1, g), (g, 1, 2), and (h, 1, 2) 
cannot be calculated. Therefore we must use a or b. But 
they can be subdivided into component electric and 
magnetic currents; reactions can be computed sepa- 
rately for them, so there are four possibilities for x and 
for y. We now fix the choices for x and y by noting that 
there are two approximations to the reaction (c, 2, h), 
obtained by substituting a and 6 for c and d in the two 
sides of (17). These are the same when we substitute 
J(o) and K(6) for x in (23), and J(a) and K(a) for y in 
(24), so this is the choice. 

When these substitutions are made into (23) and (24), 
we have 


{g, 2, J(b)) + (a, 2, J(b)) = (g, 1, J(4)) 
(g, 2, K(b)) + (a, 2, K(®)) = (g, 1, K(6)) 
(h, 2, J(a)) + (6, 2, J(a)) = (h, 1, J(a)) 
(h, 2, K(a)) + (6, 2, K(a)) = (h, 1, K(a)) 
The total fields of sources g and h are needed in the 
right-hand sides of these equations. They are unknown, 


but by the requirements implicit in (8) and (9), they 
must be proportional to: the assumed currents. That is, 


E(g) = «*J(Q) ~ « J), (26) 


and so on. In integral form, then, (25) become 


(g, 2, J(6)) + (a, 2, J(B)) = { { He y(a)-J(b)do (27) 


(25) 


(g, 2, K(b)) + (a, 2, K(b)) = { { if »K(a)-K(8)do (28) 
(h,2, Jla)) + (6, 2, J(a)) = i f [so Toe (29) 


-(h,2, Kl) + 2, K@) = [ ff 7°KO-Ke@e. (G0) 


Eqs. (27) through (30) are four simultaneous equa- 
tions, from which four parameters, one each in J(a), 
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K(a), J(0), and K(6), can be determined. By adding (27) 
to (28), and (29) to (30), it is shown that 


(a, 2, h) = (b, 2, g). 


Thus the two approximations to (c, 2, 4) are equal. 

In many cases the distributions of currents may be 
reasonably guessed at, whereas their Jevels (phases and 
amplitudes) with respect to the incident fields cannot. 
To consider this case,? let 


(31) 


J(a) = ax® 
K(a) = 6cO 


J(b) = 6c" 


K(b) = nV, (32) 


where ®, ©, W, and Y¥ are functional approximations to 
the total fields, and a, 8, 6, and 7 are adjustable con- 
stants. Then, from (27) and (28), 


[fffze-ore + ff [eons 

-{ jprermn 
[ff @-me+ ff [ one 

= [J [sovee. (9 


Here E‘(a) and H*(a) must be computed from the as- 
sumed currents J(a@) and K(a) which radiate in free space. 

Eqs. (33) and (34) are a simultaneous pair for the de- 
termination of a and @. A similar pair for the determina- 
tion of 6 and 7 are obtained from (29) and (30). The de- 
sired reaction is now given by 


(c, 2, h) = (a, 2, h) = (h, 2, a), 
= aEi(h)- — BH*(h)- dv. 35 
(h, 2, 0) fff (i) — BH'(h)-y@]dv, (35) 


The above approximation procedure differs from 
Rumsey’s “testing” technique! by the necessity for 
making the integral approximations to the right-hand 
sides of (25). This comes about because we have used 
polarization currents rather than the equivalent surface 
currents used. by Rumsey. In general there is no clear 
choice for the more useful representation, in terms of 
surface or volume currents. However, in the important 
special case when p= so the magnetic polarization cur- 
rents are zero, and we would have to make only one 
guess, for the distribution of electric field, to use the 
volume currents. On the other hand, to use free-space 
calculations, both electric and magnetic surface currents 
are needed (except for perfect conductors), and to use 
them we would have to guess at the distributions of 
tangential electric and magnetic fields. 


Higher Order Approximations 


Higher order approximations may be made by assum- 
ing that the currents a and b are combinations of ele- 


3 This is the only case considered by Rumsey. 
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mentary distributions. The phase and amplitude of each 
elementary source are to be found from (25)imlet 


J(a) = YS ane, J(b) = 2, oP 
Koy = Sie hO pe KB) =a vermense) 
j=1 T= 


Then, from (25), 
(g, 2, Ws) + (a, 2, Wx) = (2,1, We), R= 1,--- 0, (37) 


and 


dg DEN) ta, 27 Vials LAV) Pali sl aia ate: (38) 
Eqs. (37) and (38) are a set of 2” equations for the 2 
unknown coefficients a; and 8;, and the coefficients 6, 


and 7; are similarly found. 


Stationary Properties 


It can be shown that the reaction when computed by 
(35) is stationary with respect to first-order variations 
of the total fields about their true values. Eqs. (37) and 
(38) result when the fields are expanded as in (36) and 
the Ritz procedure is used to find the coefficients. 

Since the work of Levine and Schwinger,’ the varia- 
tional method has been widely used in problems of 
scattering by conductors,** and in a few problems con- 
cerning dielectrics.!7 In standard derivation, motivation 
of the manipulation of the equations to arrive at sta- 
tionary form is obscure, whereas the reaction concept 
leads one naturally and pictorially to it. It enables one 
to proceed when there is a variety of stationary forms.} 
It also is more general, since discussion is not limited to 
far fields; characteristics of transmitting and receiving 
antennas are explicitly introduced. 


APPLICATION TO A CIRCULAR CYLINDER 


As an example of the use of the above theory, we shall 
compute the back-scattering cross section of a lossless 
circular cylinder. The cylinder has a diameter equal to 
one free-space wavelength (a/Ay =4); it has unit relative 
permeability and its relative dielectric constant is e,. Its 
axis is the z-axis, as in Fig. 2. A plane wave, e*, is inci- 
dent from the left and has its electric vector parallel to 
the z-axis. By symmetry the total field also has only a 
z-component. There is therefore only one component in 


4H. Levine and J. Schwinger, “On the theory of diffraction by an 
aperture in an infinite plane screen,” Phys. Rev., vol. 74, pp. 958— 
974; October, 1948; vol. 75, pp. 1423-1432, May, 1949, 

5C. H. Papas, “Diffraction by a Cylindrical Obstacle,” Jour. 
Appl. Phys., vol. 21, pp. 318-325; April, 1950. 

§R. G. Kouyoumjian, “The Calculation of the Echo Areas of 
Perfectly Conducting Objects by the Variational Method,” Ohio 
ae University Res. Found. Antenna Lab. Rep. 444-13; November, 

953. 

7J.R. Lurye, “Electromagnetic Scattering Matrices of Stratified 
Anisotropic Media,” Rep. EM-11, N. Y. Univ. Math. Res. Group; 
May, 1951, 
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Fig. 2—Plane wave incident on a circular 
dielectric cylinder. 


the scattering cross-section matrix. Since this is a two 
dimensional problem, we redefine the cross section as 


| E*/? 


(39) 


o2 = Lim 2ar 


pis | Ei |? 
The source notation of the field vectors is omitted be- 
cause we are interested in the back-scattered cross sec- 
tion (echo width) and the same antenna is used for 
transmission and reception. In a manner analogous to 
that used for the three dimensional cross section it is 
shown that the echo width is related to the reaction by 


a = (4Ro)—'(wo)? | E*|-4| Cc, 2, g)|?. (40) 


In this formula source g is a distant line source, the 
source of the plane wave, and &» is the propagation con- 
stant of free space. The reactions are surface integrals 
containing only one term (since the magnetic currents 
are zero); €.g., 


G2) py = Ue, 2, ) & f f EX(g)J(o)ds 


= ff B@awas 


where the integration is over the cross section of the 
cylinder. When the total field in the scattering body is 
unknown the reaction in (41) may be approximated by 
using the methods outlined above in connection with 
(26) through (29), and with (35). 

The circular cylinder was chosen for illustrative cal- 
culations because the exact solution is known,’ and ap- 
proximations are not trivial. Calculations were made 
for four different assumptions for the field in the cylin- 
der. These are compared with the exact solution, and 
also with integral equation solutions through the second 
order.®'° It is hoped these comparisons will yield in- 
formation of use in other problems which are not soluble. 


(41) 


: r oe and - “oe a Be Differentialgleichungen und 
ntegralgleichungen der Mechanik und Physik, Mary S. R 
N. Y., p. 869; 1943. Bienes 
* W. Sternberg, “Anwendung der integralgleichungen auf beugung 
und eigensschwingungen in der electromagnetischen lichtheorie,” 
Z. Phys., vol. 64, pp. 638-649; September, 1930. 
“ D. R. Rhodes, “An Investigation of Pulsed Radar Systems for 
Model Measurements,” Rep. 475-6, Antenna Lab., The Ohio State 
University Res. Found., chap. 5; December, 1953. 
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Fig. 3—Echo width of circular cylinder vs relative dielectric constant 
for various approximations. The approximations are defined in 
Table I. 


The calculation procedure is discussed in the Ap- 
pendix, and the results are represented in Fig. 3. The 
various curves are identified in Table I. 


TABLE I 
IDENTIFICATION OF CURVES IN FIGURE 3 


Assumption for Labor of 
Curve Method the total field calcula- 
in cylinder tions 
Ab Boundary value (exact — First-order 
solution) 
Bb Zero-order integral Bae ZeTO 
equation 
Ge First-order integral eee first 
equation 
D> Second-order integral eke second 
equation ' 
E Reaction aneth ie first 
F Reaction ave tase ** first 
G Reaction ase™® 2 second 
Vee Reaction ase™® *+asge*,* second 
If Second-order integral ane e second 
equation 


8 ko, ki are the propagation constants of free space and of the 


medium in the cylinder, respectively. 
b Curves A, B, C, and D are reproduced from Rhodes work.” 


The labor of calculation increases with the order. For 
example, curves C and E are merely different combina- 
tions of the same integrals. Curve F requires the evalua- 
tion of a second double surface integral, so it requires 
about twice the labor of £. Curves G and H have 
k; in the exponential rather than ko, and the inte- 
grals have to be separately evaluated for every value 


Cohen: Application of the Reaction Concept to Scattering Problems 


197 


of dielectric constant; this is not the case for the other 
curves. The calculation for curve D includes the evalua- 
tion of a triple surface integral. In many cases, particu- 
larly where some of the integration has to be done 
numerically, this would be much more laborious than 
the evaluation of the double surface integrals for curves 
G and H. 

Curve J was obtained by using aye (the reaction de- 
termined field for curve £) instead of eo? as a starting 
function in the integral equation iteration (See Ap- 
pendix). The zero-order solution so obtained is just the 
reaction solution (curve £), and it turns out that the 
first-order solution is identical to the zero-order solution. 
The second order solution is represented by curve J in 
Fig. 3. This calculation consists of a recombination of 
the integrals used for curve D. It is seen that it gives a 
result which is farther from the exact curve (A) than 
curve D is (over the range of dielectric constants con- 
sidered). This is somewhat surprising, since I springs 
from an apparently better initial field assumption. 
Moreover, the zero- and first-order approximations (F) 
are closer to the exact curve than the corresponding ap- 
proximations with the simpler field assumption (B and 
C). However, neither B, C, nor E is at all representative 
of the exact curve, except in the region of very low 
dielectric constant. Conclusions based on these curves 
in other regions may be very misleading. 

Curve # is the only one which has the correct shape 
over the resonance region. The field assumption for this 
curve is presumably better than any of the others and 
it was expected that it would give a better approxima- 
tion than the other reaction calculations. It is also su- 
perior to curves D and J, although this could not have 
been guessed beforehand. 

The field assumption for curve H is 


ee rone erage 


(42) 


The first term corresponds to a forward-traveling plane 
wave; and the second, to a reflection from the “back 
face.” It is similar to the field inside an infinite slab; the 
difference of course is in the coefficients, which are de- 
termined by the reaction equations. They were sepa- 
rately evaluated and their amplitudes are represented 
in Figs. 4(a) and (b), where they are compared with the 
amplitudes of the internal waves in a slab one free-space 
wavelength thick, excited by a plane wave as shown. It 
is seen that the respective amplitudes in Fig. 4 are 
nearly equal. Hence the amplitudes of the reaction-de- 
termined fields are very close to what might be guessed 
at from elementary considerations. The phases, how- 
ever, are important, and cannot be arrived at by such 
means. 

Observe that there is a substantial difference between 
curves G and H, even though | s| >| as| . The explana- 
tion of this is that we have computed a back-scattering 
parameter, and the reaction between the forward wave 
(e*1) and the source is much less than the reaction be- 
tween the backward wave (e~1*) and the source. 
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Note that the coefficients in the field assumptions are 
functions of the aspect; for a given field assumption they 
have to be separately evaluated for every aspect of in- 
terest. Hence a complete calculation of scattered field 
would be very laborious. However, it would undoubt- 
edly yield interesting and valuable information. 
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Fig. 4—(a) Comparison of forward waves in cylinder and slab. 
) Comparison of backward waves in cylinder and slab. 


It is sometimes stated that variational calculations 
are insensitive to the choice of trial function and there- 
fore they give good results. Such remarks are misleading 
and must be interpreted with great caution. The curves 
in Fig. 3, for example, show a very marked dependence 
on the trial function; only one of four functions gives a 
good result over the resonance region. If the current 
assumption is the exact current, all calculation methods 
give the correct answer, and the simplest method is 
chosen. If the approximation to the current is poor then 
no method can be relied upon to give a reasonable re- 
sult, for a reasonable amount of labor. In principle, of 
course, one can get a formal convergent sequence for 
the echo width, either by successive iterations of the 
integral equation, or by using the Ritz method with a 
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complete set of trial functions. As a practical matter, 
though, calculations beyond the first few orders are 
generally prohibitive of labor. Higher order variational 
calculations are more feasible than successive itera- 
tions because the latter require increasing orders of 
integration; the former require only double volume or 
surface integrals, which have a chance of being done 
numerically. 


CONCLUSIONS 


The following conclusions have been established 
only for the particular field assumptions used here and 
for bodies whose dimensions are on the order of a wave- 
length, with relative dielectric constant less than 1.3. 
However, it is felt that they are valid over greater 
ranges and possibly for other field assumptions. They 
at least can serve as a guide in other problems, par- 
ticularly those where the exact solution is not known. 
Much more experience would be required to firmly es- 
tablish them over wide conditions. 


1. The reaction (variational) solution is superior to 
the iterated integral equation solution, for the 
same amount of labor. 

2. The approximation EF ~a;e*!*+-age—™1* is “close” 
to the correct field in the cylinder, in the sense that 
it gives an echo width which is close to the correct 
value over the first resonance region. This does 
not necessarily imply that it would give good an- 
gle-scattering results. 

3. It is very important to have a good current as- 
sumption. 


APPENDIX 
The formula for echo width (40) may be written as 


o2 = ue? (43) 
where 
“= eff E‘Eds (44) 
A 
and 
c= Lf ,3/2| E* }-*(e, — 1). (45) 


The different methods essentially consist in finding 
functions for E to be used in the integral in (44). 


Integral Equation Method 


From (3) one may write 


E= Ei+ ff T:xEds, 
A 


which is an integral equation of the second kind for ES 
the unknown z-component of total electric field. I is @ 
two-dimensional free-space Green’s Function;! it is 


(46) 


_ “P.M. Morse and H. Feshbach, “Methods of Th i - 
ics,” McGraw-Hill Book Company, Inc., New Vous Mercer 
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the electric field produced by a unit line source, all of 
whose elements are in phase and in the ¢-direction: 


T, = — fopoH (Ror), (47) 
where Hy (1) is the zero-order Hankel function repre- 
senting an outgoing wave. 


The iterated solution to (46) is found by making a 
zero-order guess for EF, called Ey. Then 


kk, = Ei'+ itn ToxEods 
A 

Es ==) fs a ‘hyp ToxFyds, 
A 


and so on. Different order approximations to o» are 
found when approximations to E are substituted into 
(44). Born approximations” result from initial choice 


Ey = eter, (49) 


This is the choice discussed by Sternberg,’ and by 
Rhodes;'® the first three approximations to the echo 
width of the circular cylinder using this initial choice 
were computed by Rhodes. His results are reproduced 
in Fig. 3as curves B, C, and D. 

A second series of calculations was made using 


(48) 


— ik 
Eo = aye"? 


(50) 


where a; was determined by the reaction calculation. 
With this assumption, it can be shown that uw)=%, so 
that the zero- and first-order approximations are the 
same, and are equal to the reaction, or variational, 
solution. Second order solution is in Fig. 3 as curve I. 


Reaction (Variational) Method 


The four assumptions for the total field used in the 
reaction calculations are given in Table I. The most 
‘involved calculation, for curve H, will now be dis- 
cussed. Let 


Sy se dee (51) 
and 
2 = e thie 
so that the field roles is 
= ashi + asho2. (52) 


2 1. I. Schiff, “Quantum Mechanics,” McGraw- Hill Book Com- 
pany, Inc., New York, INGIYS: pii59; 1949. 
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We now use (37) to find coefficients. Simultaneous solu- 


tion of (37) and (38) is not necessary because magnetic 
currents are zero. From (37), then, we have 


{Hi pye'**ds ob ff ff ul 2k (arsgr’ -- aghy’)ds'ds 
A A A 


= ff dulash: + acde)ds, w=1,2. (53) 
a 


Let 
=f f seid, w= 1,2, (54) 
A 
and 
Dy = =D, = fhe dud ds 
A ale 
— ul fafi ul 2x¢,'ds'ds, (55) 
hi psail 2 
Then the coefficients are given by 
a5 = (N Dee == N2Dy2)/(Di1D22 5 D,2”), (56) 
and 
ete (NDu = N1Do1)/(DiuD22 mee D2”). 
Thus 
Uu = c(asNy + agNo), (57) 


where as and ag are given in (56). 

The evaluation of the double surface integrals is 
rather tedious. The procedure used was to express the 
Hankel function and the exponentials as sums of Bessel 
functions. Their products are integrable and the result- 
ing series converge fairly rapidly.!° 
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Radiation Patterns of Slotted-Elliptic Cylinder Antennas” 
J. Y. WONG 


Summary—An expression is developed for the radiation pattern 
of an array of thin axial slots arbitrarily disposed around the circum- 
ference of an elliptic cylinder. Calculated azimuthal patterns are pre- 
sented for single slots on elliptic cylinders of vanishing thickness. 
The results can be applied to approximate the pattern of practical 
cylinders having very thin cross sections. 


INTRODUCTION 


HE PATTERN of slots located on cylinders of 
[[ eniotic cross section has become an increasingly 

important problem in the design of slotted- 
cylinder antennas. A paper directly applicable to this 
problem was given by Sinclair! who calculated the far 
fields of small dipole and loop antennas in the vicinity 
of elliptic cylinders. In a previous paper? an expression 
for the radiation field of axial and transverse slots in 
cylinders of elliptic cross section was derived by a 
method based on the application of diffraction theory 
and the principle of reciprocity. Attention was de- 
voted primarily to the radiative properties of trans- 
verse slot radiators. A different approach to the elliptic 
cylinder problem was extended by Wait.’ The tangen- 
tial component of the electric field in the slot was as- 
sumed to be a prescribed function and the integrals in 
the formal solution were evaluated by the saddle point 
method for the far-zone fields. 

The material contained in this paper is an extension 
of the work previously carried out for the axial slot 
problem. A general expression is developed for the 
pattern of an array of axial slots located arbitrarily 
around an elliptic cylinder. Calculated azimuthal pat- 
terns are presented for single slots and for special cases 
of dual slots located on elliptic cylinders of vanishing 
thickness. The results can be applied to approximate 
the pattern of slots in very thin cylinders such as wing 

_and tail surfaces of aircraft. 


PATTERN OF AN ARRAY OF AXIAL SLOTS 


The expression* for the radiation field in the hori- 
zontal plane, @=90°, produced by an axial slot located 
at (%o, v, 0) on the surface of an elliptic cylinder defined 
by uo is given by: 


* Original manuscript received by the PGAP, February 9, 1955; 
revised manuecript received, June 22, 1955. 

{ Radio and Electrical Engineering Division, National Research 
Council, Ottawa, Canada. 

1G, Sinclair, “The pattern of antennas located near cylinders of 
elliptical cross section,” Proc. IRE, vol. 39, pp. 660-668; June, 1951. 

2J. Y. Wong, “Radiation conductance of axial and transverse 
slots in cylinders of elliptical cross section,” Proc. IRE, vol. 41, 
pp. 1172-1177; September, 1953, 

§ J. R. Wait, “Fields Produced by an Arbitrary Slot on an Elliptic 
Cylinder,” Defence Res. Telecommun. Establishment, Ottawa, 
Canada; Rep. 19-0-8; June, 1954. 

4 Wong, op. cit., p. 1173, Eq. (7). 


Es = A >> BmSem(c, C08 0) S€m(¢, COs ¢) 


m=0 
+ Bi S0m(C, COS V)SOm(c, COS ¢), (1) 
where 
4™ 
Bie = 5 and 
Nm? Hem’ (c, cosh uo) 
Bat == oe 


N m° Hom?’ (c, cosh uo) 


The primes on the functions Hem?’ (c, cosh uo) and 
Ho,’ (c, cosh uo) denote derivatives with respect to 
the variable uw. In (1), the parameter A depends upon 
the type of excitation, but is independent of the azi- 
muthal angle @. Consequently, it is not necessary to 
know A in order to compute relative radiation pattern. 

Consider the case of an array of “p” axial slots lo- 
cated arbitrarily on the periphery of the cylinder. The 
geometry of such an array is illustrated in Fig. 1. The 


ELLIPTIC CYLINDER U, 
OF CHORD LENGTH C 


270° 


Fig. 1—Array of # axial slots. 


total field of the array can be computed by adding the 
fields contributed by each slot, taking into account the 
magnitude and phase of the exciting voltages of the vari- 
ous slots. Examination of the derivation reveals that 
the reference for the phase of each slot is at the axis of 
the elliptic cylinder, consequently it is not necessary to 
introduce any phase factor involving the distance of 
each slot from the axis. The expression for the total 
field of an array of p slots becomes 


Dp e) 
Ln 3 A xei¥t{ Bn’Sem(c, COS Vz) SEm(C, COS P) 


k=1 m=0 


+ Bn°Som(c, Cos 0x)SOm(c, CoS $) (2) 


where A; is the ratio of the voltage fed to the kth slot 
and the voltage fed to slot number 1; 7.e., A,= V;./ Vi. 


Yi is the phase angle by which voltage V; leads the 
reference voltage Vj. 
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For the particular case of two slots located on the 
major or minor axis of the cylinder, (2) can be simpli- 
fied considerably because advantage can be taken of the 
special properties of the angular Mathieu functions, 
Sem(C, cos v) and So,(c, cos v). From the equivalent 
series expansions of the angular functions, it is readily 
apparent that the following relationships exist: 


Sém(c, cos (—v)), 


— Som(c, cos (—v)). 


S€m(c, COS 2) 


Sones COs v)) = 


Values of the functions for particular angles of v are 
listed in Table I. 


TABLE I 


SPECIAL VALUES FOR THE EVEN AND Opp 
ANGULAR FUNCTIONS 
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v m S€m(C, COs 2) SOm(¢, Cos v) 
0 even 1 0 
odd 1 0 
90° even Sem(c, cos 90°) 0 
odd 0 Som(c, cos 90°) 
180° even 1 0 
odd —1 0 


Expressions for the azimuthal pattern for some spe- 
cial cases of single slots and dual slots are shown below: 


1. Single slot at v=0° (major axis), 
Es = >> BnSem(c, cos 4). (3) 
m=0 
2. Single slot at v=90° (minor axis), 


Es = >> Bom*Seom(c, cos 90°)Se2m(c, Cos ) 


m=0 
+ Bee Sore c, cos 90°) S02m+41(6, COs d). (4) 
3. Two slots at v7; =0° and v2 = 180°, 
a. Equal in-phase voltages (Vi= V2, x=0°), 
Ey = De Bom S€om(C, cos d). (5) 


m=0 


b. Equal out-of-phase voltages (Vi= V2, y= 180"), 

Ey = x Bom+1%S€2m+1(C, COS @). (6) 
m=0 
4. Two slots at 7: =90° and v2= 270°, 
a. Equal in-phase voltages, 
Es, = so Bam’S€am(C, COS 90°)S€am(C, COS $). (7) 

m=0 

b. Equal out-of-phase voltages, 

E, = 3 Bomy19S0am+1(6, COS 90°)S'am41(6, CoS @). (8) 


m=0 
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DISCUSSION OF RESULTS 


In order to facilitate the computations of the radia- 
tion patterns, only the case of a degenerate elliptic 
cylinder (uw —0) is considered. Simplification of the 
expressions result, because advantage can be taken of 
the joining factors associated with the various radial 
and angular Mathieu functions.® Patterns are com- 
puted for three different sizes of cylinder. The values 
chosen for the major axis or chord dimension C, are 
0.55, 1.0, and 1.42d. The results of these calculations 
are shown in Figs. 2, 3, 4, 5, and 6. The elliptic cylinder 
is oriented so that the major axis coincides with the 
ordinate axis of the coordinate system (0°-180° line). 

Eq. (3) represents the pattern of a single axial slot 
located on the major axis of the cylinder. It is not dif- 
ficult to demonstrate that (3) yields a circular pat- 
tern for this particular case. Simplification of (3) re- 
sults by making use of the following identities :® 


1 


Hen,’ (¢, 1) = — ———— 9 
(ert) Tey. (9) 
and 
egic(cosh u cos v cos @) — (Sar) 1/2 Se Teme, cosh 2) 
m=0 mo 
‘S€m(C, COS V)S€m(6, COS d). (10) 


Substituting (9) and (10) in (3) and simplifying, yields 
the following expression: 


eilka cos ¢—(m/2)) 


(11) 


POS cei 
Eq. (11) can be recognized as the equation of a circle, 
hence the pattern of a slot located on the edge of a de- 
generate elliptic cylinder is circular. Furthermore the 
pattern is independent of the size of cylinder. This re- 
sult could have been inferred from physical reasoning. 
A magnetic dipole located parallel to the edge of a de- 
generate elliptic cylinder does not excite any current on 
the cylinder, consequently the presence of the cylinder 
can in no way affect the radiation. 

Patterns of a single axial slot located on the minor 
axis are illustrated in Fig. 2. Most of the energy is con- 
centrated in a single lobe in the forward region, al- 
though some diffraction of the energy occurs around the 
edges of the cylinder. 

The effect of offsetting the slot from the center line 
for the case C/A=0.55 is shown clearly in Fig. 3. The 
parameter x refers to the amount of offset and a refers 
to one-half the chord dimension. As the slot approaches 
the edge of the cylinder, two effects occur. The pattern 


5 “Tables Relating to Mathieu Functions,” prepared by the Com- 
putation Lab. of the Nat. Appl. Math. Lab. Nat. Bur. Standards, 
Columbia University Press, New York, N. Y.; 1951. 

6 P. M. Morse and H. Feshbach, “Methods of Theoretical Phys- 
ics,” Part II, McGraw-Hill Book Co., Inc., New York, N. Y., pp. 
1572-1573; 1953. 
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Fig. 2—Pattern of single slot located at v=90°. 
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Fig. 3—Effect of offsetting slot from the minor axis. (a) Pattern of slot located at v=80° (x=0.17a). (b) Pattern of slot 


located at v=60° (x =0.50a). (c) Pattern of slot located at v=40° (x =0.76a). 
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Fig. 5—Pattern of two slots located at 1; =0° and v?=180° fed with equal out-of-phase voltages. 
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Fig. 6—Pattern of two slots located at 71=90° and v2 =270° fed with equal out-of-phase voltages. 


of the main lobe becomes more asymmetrical with the 
maximum directed away from the side containing the 
slot, and secondly more radiation is encountered in the 
shadow region. 

The patterns of some special cases of dual slots are 
shown in Figs. 4, 5, and 6. Applying the results ob- 
tained for a single slot an an edge, it is apparent that 
for the case of a pair of axial slots located on the major 
axis of the cylinder, the patterns are the same as two 


free-space dipoles placed at a distance C apart, either 
in phase or 180° out-of-phase. 


It was found that the case of two slots located on the 
minor axis excited with equal in-phase voltages yielded 
an omnidirectional pattern. If one examines the electric 
field lines produced by a symmetrical balanced antenna 
located on either side of a thin conducting plane, it can 
be seen that the fields at the edges combine in such a 
manner as to eliminate any diffraction effects. 
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The Transmission-Line Properties of a Round 


Wire between Parallel Planes* 
H. A. WHEELER} 


Summary—A transmission line is useful as a primary standard of 
impedance if there is a formula relating its wave resistance to its 
cross sectional shape. For this purpose, the round wire between 
parallel planes is computed to a very close approximation. The 
_ method is based on dividing the usual single pole into a pair of spaced 

half-poles, between parallel lines on the plane of cross section. The 
resulting formula covers all shape ratios to a high degree of approxi- 
-mation. Simple explicit formulas are given for the shape ratio and the 
' wave resistance, each in terms of the other, sufficiently close for 
_-practical purposes. For example, a line made of perfect conductors 
(assumed) and air dielectric (laboratory conditions), has a wave re- 
“sistance of 50.000+.003 ohms if the ratio of radii is 0.54850 
=1/1.82315. 


* Original manuscript received, February 18, 1955, revised manu- 
script received, May 20, 1955. This paper presents a small part of the 
treatment of the same subject in Wheeler Monographs, No. 19, 
June, 1954. Reprints of this complete paper are available from the 
_ .author on request. 

+ Wheeler Laboratories, Great Neck, Noy: 


I. INTRODUCTION 
Meso atie transmission lines of certain shapes of 


cross section are useful as primary standards of 

impedance because it is possible to compute the 
relation between their wave resistance and their cross- 
sectional dimensions. The simplest case is the coaxial 
line whose inner and outer conducting surfaces have 
circular cross section. 

Another shape of coaxial line has come into use be- 
cause it offers some advantages of mechanical con- 
struction and of accessibility for standing-wave de- 
tection. This shape is the round wire between parallel 
planes [4, 5, 13, 15, 20, 21]. It offers a problem in ap- 
proximation because the relation between its wave re- 
sistance and its cross-sectional shape ratio is not ex- 
pressed exactly by any of the usual tabulated functions. 
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The formulation and computation of this relation is 
the subject of this paper. 

The method of computation on the plane of cross sec- 
tion is a change from the usual single pole between paral- 
lel lines, to a pair of spaced half-poles giving an equi- 
potential line closely approximating circle represent- 
ing the inner conducting surface. This method places 
the wave resistance between upper and lower limits 
which are close together for the usual cases. The closest 
approximation is a mean value between these limits, 
which is found to be closer than required for practical 
applications. This complete formulation requires elab- 
orate computation, though relying only on hyperbolic 
and more elementary functions. 

For practical purposes, a simple series is derived for 
explicit computation of either the shape ratio or the 
wave resistance in terms of the other. Each series con- 
verges so rapidly that the first two terms are sufficient 
for the usual requirements. 

For use as a primary standard of impedance, the 
50-ohm line is computed to 5 places. The correction is 
applied for air dielectric in laboratory conditions. The 
remaining corrections are those for the skin effect in 
the conductors and for the curtailment of the parallel 
planes down to finite width; both are given and are 
shown to be small. 

The present formulas are shown to be a great ad- 
vance over the formula of Wholey and Eldred [13] 
which is the most advanced solution previously pub- 
lished. 

Only the simplest and most useful formulas are pre- 
sented herein. The complete treatment has been given 
by the writer [21]. 


II. SymBoLs 


ry =radius of round wire. 
c=distance from center of wire to each plane. 
c’ = (4/m)c=radius of round shield equivalent (for 
small wire) to parallel planes. 
R=wave resistance of line (ohms). 
R,=120 wr (approximate) =wave resistance of free 
space (ohms). 
R,/2m = 60 ohms (approximate). 
a=offset of half-pole. 
p= Geis) ee oer cl. 


ters (radians). 


a=(m1/4) a/c=a/c'; parame- 


III. THe BAsIs FoR THE FORMULAS 


Fig. 1 shows the cross section of a round wire be- 
tween two planes. The outer contour of the wire is the 
circle of radius ry and the inner contour of the planes is 
the pair of parallel straight lines at a distance c from 
the center of the circle. Perfect conducting material 
(free of resistance) is assumed for all contours. The 
wave resistance of the space between the conductors 
is taken to be 120 7 ohms, which is approximately that 
of free space or air dielectric. 


The outer circle, of radius c’=(4/m)c, is the circular 


IRE TRANSACTIONS—ANTENNAS AND PROPAGATION 


October 


outer contour which is equivalent to the parallel planes 
if the wire is small (r/c « 1) [15]. To simplify some of the 
formulas, they are expressed in terms of the ratio, 


which approaches the limiting value 7/4 as the wire 
approaches contact with the planes. 

The formulas to be presented are based on a refine- 
ment beyond the usual formula for a small wire be- 
tween parallel planes. This refinement preserves the 
parallel planes and enables closer computation of the 
potential on the wire. The usual single pole in the center 
of the wire is divided into a pair of half-poles located at 
a distance of a on opposite sides of the center as shown 
in Fig. 1. 


_—— | —_- 


! 
‘PAIR OF | 
HALF-POLES | 


Fig. 1—Cross section of round wire between parallel planes, 
approximately to scale for R=50 ohms. 


If the circular contour were replaced by the pair of 
half-poles, the potential on the locus of the circle would 
have a cyclic variation around the circumference. The 
distance of the half-poles from the center is chosen so 
that the potential on the circle has equal maximum 
values on four points (the two points nearest to the 
planes and the two points farthest from the planes). 
This gives nearly uniform potential around the circle. 

By the usual methods, it is then possible to compute 
the maximum potential at these four points and the min- 
imum potential at intermediate points on the four 
quadrants of the circle. A circular conductor at this 
locus would assume a constant potential having a value 
which is a mean between the computed maximum and 
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minimum values. Since these values are close together, 
the mean value gives an extremely close approximation. 
This rather elaborate computation is formulated and 
exemplified in Reference 21; it will not be given here 
but some of its results will be utilized. 

The evaluation of the proper distance to each half- 
pole has been reduced to a simple formula by consider- 
able effort. It is given here, after which the remaining 
derivations become routine. It is expressed in terms of 
the angular parameter, 


a 1 
== a anticos +/cosh 2p cos 2p 
¢ 


= 3 antisin 2\/sinh? p sin? p — (sinh? p — sin? p) 
= V2/3 p> +---. (1) 
IV. FORMULAS 
Simple Logarithmic Formula 


The simple logarithmic formula for a small wire be- 
tween two planes has been known for many years [4, 5, 
6, 9, 10, 11]. It is the first term in a series expansion for 
a small wire. By the present procedure (mean value), 
the second term has been added. The resulting for- 
mula is 


1 2 
R= 60(In— — = pt — 
p 9 


) ohms. (2) 


The coefficient of the second term is the mean of 23/90 
and 17/90 for the maximum and minimum values of 
potential; these values are only +1/30 from the mean 
value 2/9. 


New Logarithmic Formula 


The best simple formula that has been derived is the 
. following, an outcome of the present study. 


1/1 + 2 sinh? p + /1 — 2 sin? p 
4/2(sinh? p + sin? p) 


R= 60 E 


1 
zed 014p? + --- | ohms. (3) 


The main term in this formula is based on the maxi- 
mum potential on the locus of the circle. It alone is 
closer than any previously published formula, for all 
shapes from a small wire to a large wire near contact, 
Even in latter extreme case, its relative error is less than 
0.12 (as compared with nearly unity for Reference 13) 
' The second term in (3) is the correction to obtain 
‘the mean value for small wires (as mentioned above). 
The third term is a further correction to obtain the 
mean value for a rather large wire (20-ohm line). Any 
smaller wire has a remaining correction much less than 
the third term. For a 50-ohm line, (2) is within 1/4000 
and (3) is within 1/100,000 (the third term being about 
-1/50,000), whereas the closest previous formula is off 
by 1/300 [13]. 
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Explicit Formula for Shape 


The most remarkable outcome of this derivation is a 
simple explicit formula for the ratio of radii in terms of 
the wave resistance. This formula is in the form of a se- 
ries which, for a small wire, converges so rapidly that the 
first two terms meet all the usual needs. It gives a close 
approximation for all shape ratios, even near contact. 


Yr 4 asi 2 
= — exp (1 


ea etl: fs te 8 Rane 
(3:4) 0) tg (S = 
ee CEeeaG 


ce 60 9 
mK 
eo ULoozeexy 
riveree ass oe ) (4) 
aire! ex aie Shae ; 
ae 


In the parentheses, the coefficients of the third, fourth 
and fifth teriis are approximated on a logical basis in 
Reference 21; the remaining error is less than these 
terms. For example, if R=50 ohms, the first two terms 
leave a relative error of about 1/300,000. The first two 
terms are more than adequate for slide-rule computa- 
tions and all practical purposes if R>20 ohms, leaving 
a relative error less than 1/5000. 


Precise Properties of Free Space 


Very close computation requires that the wave resist- 
ance (R.) be evaluated closer than the assumed value 
(1207 =60 X27). This requires small corrections for 
free space and air. According to the latest authority 
[18], the wave velocity in free space is 


299.793 = 300 — 0.207 = 300/1.00069 meters/ysec. (5) 


Since the magnetic constant in free space is exactly 
0.4 + microhenries per meter, this correction of velocity 
is caused by electric constant being slightly greater than 
that corresponding to 300 meters per psec. 


Precise Properties of Air Dielectric 


The dielectric constant in air is slightly greater than 
that in free space, in proportion to the density of dry 
air and, in greater degree, to the density of water vapor. 
These factors have been formulated and a standard at- 
mosphere has been specified as typical of laboratory 
conditions [7]. This atmosphere is defined: 


Condition Standard 
Temperature ==025 -C 
Relative humidity = 0.60 
Altitude = 0 (sea level) 
Total pressure = 1 atmosphere 


For present purposes, some limits of variation are stated: 


Condition Limits 

Temperature °C 

Relative humidity +0.20 

Pressure variation with time +1/60 atmosphere 

Pressure variation with —1/25 atmos (+1000 
altitude feet) 
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The standard conditions, modified to the extent of 
these variations, give a dielectric constant of 


1.00070 + .00010. (6) 


Dividing the free-space velocity (5) by the square root 
of this dielectric constant, 


300/(1.00104 + .00005) = 300 — 0.31 + .02 


= 299.69 + .02 meters/usec. (7) 


This value is taken to be typical of laboratory space in 
regions of low altitude. 

The wave resistance in this air dielectric is 0.4 7 
times the velocity: 


R, = 120r/(1.00104 + .00005) 


= 376.60 + .02 ohms, (8) 
R,/2e = 60/(1.00104 + .00005) 
= 59.938 + .003 ohms. (9) 


For very close computation, the latter value is used 
instead of 60 ohms in various formulas such as (2) to 
(4). The tolerances of variations in the atmosphere leave 
an uncertainty of +1/20,000 of the wave resistance of 
the transmission line. For a 50-ohm line, this is only 
+ .0025 ohms, which is negligible in practice. 


Skin Effect 


A transmission line to be used as a standard of im- 
pedance is usually in the form of a slotted line for stand- 
ing-wave observations. In this application, its length 
should be at least one-half wavelength and preferably 
one wavelength. Since a slotted length longer than one 
meter is inconvenient for laboratory use, the slotted 
line is most useful at wavelengths less than 2 meters, 
or frequencies above 150 mc. In this case, the skin depth 
in copper or silverplate is less than 5.4 microns or 0.21 
mil-inch [17]. After computing the dimensions for a 
perfect conductor, each surface exposed to the field in 
the line should be raised by + the skin depth. Such a 
correction can be computed for any one frequency, or a 
mean value taken for a frequency band. However, it is 
so small that it is usually negligible in any of the non- 
magnetic metals used for structure or plating, such as 
silver, copper, gold, aluminum, or brass. 


Fimite Parallel Planes 


The computations are based on parallel planes ex- 
tending indefinitely far from the wire between them. 
In practical construction, the use of a finite width of 
both planes introduces some uncertainty in the wave 
resistance. For any width, the edges may be closed 
(except for a slot) to form a shield of rectangular cross 
section, in which case the change of wave resistance 
can be computed for a small wire (by the method of 
images, from the formulas in References 4, 5, 11). This 
change is tabulated for various ratios of the width of 
each plane, over the separation (2c) between planes. 
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Width Fraction Reduction of 


Separation Wave Resistance 
2 0.45 ohm 
2.5 .094 
3 .020 
3.5 .0041 
4 .00084 


From this table, the width of the planes can be made 
sufficient to subordinate this correction. Leaving the 
edges open has the opposite effect, increasing the resist- 
ance, but less than 2 as much. However, the open edges 
leave some uncertainty depending on the external con- 
tours and environment. 


End Connections and Supports 


One end of the slotted line should be precisely matched 
to an extension line of equal wave resistance, usually of 
coaxial circular cross section, which connects with the 
impedance under test. The matching of this connection 
is not adapted for computation, but can be accomplished 
by experimental adjustment of the contours. If the line 
is long enough to need supports at both ends, and per- 
haps also distributed along its length, there remains a 
difficult problem of providing a dielectric support in a 
form that is reliable and computable. 


Computation of 50-ohm Line 


A 50-ohm line is of particular interest for use as a 
slotted line for standing-wave tests of reflection, so it 
has been computed closer than the variation of air di- 
electric in laboratory conditions. From (4), the approxi- 
mate value of the shape ratio is r/c =0.54896; this has 
been verified by mean-value computation and (3). 
Applying the corrections for free space and air dielec- 
tric, the final computation gives 


r/¢ = 0.54850 = 1/1.82315 (10) 


for 


R = 50.000 + .003 ohms. 


Comparison with Published Curves 


Referring to the writer’s previously published curves 
[15], the present derivation indicates that, for the round 
wire between parallel planes, the previous curve is about 
1 ohm too low in the transition region of 20 to 50 ohms. 


V. CONCLUSIONS 


The round wire between parallel planes has yielded 
to a simple and direct computation of its properties, to 
a degree of approximation more than adequate for the 
usual applications. Simple explicit formulas are given 
for the shape ratio and the wave resistance, each in 
terms of the other. The complete interpolation formula 
covers the entire range of shape ratio and wave resist- 
ance. 
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Current Distribution on Wing-Cap and Tail-Cap Antennas* 
IRENE CARSWELL} 


Summary—Surface current distributions on aircraft are examined 


for wing-cap and tail-cap antennas in the hf band using scale models. 


Current paths involved in structural resonances are investigated, and 


the effect of such resonances on antenna impedance is discussed. 


Measured results are also applied to the interpretation of radiation 


patterns. 


A technique is described for measuring the relative amplitude and 
phase of the surface current, and for determining the direction of 
current flow. In this method the detector is placed in the model and 


the probe is driven by an external source of rf power. This technique 
permits the frequency to be changed easily. The system is more 


sensitive than others described in the literature, since a higher power 
source can be used when it is not required to be enclosed in the air- 


craft model. 


INTRODUCTION 
[: THE hf frequency range, where the principal 


dimensions of an aircraft are comparable to the 
wavelength, current paths composed of portions 


' * Original manuscript received by the PGAP January 10, 1955; 
revised manuscript received June 22,1955. 
+ Stanford Research Inst., Menlo Park, Calif. 


of the fuselage, wings, and tail surfaces are known to be 
responsible for resonance effects in the input impedance 
curves of wing-cap and tail-cap antennas.’ Frequencies 
at which the principal current paths on an aircraft may 
become resonant have been predicted by a backscat- 
tering method.? Determination of the nature of these 
resonances, and verification of the frequencies at which 
they occur, is of interest in connection with the excita- 
tion of an aircraft as an antenna, since both the input 
impedance and radiation patterns are affected.* Cur- 
rent distribution is a fundamental electrical property for 
investigation, since the other quantities are dependent 
on it in the usual antenna analysis procedures. 

1 J. V. N. Granger, “Wing-Cap and Tail-Cap Aircraft Antennas,” 
Tech. Rep. 6, Contract AF 19(122)-78; SRI Proj. 188, Stanford Res. 
Inst., Menlo Park, Calif.; March, 1950. 

2 A. S. Dunbar, “Electromagnetic Resonance Phenomena in Air- 
craft Structures,” Tech. Rep. 8, Contract AF 19(122)-78; SRI Project 
188, Stanford Res. Inst., Menlo Park, Calif.; May, 1950. 

3 J. V. N. Granger and T. Morita, “Radio frequency current dis- 


tributions on aircraft structures,” Proc. IRE, vol. 39, pp. 932-938; 
August, 1951. 
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When an antenna is mounted on the surface of an 
aircraft structure, rf currents are not confined to the 
antenna alone, but flow on the conducting surface of 
the aircraft as well. If the relative amplitude and phase 
of the surface current could be determined experimen- 
tally, it would be possible to locate resonant path-lengths, 
and to discuss, at least qualitatively, their effect on 
input impedance and radiation patterns. 

This paper describes a measuring technique for ob- 
taining the distribution of surface current density. 
The method was applied to determine the current dis- 
tributions on a 749 Constellation aircraft. All measure- 
ments contained in this report were made on a 1/43- 
scale model aircraft, constructed of cast aluminum. 
The measured distributions are given and the relation 
of these distributions to the measured input impedance 
and radiation patterns is discussed. 
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Fig. 1—Orientation and response of loop probe. 


METHOD oF MEASUREMENT 
Loop Probe 


If the conductivity of the model is high, and skin 
depth is small, it is known from Maxwell’s equations 
that the current vector and magnetic field vector are 
perpendicular to each other at the conducting surface, 
and mutually perpendicular to a vector directed out- 
wardly normal to the current-carrying surface. This re- 
lationship is expressed by the boundary condition, 


AxXxH=—J, 
where 7 is the unit vector normal to the current carry- 
>} 


ati 
ing surface, H is the magnetic vector, and J is the sur- 
face current density. Therefore a loop probe held per- 
pendicular to the surface of the model will indicate the 
magnitude and direction of surface current density by 
virtue of the proportionality of its output to the num- 
ber of linkages with the magnetic flux.4 The ideal re- 
sponse of the loop, as illustrated in Fig. 1, is a figure 
eight having maxima when the loop plane coincides with 
the direction of the surface current at a given point; 
nulls when the loop plane parallels the magnetic vector. 


4J.V. N. Granger and T. Morita, “Current Distribution on Air- 
craft,” Tech. Rep. 78, Cruft Lab., Harvard Univ., Cambridge, Mass.; 
November, 1949, 
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Ideally, the shielded-loop probe should not respond 
to the codirectional electric field, or “dipole mode.” 
Practically, however, some response to this mode will 
exist, due to slight mechanical dissymmetries in the 
structure. Response of the loop probe to electric field 
pickup is indicated by the variation of the output of 
the measuring system from the desired cosine form, 
which is illustrated in Fig. 1. 

When the electric field has components that are not 
parallel to the line of probe symmetry, such as occur 
near cusps in the conducting surface, and near regions 
of impressed field, errors in the measurement will be 
unavoidable. This fact must be considered in inter- 
preting the data. 


Amplitude Measurement 


The measuring technique utilized in this investiga- 
tion differs somewhat from others described in the liter- 
ature. Instead of driving the antenna terminals and 
measuring the current picked up in the loop probe, the 
loop is driven and the current in the antenna terminals 
is measured. By locating the detecting equipment in- 
side the aircraft model, as illustrated in Fig. 2, and 
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Fig. 2—Block diagram of apparatus for determining 
amplitude of current. 
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bringing out the detected audio component on high- 
resistance leads, the rf field around the model is not 
disturbed.® This audio signal is amplified and presented 
to the observer on a vacuum-tube voltmeter. The 
source of power is an rf oscillator modulated at 1,000 
cps, whose output is carried to the loop probe by double- 
shielded coaxial cable. The direction of the loop plane 
at maximum output coincides with the direction of 
current flow, and the amplitude of maximum at that 
point is relative amplitude of surface current density 
with respect to an arbitrarily chosen reference. 
Elliptically-polarized currents can occur on portions 
of the surface where the lines of current flow change 
direction during the rf cycle. Their presence is indicated 
by a “filling in” of the nulls in the output as the loop 


*R. M. Hatch, Jr., “Improvements in Instrumentation for the 
Investigation of Aircraft Antenna Radiation Patterns by Means of 
Scale Models,” Tech. Rep. 26, Contract AF 19(604)-266, SRI Project 
591, Stanford Res. Inst., Menlo Park, Calif.; August, 1952, 
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probe is rotated. Regions of elliptical polarization are 
shown in the current plots as areas where no current 
distribution data was obtained. 


Phase Measurement 


, The phase of the surface current density is measured 
by comparing the signal from the probe with the signal 
received from a dipole illuminating the model and con- 
nected to the same source. Reference to Fig. 3 shows 
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Fig. 3—Block diagram of apparatus for determining 
relative phase of surface current density. 
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that the dipole is mounted in such a way that its dis- 
tance from the model is variable. Since the phase of the 
signal in the far field is a known linear function of the 
distance from the antenna under test, the dipole is an 
rf phase shifter. When this signal is combined with the 
signal from the loop probe and impressed on a bolom- 
eter detector, the output is presented to the observer 
as a succession of maxima and minima as the dipole 
is moved. It is important that the amplitudes of the 
two signals be of the same order of magnitude, in order 
for the minima to be well defined. It is convenient to 
define a position of the dipole at which a sharp null is 
observed when the loop is in a reference position on 
the aircraft model, as zero relative phase. The relative 
phase of the current at any other point on the model 
is then computed from the distance in wavelengths 
that the dipole must be moved to produce a new mini- 
mum voltage. 


MEASURED RESULTS 


Fig. 4 shows a contour plot of the amplitude of the 
current distribution on a 749 Constellation aircraft 
excited by a wing-cap antenna at a full-scale frequency 
of 11.2 mc. At this frequency a pathlength L composed 
of the leading edge of the driven wing and the side of 
the fuselage forward of the wing, is one wavelength 
long. Examination of Fig. 4 shows that the level of cur- 


Fig. 4—Contour plot of amplitude of surface current on constellation aircraft excited by a wing-cap antenna at a full-scale frequency of 11.2 me. 
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rent amplitude along this path is high when compared 
to relative amplitude elsewhere on the model surface. 
Reference to the input impedance curve of Fig. 5° re- 
veals a sharp peak in the curve of input resistance at 
11.2 mc.? Resonance of the wing and fuselage current 
path appears to be responsible for this effect. 

A contour plot of relative amplitude at 7.1 mc is 
shown in Fig. 6 for the same antenna. Although the 
wingspan is approximately a wavelength long at this 
frequency, no strong resonance is apparent, either in the 
plot of current distribution or in the input impedance 
curve of Fig. 5. Inspection of Fig. 6 reveals that undriven 
wing is only weakly coupled to radiating structure. Cur- 
rents flowing on driven wing tend to divide principally 
between fore and aft portions of fuselage where no reso- 
nant pathlength exists at this frequency. 

In the vicinity of 24 mc, another resonant peak occurs 
in the resistance of the wing-cap antenna. Fig. 7 shows 
the distribution of the current amplitude at this fre- 
quency. It should be noted that at 24 mc, the wing chord 
is no longer small compared with the span and the 
wavelength, and the diameter of the fuselage is an ap- 
preciable fraction of the wavelength. The quasi-static 

§ Presented here through the courtesy of the Lockheed Aircraft 
Corp., Burbank, Calif. 

‘J. Taylor, “Flush-Mounted H-F Antennas for the 1049 Con- 


stellation Aircraft,” Final Report, SRI Project 409, Stanford Res. 
Inst., Menlo Park, Calif.; February, 1952. 
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Fig. 5—Impedance of cap-type antenna on wing tip of 
constellation aircraft. 
assumptions of longitudinal current flow and static 
amplitude distribution at any cross section, which are 
valid at lower frequencies, are not applicable here.® 
At 24 mc, the wingspan is nearly three wavelengths long, 
as is a path involving the driven wing, the leading edge 
of the horizontal stabilizer, and the slide of the fuselage 
aft of the wing root. At the junction of the wing and 
fuselage, the wing chord is nearly a half-wavelength 


8 J. V. N. Granger and T. Morita, op. cit. 


Fig. 6—Contour plot of amplitude of surface current on a constellation aircraft exited by a 
wing-cap antenna at a full-scale frequency of 7.1 mc. 


CURRENT FLOW LINE 


— —— CONTOUR OF CONSTANT 
AMPLITUDE 


LD» 


— 20 


Aa 68 ae Z 
Tae —+~ 7 
5 


ieee = 


Fig. 7—Contour plot of amplitude of surface current on a constellation aircraft excited by a 
wing-cap antenna at a full-scale frequency of 24 mc. 
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Fig. 8—Contour plot of surface current amplitude and phase on a constellation aircraft excited by a 
tail-cap antenna on the center fin at a full-scale frequency of 15 me. 
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long, and the chord of the isolated wing-tip is a quarter- 
wavelength long at its base, where the driving voltage 
is impressed. 

Amplitude and phase of surface current on Constella- 
tion with center fin isolated are in Fig. 8 on the previous 
page. Since system possesses symmetry about center 
line of fuselage, only half of each contour is plotted. 
The measurements were made at 15 mc, where the 
span of the horizontal stabilizer is a half-wavelength, 
and where a peak in the input resistance is known to 
exist. Inspection of the amplitude distribution reveals 
that the currents on the horizontal stabilizer are high, 
whereas the coupling to the rest of the structure is low. 
The phase angle increases with distance from the 
source, as would be the case for a lossy transmission 
line. It lags with respect to the antenna terminals, and 
is not constant at the cross sections. 
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Fig. 9—Field pattern of a tail-cap antenna on center fin of a 
lockheed constellation at 15 mc. 


The field pattern of a tail-cap antenna on the center 
fin of the Constellation is illustrated in Fig. 9, where 
the standard spherical coordinate system has been used 
to define the principal planes. The azimuth angle ¢ is 
measured in the counterclockwise direction from the 
nose of the model. The angle of elevation @ is measured 
from a point directly above the model and is 90 degrees 
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at the horizontal plane. Some of the ways in which the 
measured current distribution at 15 mc affects the 
radiation pattern at this frequency may be pointed out. 
It must be emphasized, however, that exact interpre- 
tations must await the completion of further work. 
The transverse-plane Hs pattern appears to be the re- 
sult of the current on the center fin, modified by the 
effect of current excited on the horizontal stabilizer. 
The longitudinal pattern is approximately that of a 
bent dipole consisting of the center fin plus the top of 
the fuselage aft of the wing, modified by the contribu- 
tion of currents on the sides and bottom of the fuselage, 
which are lagging in phase, and are of lower amplitude. 
It should be noted that the currents on the fuselage for- 
ward of the wing are too low to affect the major lobes 
of the pattern. The £, azimuth pattern exhibits nulls 
fore and aft, and off the wing tips, as would be expected. 
It approximates the radiation from a horizontal dipole 
about one wavelength long, with modifications due to 
the currents on the statilizer and wings. 


SUMMARY 


This paper has described a technique for measuring 
the relative amplitude and phase of the surface current 
excited on the skin of an aircraft by antennas in the 
hf range. Measured results have been obtained at sev- 
eral frequencies for wing-cap and tail-cap antennas, 
using scale models of the aircraft, and proportionately 
increased operating frequencies. An attempt has been 
made to establish relationships between the current dis- 
tribution data and the appropriate radiation. pattern 
and input impedance information, at least to a first ap- 
proximation. Further work will be necessary before 
more exact interpretations are possible. 
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Note on a Method for Calculating Coupling Coefficients 
of Elements in Antenna Arrays’ 
V. T. NORWOODt 


N EXPLICIT expression for calculating cou- 
A pling coefficients for each element of a traveling- 
wave array is derived below. The coupling 
coefficients are shown to be completely determined by 
the desired current distribution across the aperture, 
the waveguide attenuation, and the percentage of the 
total power which is jettisoned into the load. The ad- 
vantages of the explicit expression over the difference 
equation method previously used are (1) facility in de- 
termining the proper balance between the amount of 
power assigned to the terminal load and the maximum 
permissible coupling between an element and the 
array, and (2) computing simplification; i.e., a com- 
_ puting error affects only one coefficient without causing 
errors in every successive coefficient. The importance 
of these points is dramatized by considering the work 
involved in the design of a 1,000-element Dolph-Tche- 
byscheff array. 
The value of the power at various points in the array 
is indicated in Fig. 1. The subscript in the first position 
refers to the power into a kth element (Pix), radiated 


normalized to 1 for the last element. We shall call these 
coefficients C,. 


ity ae Or ate (A) 
Combining (2) and (3) yields 
Pity = "(Pa Pe * (S) 
or 
= (Pin — CrPenlem@. 


Starting with P;;=1 we may deduce the kth term 
after writing a few terms. 


IP = 1 
Pi = (1 — Cure. 


1D fel — Ci) Cae — CoP rn |e-% 
Le= | [e-2a — CyP pre?! — CoPrne~@ |CsP en} e- 
k 
Piceps) a lage De Cp SG awe (6) 


d=1 


Fig. 1—Schematic representation of power in a slot array. 


by a kth element (P,x), or just beyond a kth element 


(Pox). 
Assume that 


Pj; = 1 (total power in). (1) 

The following relationships can be readily deduced: 
Pow = e'Picny1y (a attenuation in guide) (2) 
Pix = Pow + Pre: (3) 


The power radiated from each slot is proportional to 
the desired aperture distribution (in terms of power), 


* Original manuscript received May 3, 1955. | 
if ice and Decco Labs., Hughes Aircraft Co., Culver 


City, Calif. 


Eq. (6) is then written for k=n, and the expression 
P.in—P 1 is substituted for P,,. The resulting expression 
can be solved to show the relationship between P,, and 
P, (the power into the load). (Alternatively everything 
could have been written in terms of Pp.) 


n 
Led) ag tes [aps bab Ca Oe (7) 
d=1 
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Pp = — (8) 
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We may now write the coupling coefficient 


Cal 
cs (9) 


k 
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Prk 
Por 


which can be expressed in terms of /, Pz, eyousl ox, Witt 
should be noted that the summation terms are inde- 
pendent of load considerations and these should be 
computed first in the design of an actual array. With 
these values at hand a trial value for Pz, may be sub- 
stituted into the expression of y and a little practice 
quickly enables one to determine the maximum con- 
ductance value by computing conductances for about 
three consecutive elements. Two or three values for 
P, usually result in appropriate maximum conductance, 
set as a design criterion. An analytic method for solving 
this extremal problem is easily derived, but expressions 
are unwieldy except for a few trivial aperture distribu- 
tions. After 7, Pz, and a are determined, P,, becomes a 
constant and computing is straightforward. 
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Fig. 2—Equivalent circuit of a slot in an array. 


Fig. 2 shows an equivalent circuit of an element and 
the lumped conductance of the remainder of the array 
(for example, a waveguide slot array). By conventional 
circuit relationships, 


Peg Gi. 


; (10) 
Por , 


For G, equal to unity (a common assumption in the 
design of linear arrays) the ratio P,«/Po. reduces to the 
conductance of the element as it is measured in a trans- 
mission line. 


On-Axis Defocus Characteristics of the 
Paraboloidal Reflector* 


D. K. CHENG} anp S. T. MOSELEYt 


Summary—Lack of an adequate test site often requires the use of 
on-axis defocusing to obtain simulated Fraunhofer radiation patterns 
in the optical-Fresnel region from microwave reflector-type antennas. 
This note presents a way of determining the required amount of de- 
focus from the viewpoint of aperture phase correction. It is shown 
that the important part of the main lobe of the radiation pattern can 
be approximately reproduced in the optical-Fresnel region with proper 
amount of defocus. For a given amplitude illumination function, the 
change in beamwidth with the position of the primary source can 
also be predicted. 


INTRODUCTION 


HE TREND toward higher operating frequen- 
cies and larger aperture dimensions for micro- 
‘ wave reflector-type antennas has brought about 
the need for taking far-zone (Fraunhofer region) radia- 
tion patterns at distances far too great for the average 
available test site. This has led to a defocusing practice 
whereby the antenna feed is displaced along the prin- 
cipal axis away from the reflector in the attempt to 
simulate far-zone patterns in the quasi-near zone (op- 
* Original manuscript received by the PGAP, March 16, 1955. 
The work described in this note was done under Contracts AF 
30(602)-925 and AF 30(602)-1360 sponsored by the Rome Air Dey. 
Center, USAF. Presented at the IRE Western Convention, Los 
Angeles, Calif., August 25, 1954. 
+ Elec. Engrg. Dept., Syracuse University, Syracuse, N. Y. 


t Signal Corps, U. S. Army Pacific Headquarters, Honolulu, 
Hawaii. 


tical-Fresnel region). The amount of defocus required 
for a given test site is usually determined somewhat 
empirically by an approximate geometrical optics 
method. 

This note presents an analytical method of deter- 
mining the required amount of defocus from the view- 
point of aperture phase correction for the circular aper- 
ture diffraction problem. Proper defocusing calls for the 
cancellation of the effect of the Fresnel phase term in 
the diffraction integral by a phase term as a result of 
defocusing. The required amount of defocus is also 
determined from the criterion that the beamwidth of 
the main lobe in the simulated pattern be a minimum. 
It will be shown that the minimum beamwidth criterion 
coincides with the condition for maximum radiation 
amplitude along the reflector axis. 


THE DeEFocus INTEGRAL 


If it is assumed that the aperture field is linearly 
polarized, that the aperture dimensions and the radii 
of curvature of the reflector are large compared with a 
wavelength, and that the reflector is in the far-zone of 
the primary source, the diffraction problem reduces to 
a scalar one and the normalized diffraction integral, 
which gives the Fresnel pattern in a horizontal plane, 
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can be approximated as! 


1 
Ia) = if F(r)e~ikD*r?/8Ry J(yr) dr, (1) 
0 
In (1), D is the aperture diameter: R, the distance of 
measurement; 7, the radial dimension in the aperture 
plane normalized with respect to D/2; u=(rD/X) sin 
6, @ being the azimuth angle; k =27/: and F(r), the cir- 
cularly symmetrical amplitude illumination function 
over the aperture. The exponential term is the Fresnel- 
zone contribution; terms above the second order are 
neglected. 

When the primary source is displaced from the focus 
of a paraboloidal reflector by a distance € along the 
reflector axis in the direction away from the reflector 
with a view to simulating far-zone patterns in the 
quasi-near zone, there will be a relative phase variation 
over the aperture. The phase variation referred to the 
center point is found with good approximation to be? 

2e fe 


iss = — 


mest: <; [al | 


The normalized diffraction integral now becomes 


1 
I(u) = f F(r)ei*(8-D*r?/8B)y J (ur) dr. (3) 
0 
This integral is difficult to evaluate. However, exam- 
ination of (3) reveals that a far-zone pattern will be 
obtained if € is chosen to make the exponent vanish or 
equal to a constant. Unfortunately this cannot be done 
exactly for all values of 7. It has been found? that an 
acceptable and convenient approximation for 6 is 


72 ] 

Af? 

fe | 
D 

Eq. (4) is exact for r=0 (center) and r=1 (edge of 

aperture). For other values of r, the | 5| given by (4) is 

slightly larger than that given by (2); the error de- 


creases when the (f/D) ratio of the reflector increases. 
Using (4), the required amount of defocus can be read- 


ily determined: 
ane)! 
Sanna af) 


6 = —2e}/1— (4) 


(5) 


‘If it is desirable to express R in multiples of D?/h, 


_ R=nD?/y, (5) reduces to 


B 1 2 LA\2 
*-f1(Z) +() | (6) 
nN nN D 4 
1S, Silver, “Microwave Antenna Theory and Design,” M.I.T. 
Rad. ee aL 12, McGraw-Hill Book Co., Inc., New York, 


N. Y., ch. 6; 1949. 

2S. T. Moseley, 
loidal Reflector,” Final Report for Contract No. 
Syracuse Univ. Res. Inst.; August, 1954. 
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Defocus Distance e€ in Multiples of A 


.20 .50 1.0 2 5 10 20 50 
Distance from Antenna in Multiples of D@/r 


Fig. 1—Defocus chart for paraboloidal reflectors. 


Fig. 1 shows a plot of (€/A) vs x for five values of 
f/D in log-log scale; they are parallel straight lines. It 
is noted that for n=2 (R=2D?/X, which is normally 
considered as an adequate far-zone distance), appreci- 
able defocus is still necessary, the required amount 
being larger for larger f/D ratio. 


THE MINIMUM BEAMWIDTH CRITERION 


The previous derivation for the required amount of 
defocus is independent of the amplitude illumination 
function F(r). In this section, the criterion that proper 
defocusing give a minimum 3-db beamwidth will be 
used; the problem of simulating far-zone patterns in the 
quasi-near zone then reduces to that of finding a rela- 
tion between € and 7 such that w; is a minimum for 
I(0) 


I (uw) 


In order to evaluate the diffraction integral, an ampli- 
tude function of the type 


no =) 


will be assumed. Substitution of (4) and (8) in (3) re- 
sults, with y=27/D, 


20 logio = 3. (7) 


(8) 


| T(w)| = f GL yperieyI(ws)dy | (9) 


where 


(10) 
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Eq. (9) can be integrated through repeated integration 


by parts to give 
Se (bere 
m=0 “ 


In order to determine the values of uw, the infinite 
summation in (11) is separated into real and imaginary 
parts and a table is constructed of 20 log 10| 1(0)/I(u)| 
vs u for various values of B. The values that sat- 
isfy (7) are then determined by interpolation. Fig. 2 isa 
plot of half-power beamwidth, 2u,=kD sin 61, vs €/A with 
nm asa parameter for a paraboloidal reflector with f/D 
=0(.35. It is distinctly seen from the curves that for 
every value of , as far as the Fresnel approximations in 
(9) hold, there exists a proper defocus distance which 
gives a minimum main-lobe 3-db beamwidth; this de- 
focus distance corresponds to the condition B=0. The 
B=0 points are shown by crosses in Fig. 1. Setting B in 
(10) equal to zero, a relation identical to (6) is obtained. 


I m3(t) 


ae 
) QjBx)* Piro 5 (11) 


| 1(w) | =" 


Derrocus CHARACTERISTICS IN THE FRESNEL REGION 


The curves in Fig. 2 are useful not only because they 
give the proper amount of defocus in order to simulate 
far-zone radiation patterns at distances of R=2D?/X, 
1D?2/\, 4D?2/r, and 4D?/)\ for an antenna with (1—-y”)? 
illumination having f/D=0.35, but also because they 
indicate the extent to which the 3-db beamwidth is in- 
creased if the primary feed is moved by a given amount 
on either side of the proper defocus position. The fact 
that for different values of m the minimum 3-db beam- 
width remains unchanged when the primary feed is de- 
focused properly is plausible inasmuch as proper de- 
focus will simulate the same far-zone pattern. It is also 
clear that the initial minimum beamwidth criterion 
could have been made with reference to other decibel 
levels, instead of —3 db, without change in the out- 
come of the required amount of defocus. 

Referring again to (9) one observes that for wu 2.405, 


[ru B)| sf (1s) a(uydy = T(u,0). (12) 


Eq. (12) states that the condition B=0 gives the least 
upper bound for the magnitude of J(u). From this, it 
follows directly that the magnitude of the defocus in- 
tegral is maximum with B=0 (proper defocus) and, 
in particular, | 7(0) | is maximum. There is no danger 
of main-lobe deterioration when the criterion of mini- 
mum beamwidth is met. 

The curves in Fig. 2 are symmetrical with respect to 
the minimum beamwidth point (B =0). This can be seen 
from an examination of (11), the magnitude of J(u) 
being independent of the sign of B. Numerical evalua- 
tion of theoretical radiation patterns in the Fresnel 
region for various defocus distances (performed on the 
Mark IV Computer at Harvard Computation Labora- 
tory) checked well with the predicted results. 
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Fig. 2—Beamwidth characteristics for a defocused paraboloidal 
reflector, f/D =0.35. (1 —y*)? illumination. 


CONCLUSION 


The defocusing practice of simulating far-zone radia- 
tion patterns in the quasi-near zone has been evaluated 
analytically and validated for the paraboloidal reflector. 
Using the minimum main-lobe beamwidth criterion, the 
defocus characteristics for an amplitude illumination 
function of the (1—y?)? type have been obtained, which 
enable one to determine the proper amount of defocus 
as well as to predict the increase in beamwidth when the 
antenna system is not properly defocused. It has also 
been shown that the minimum beamwidth criterion 
coincides with the condition for maximum radiation 
amplitude along the reflector axis. 

For a given paraboloidal reflector, the proper amount 
of defocus in multiples of wavelength (€/A) is approxi- 
mately inversely proportional to the distance at which 
measurements are to be made; and for measurements at 
a given distance, the required amount of defocus in- 
creases with the f/D ratio of the reflector in a parabolic 
manner, as shown by (6) and Fig. 1. 

The fundamental requirement for simulating far-zone 
radiation patterns in the quasi-near zone is equivalent to 
that for the cancellation of the effect of the additional 
phase term in the diffraction integral for the Fresnel 
region by the phase term introduced by defocusing. 
Neither a change in the amplitude illumination function 
nor a change in the integration limits would alter this 
condition. Actual pattern characteristics, such as beam- 
width, radiation intensity, etc., would of course depend 
upon both the amplitude function and the integration 
limits. 
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Available Bandwidth in 200-Mile VHF 
Tropospheric Propagation * 
L. A. AMES} anv T. F. ROGERS} 


spheric paths extending far beyond the radio 

horizon indicated that the propagation medium 
might support a relatively distortion-free bandwidth of 
several megacycles.! More recently, the direct trans- 
mission of television signals to distances beyond the 
radio horizon on both vhf? and shf* has been reported, 
with no evidence of degradation in signal quality due 
to the propagation medium. Reception of the shf signals 
over a 180-mile path took place well within the extra- 
diffraction field strength region, but measurements were 
confined to July and August; the year-round vhf meas- 
urements were made at 120 miles, and because of the 
relatively long wave-length and short path length, were 
influenced by diffraction propagation. In order to de- 
termine whether wide-band signal transmission is gen- 
erally possible over distances far in excess of the radio 
horizon, measurements must be made both under cold 
weather propagation conditions and at such a distance 
that the contribution of the diffraction field—already 
‘established as capable of supporting wide-band signal 
transmission—is negligibly small. 

During the coldest winter months extra-diffraction 
tropospheric signal levels are appreciably weaker than 
the yearly median level, and much weaker than the 
summer maximum.’* This cold weather reduction in 
median signal level may be attributed in part, to the 
near absence of tropospheric refractive index gradients 
or layers strong enough to cause super-refraction or re- 
flection propagation. Therefore, to the extent that 
scatter propagation® contributes to extra-diffraction 
field signal levels, this mode of propagation should be 
observable in its purest form in winter, and available 
bandwidth should be measured under such conditions. 

In our experiments we have utilized the signals trans- 
mitted by WATV which are now radiated from the Em- 
pire State Building in New York City. These television 


HE sober SHORT-PULSE experiments on tropo- 


* Original manuscript received, June 1, 1955. 

+ AF Cambridge Research Center, Cambridge, Mass. | ; 

1K. Bullington, “Radio propagation beyond the horizon in the 
40-4000 mc band,” Proc. IRE, vol. 41, pp. 132-135; January, 1953. 

2 J. C. Simon and V. Biggi, “Un nouveau type d’aérien et son 
application 4 la transmission de télévision a grande distance, 
L’Onde Elect., vol. XXXIV, pp. 883-896; November, 1954. 

3K. Bullington and W. H. Tidd, Panel Discussion: “Extended- 
range vhf and uhf propagation,” PRoc. IRE, vol. 43, p. 365; March, 
1955. : Ay 
4 For the New England Area see, for instance, G. S. Wickizer and 
A. M. Braaten, “Field strengths recorded on adjacent FM channels 
at 93 mc over distances from 40 to 150 miles,” Proc. IRE, vol. 40, 

p. 1694-1699; December, 1952. Our own measurements (unpub- 

fished) on a 220 mcps 200-mile overwater path agree in this particular 
with these lower frequency measurements. f , 

5 H. G. Booker and W. E. Gordon, “A theory of radio scattering 
in the troposphere,” Proc. IRE, vol. 38, pp. 401-412; April, 1950. 


transmissions cover the frequency range 210-216 mcps, 
the effective radiated peak power is 316 kw, and the 
horizontal radiation pattern is nearly omnidirectional. 
The signals are received 196 statute miles away in 
Scituate, Massachusetts, where the smooth earth dif- 
fraction field is calculated to be some 135 decibels below 
free space. A large 40’X70’ dipole receiving array 30’ 
above the immediate terrain is used in conjunction with 
a sensitive high-quality television receiver. The trans- 
mitting and receiving parameters indicate a minimum 
detectable signal level some 82 db below free space, 
thereby providing sufficient sensitivity to permit cold 
weather observations to be made. 

Photographs of the WATV test pattern were made 
during the entire first quarter of 1955, intermittently at 
first, and systematically after mid-February. An ex- 
posure was made every five minutes during the hour in 
which the test pattern was displayed—usually in the 
early morning—each day of the week. Over 600 photo- 
graphs were obtained in this manner, and distinguish- 
able test patterns can be seen in about 90 per cent of 
them. The median signal level for this period was close 
to 7 decibels above the minimum detectable signal, ‘in 
excess of 20 db about 5-10 per cent of the time and 
maximum values exceeded 30 db. 

Close examination of these test patterns indicates 
that, with one exception, there is no discernible distor- 
tion or loss of picture detail except that caused by poor 
signal-to-noise ratio; certainly no marked distortions 
have ever been seen. Fig. 1 presents two photographs, 
representative of signals approximately 36 db and 
15 db, respectively, above the minimum detectable and 
it may be noted that, except for the deterioration 
caused by receiver noise in the weaker signal case, the 
picture quality is maintained; this is true for all the 
photographs taken. The conclusion may be drawn, 
therefore, that the propagation medium allows an effec- 
tive bandwidth of some 3-4 mcps at vhf for 200 miles 
even during the weaker signal level periods predomi- 
nantly characterized by the extra-diffraction type prop- 
agation mode. While an “effective” bandwidth, meas- 
ured in this manner, cannot be compared directly with 
the bandwidth predictions of Rice’ or Gordon’ such a 
conclusion is in essential agreement with their calcula- 
tions which are based upon the tropospheric scatter 
propagation theory. 


6S. O. Rice, “Statistical fluctuations of radio field strength far 
beyond the horizon,” Proc. IRE, vol. 41, pp. 274-281; February, 
1953. 

7W. E. Gordon, “Radio scattering in the troposphere,” Proc. 
IRE, vol. 43, pp. 23-28; January, 1955. 
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Fig. 1—(a) s/Smin %36 db, (b) 5/smin ¥15 db. 


The one exception observed to the above conclusion 
is that of multipath propagation caused by signal re- 
flection from aircraft flying across the propagation path. 
At times, a characteristic’ rapid and intense variation 
occurs in the received signal, resulting in a fleeting 
“shost” time-displaced image. The use of narrower an- 


8 C. R. Englund, A. B. Crawford, and W. W. Mumford, “Ultra- 
short-wave transmission and atmospheric irregularities,” Bell Sys. 
Tech. Jour., vol. 27, pp. 489-519; October, 1938. (See esp. Fig. 19.) 


tenna beams should reduce the incidence of such oc- 
currences. 

Interesting discussions have been held with J. C. 
Simon of Compagnie Générale de Télégraphie Sans Fil, 
K. Bullington of the Bell Telephone Laboratories, and 
J. H. Chisholm of the Lincoln Laboratory, MIT. Thanks 
are due J. Frazier, E. Hodges, and D. Goode for assist- 
ance in making the measurements, and to E. Martin in 
reducing the data. 


Bibliography of Nonuniform Transmission Lines* 
H. KAUFMAN} 


lines has grown slowly but steadily since original 
papers by Heaviside appeared at close of the nine- 
teenth century. (Heaviside’s contributions have been 
almost completely overlooked by later writers.) The 
number of papers, although not overwhelming, is now 
large enough to warrant preparation of a bibliography. 
References to periodicals are in Part I, to texts in Part 
II, and to a selected group on horn theory in Part III. 
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Part I—Report on Microwave Optics 
R. C. SPENCER} 


HE ANTENNA or propagation engineer has long 
‘es familiar with antenna theory and electromag- 

netic theory but of late he has been hearing the 
term microwave optics. This is a term introduced during 
World War II to cover those aspects of microwaves that 
could readily be studied by means of optical analogies. 
That there is no sharp demarcation between these fields 
can be gathered from the fact that a symposium on any 
one of them inevitably includes papers relating to the 
other two. 

In its narrow sense, microwave optics makes use of 
analogies in geometric and physical optics. Thus de- 
fined, microwave optics is an interesting point of view, 
useful pedagogically for teaching either students or lay- 
men, and useful scientifically for obtaining quick if 
limited answers to problems. In its broader sense, micro- 
wave optics draws on the more rigorous techniques of 
diffraction and electromagnetic theory. Operational or 
Fourier techniques are used as tools to exploit the meth- 
ods of physical optics and electromagnetic theory, and 
also to analyze the role of the microwave optical system 
as a link in a communication chain. At this point certain 
formal aspects of communication theory or information 
theory can be applied. 

Commission VIc has been concerned with the broader 
definition—in particular, with the inputs from Fourier 
transforms and information theory. 


* Original manuscripts received, April 25, 1955. A report by 
URSI Subcommission VIc at the XIth General Assembly, The 
Hague; August 23-September 2, 1954. 

} Chairman, URSI Subcommission VIc, Microwave Optics. Chief 
Antenna Lab., AF Cambridge Res. Center, Cambridge, Mass. Now 
with Sylvania Electric Products, Inc., Waltham, Mass. 


HIstTory OF SESSIONS ON MICROWAVE OPTICS 


During World War II, microwave equipment was 
mass-produced for various radar and communication 
applications. In most of this equipment the antenna was 
designed to control the shape of the transmitted or re- 
ceived beam as well as to effect an impedance match. 
The use of microwaves stimulated studies in antenna 
theory and design, artificial dielectrics, diffraction, and 
electromagnetic theory, and interest in the new field led 
to a number of symposia on microwave optics. The fol- 
lowing list of meetings held is exclusive of the various 
radio engineering conventions at which papers on the 
subject were given in the sessions on antennas. 

Oct., 1947. Optical Society of America, Cincinnati, 
Ohio. Symposium on Microwave Optics 
(one-half day) 

Apr., 1951. Joint URSI-IRE meeting, Washington, 
D. C. Symposium on Microwave Optics 
(one-half day) 

Oct., 1951. National Bureau of Standards, Washing- 
ton, D. C. Symposium on Quality of the 
Optical Image (three days). Although this 
symposium did not concern itself with 
microwave optics, it proved to be stimu- 


lating to those in attendance who were in- 


terested in this analogous field. (The Pro- 
ceedings have recently been issued as Na- 
tional Bureau of Standards Circular No. 
526. See papers by Zernike, Marechal, 
Osterberg, McDonald, and Scade; see also 
discussions by Spencer, pp. 38, 249; and 
by Grey, p. 251). 


1955 


June, 1952. International Physics Union, Milan, Italy. 


(three days). Symposium on Microwave 
Optics. This was the first international 
symposium devoted entirely to the subject 
of microwave optics. (Proceedings have 
been published in Supplemento No. oF 
Nuovo Cimento, vol. 9, 1952). 
Tenth General Assembly of URS] ‘Syd- 
ney, Australia. Several papers on micro- 
wave optics were presented to Commis- 
sion VI. Subcommission Vic, Microwave 
Optics, was formed particularly to study 
inputs from other fields such as Fourier 
techniques and information theory. 
Symposium on Microwave Optics, McGill 
University, Montreal, P.Q., Canada. This 
symposium was planned to solve some of 
the problems of Commission VIc. Note re- 
views of the Symposium by E. Wolf in 
Nature, October 3, 1953, pp. 615-616, and 
by Morris Kline in Physics Today, March, 
1954, pp. 30-32. The Proceedings will be 
published in two volumes by the Air 
Force Cambridge Research Center, Air 
Research Development Command (224 
Albany Street, Cambridge 39, Massa- 
chusetts, U.S.A.), in the fall of 1955. 

The following symposia were also scheduled during 
the summer of 1954: 


Aug., 1952. 


June, 1953. 


Aug., 1954. Report of URSI Commission VIc, Micro- 
wave Optics, to 11th General Assembly, 
The Hague, The Netherlands. 

Sept., 1954. Florence, Italy, Meeting on Problems in 


Contemporary Optics, Session on Non- 
classical Image- Forming Devices. See “Pre- 
liminary Remarks on Microwave Optics,” 
by R. C. Spencer, to appear in the PRo- 
CEEDINGS. 


INTERACTION OF VARIOUS FIELDS WITH 
MICROWAVE OPTIcs 

You may be interested in a topologic diagram (Fig. 1) 
that illustrates the dependence of microwave optics on 
other fields. In the case of microwave optics shown at 
the bottom of the diagram, we are indeed grateful to the 
inputs from Physics on the left which include, in order of 
ascending difficulty, geometric optics, physical optics, 
and electromagnetic theory. On the other hand, micro- 
waves have both drawn from and contributed to the 
Electronics fields of circuit theory and communication, 
shown on the right. At the top of the diagram are the 
operational methods. These usually lead to Fourier 
transforms for problems in physics and microwave 
optics, and to the Laplace transform for problems in 
electronics. More recently, information theory has been 
applied to optics and microwave optics. 


Report or URSI Suscommission VIc 


In addition to the preceding historical remarks and a 
short discussion of Fourier methods, the report to Com- 
mission VI included the following survey papers: 
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MICROWAVE 
OPTICS 


Fig. 1—Relationships between microwave optics and other fields. 


1. Professor E. Wolf, “Relationship Between Geo- 
metric Optics and Diffraction Theory.” This included 
certain energy flow relationships which hold both for 
geometric optics and electromagnetic theory. 

2. Dr. H. Bremmer, “Diffraction Problems of Micro- 
wave Optics.” This summarized recent contributions, in- 
cluding monograph by Bouwkamp and papers presented 
at McGill Symposium on Microwave Optics. 

3. Professor A. Blanc-Lapierre, “Information Theory 
and Optics.” This paper will be published in French in 
Optica. It is significant that although each of these three 
papers treated optical problems, they were presented by 
request to the International Scientific Radio Union be- 
cause of the direct application of the ideas to microwave 
optics problems. 


CURRENT PROGRAM AND MEMBERSHIP 


The original membership of URSI Subcommission 
Vic as formed in August, 1952, included Professor G. A. 
Woonton of McGill University, P.Q., Canada, and 
P. M. Woodward of TRE, Malvern, England, but pres- 
sure of other duties forced their withdrawal. The present 
membership consists of A. Blanc-Lapierre, H. Bremmer, 
E. Wolf, and R. C. Spencer. 

At the 11th General Assembly, Subcommission VIc 
was reconstituted as a “Working Group” on Microwave 
Optics under Commission VI and has been asked to pre- 
pare a monograph on Microwave Optics and Its Relation 
to Other Fields. This will be presented at the 12th Gen- 
eral Assembly of URSI in Denver, Colorado, USA, in 
August, 1957. 

In ending these preliminary remarks I would like to 
point out that microwaves, even more than the longer- 
wavelength radio waves, have stimulated the applica- 
tion of optical techniques to laboratory or field condi- 
tions. This is because the average laboratory experiment 
can both measure the phase and amplitude of points in 
microwave fields and yet deal with a sufficiently large 
number of waves so as to approximate optical condi- 
tions. 

Until now the flow of ideas has been mainly from the 
other fields represented in the chart to microwave 
optics. I believe that now some of the improved tech- 
niques developed for microwave optics can contribute 
to these other fields. 
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Part Il—Diffraction Problems of Microwave Optics 
H. BREMMER?} 


INTRODUCTION 


rently applied in diffraction theories of micro- 

waves. It is not a total summing up of the com- 
plete literature, references being restricted mainly to 
papers given in the Proceedings of the Symposium on 
Microwave Optics held at McGill University, Montreal, 
June 22-25, 1953.1 These Proceedings, to be published by 
the Air Force Cambridge Research Center in the fall of 
1955, give an excellent survey of recent theoretical and 
experimental results. Bouwkamp’s paper “Diffraction 
Theory,” in volume 17 of Reports on Progress in Physics, 
is also highly recommended for a study of diffraction 
theories in general. 

Geometric optics concerns mainly ray trajectories 
that satisfy the reflection law and Snell’s refraction law. 
The transportation of energy along the rays implies that 
the field amplitudes vary in inverse proportion to the 
square root of the cross section of a pencil of neighboring 
rays, that is, in direct proportion to the square root of 
the local Gaussian curvature K, of the wavefront 
passing through the point of observation. The geometric 
optics approximations are sufficiently accurate for the 
design of microwave lenses even when aberration effects 
are taken into account (see SM 34, 36, 37, 39). How- 
ever, geometric optics breaks down near caustic sur- 
faces and in shadow regions, the field there becoming 
infinite and zero respectively. In these cases a wave 
treatment becomes inevitable. Diffraction describes the 
wave phenomena in this domain. 


lke STATUS report deals with methods cur- 


A. SCALAR THEORIES 


In free space the single scalar wave equation 
(A+k?)u=0 holds independently for every rectangular 
field component. The complications in the problems 
arise from the coupling of these components by the 
boundary conditions. The effect of this coupling is there- 
fore significant first of all near the boundaries; where 
microwaves are concerned, however, the field at the 
boundaries can be checked experimentally since these 
waves admit field measurements up to a distance of one 
wavelength from the boundaries. The influence of the 
boundary conditions explains why the analogy between 
acoustic and electromagnetic diffraction through aper- 
tures and around disks (see Bekefi SM 53) holds in the 
far field but fails in the near field. In view of these 
circumstances, and discounting the fact that the field 
across the aperture was incorrectly assumed as equal to 


} Philips Res. Labs., N. V. Philips’ Gloeilampenfabrieken, Eind- 
hoven, Netherlands. 

1 Such references are identified by the initials SM after the 
author’s name, followed by the publication number of the paper. 


that of the incident wave, it is remarkable that a scalar 
treatment of the dominant electric field component 
could nevertheless account for the empirical field values 
observed by Andrews (.SM 49) along the axis of a circu- 
lar aperture. 

The most satisfactory diffraction theories start from 
Maxwell’s equations. Even then, a scalar treatment is 
not impossible in principle since the fields satisfying 
these equations are expressible in terms of some complex 
scalar quantity. Examples are: (1) a complex combina- 
tion of the two scalars describing the field of moving 
charges in a theory of Whittaker’s,? (2) the scalar intro- 
duced in free space by Wolf (SM 17) in order to arrive 
at expressions for the energy and momentum densities 
similar to those in quantum mechanics, and (3), a com- 
plex combination of two scalars describing the field in 
terms of multipole contributions in a recent paper by 
Bouwkamp and Casimir. A drawback in all these 
scalars being the inconvenient form of the corresponding 
boundary conditions, an explicit vectorial treatment 
becomes more convenient. In problems dealing with 
phases (as defined in geometric optics), which are the 
same for all the field components, the inconvenience is 
slight and aberration effects of microwave lenses can 
therefore be investigated very simply by scalars. Dun- 
bar (SM 40) thus found that the maximum tolerable 
spherical aberration producing a confusion equal to that 
caused by diffraction (Airy disk) is one-fourth as large 
as the one imposed by the Rayleigh criterion; in prob- 
lems of astigmatism and defocusing the two correspond- 
ing tolerances proved to be of equal magnitudes. 

The useful vector diffraction theories can be sepa- 
rated according to those dealing with corrections for the 
geometric optics approximation (see Section B), and 
those starting from very general properties of Maxwell’s 
equations without any explicit connection with geo- 
metric optics (see Section C). In the latter case, geo- 
metric optics is recognized as the limit for A\—-0 by 
applying saddlepoint approximations. 


B. VECTOR DIFFRACTION THEORY AS A CORRECTION 
FOR GEOMETRIC OPTICS 


The Maxwell equations for a monochromatic field 
(time factor e~‘* in a source-free inhomogeneous me- 
dium (dielectric constant €, permeability u, refractive 
index n=(eu)'? are comprised in the single vectorial 
equation 


curl M + knM = 1/4 grad log — XM* (1) 
€ 


? Whittaker, Proc. London Math. Soc., vol. 2, pp. 367, 1904. 
* Bouwkamp and Casimir, Physics, vol. 20, pp. 539; 1954. 
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| 
‘or the complex vector (with imaginary unit /) 
M = (u/e)4H + i(e/u)'/4E, 


which was introduced by Bateman.‘ Luneburg and 
Kline’s theory® amounts to the following in general 
asymptotic expansion 


mn (Vey, 
r=0 ke 


in which the function S of the space variables should be 
independent of k. The two independent functions for H 
and E result from the real and imaginary parts with 
respect to j while treating H, E, and the other imaginary 
unit 7 as real quantities. A substitution of this expansion 
into (1) leads to the following relation resulting from the 
terms of order zero in k: 


grad SXmp = jnmo, 


whose vector components prove to be compatible only 
if S satisfies the conventional eikonal equation 
| grad S| =n. A further relation, however, resulting 
from the compatibility of the first-order terms, shows 
that the mp field should vary proportionally with n./K, 
along a ray trajectory (defined as a curve perpendicular 
to the wavefronts S=constant), K, being the local 
gaussian curvature of the wavefront through the point 
of observation. In addition, the variation in the direc- 
tions of the field vectors contained in m is fixed by the 
property that these vectors suffer a parallel displace- 
ment in the Levy-Civita sense, the trajectories being 
geodesics in the non-Euclidean space defined by the line 
element 


ds? = n%(dx? + dy? + dz’). 


The mp field thus obtained as the limit for Rk or 
A—0 represents the geometric optics approximation. 
The relations resulting from the higher-order contribu- 
tions in (1) determine the variation of the vectors 
m;, m2, - ~~ along the trajectories, and thus fix the cor- 
rections for the geometric optics field. 

The splitting of the diffraction field as established 
here can be connected with pulse solutions, that is, with 
solutions describing the field that exists if all the sources 
simultaneously undergo sudden identical relative 
changes in intensity instead of varying oscillatorily. The 
resulting field discontinuities then pass in succession 
along the wavefronts otherwise associated with the sta- 
tionary solution. The vectors m, then prove to be pro- 
portional to the discontinuities of 6"m/t* for the pulse 
solutions at the moment the discontinuity passes 
through the point of observation. In particular, the 
seometric optics field is linked with the discontinuity in 
the field described by the pulse solutions. 


4 Bateman. 
5 Luneburg and Kline, Comm. Pure Appl. Math., vol. 4, pp. 225; 


1951 
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C. THEORIES BASED ON GENERAL SOLUTIONS 
OF THE MAXWELL EQUATIONS 


General solutions of the Maxwell equations may be 
divided into three groups according to their unique 
suitability in the special problems involved in: (a) mi- 
crowave lenses; (6) diffracting screens and disks, reflec- 
tors, and scatterers; and (c) waveguides. 


a. Solutions Usable for Microwave Lenses 


For the case of negligible reflection losses, Luneburg® 
arrived at the following vectorial solution of Maxwell's 
equations in the image space (e=u=1) beyond the lens: 


Lk Latin 415 
MP = J f Ocv {Rs(O) mao QetH wor ceruer, (2) 


Ray trajectorics are produced beyond the image 
space. The integration extends over a wavefront that 
cuts these at right angles, far away from the lens; K, is 
the Gaussian curvature of the wavefront (with surface 
element d0,); W(Q) is the mixed eikonal; r, the radius 
vector from the origin to the point of observation; and 
U,, the unit vector along the normal of the wavefront. 

Eq. (2) includes all information about optical instru- 
ments. If the direction of polarization is constant 
throughout space, the field derived from a scalar F, repre- 
sented by an integral similar to' that of (2), can serve as 
a basis for the investigation of optical aberrations (see 
the derivation of the spherical aberration of Nijboer and 
Chako, SM 21). The case of a nonpolarized source can 
be treated by averaging over all polarization directions; 
it will be found that the propagation of energy is along 
curves parallel to the gradient of the phase of the com- 
plex function F. Unlike the ray trajectories, these 
curves may penetrate into shadow regions which would 
explain the establishment of small field intensities there. 


b. Solutions Usable for Diffracting Screens and Disks, 
Reflectors, and Scatterers 


The field inside the space V enclosed between the sur- 
faces of these elements and infinity can be expressed in 
terms of the distribution of the electric and magnetic 
fields e and A across the finite parts of the boundaries. 
The formulas for the field inside V, in combined form, 
read: 


1 
6, Fi P= Crp, Nod) ate re curl eu dO,(e, hXu),e°*?2/P2 
T 2 


1 
+ — curl curl ff dou, eXu),e1#PelP2, (3) 


4k 


the first term of which represents the primary field that 
would be produced by the sources inside V if the enclos- 
ing walls constituting 2 were absent, while u is the unit 
vector along the outward normal of 2. 


6 Luneburg, “The Mathematical Theory of Optics,” Brown Uni- 
versity Press, Section 47; 1944. 
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The expression in (3) can be interpreted as the linear 
superposition of contributions produced by infinitesimal 
electric and magnetic dipoles on 2, the electric dipoles 
depending on the tangential components of the mag- 
netic field on =, and the magnetic dipoles depending on 
those of the electric field on ©. The actual field results 
from the radiation from these dipoles in accordance with 
Huygens’ principle (see Possel, SM 50). The field of the 
elementary dipoles can be represented conveniently by 
Green tensors, whose elements connect the vector for a 
dipole of unit moment, located somewhere on 2, with its 
field observed at a point P, taking into account the in- 
fluence of ©. Such a Green tensor appears to have been 
evaluated for the cases of free space and of a half-space 
limited by a plane (see Severin, SM 54). 

The difficulties encountered in diffraction problems 
are due to the interdependence of the electric and mag- 
netic fields tangential on Y. Eq. (3) constitutes the basis 
for the deduction of integral equations fixing one of 
these quantities. In the usually considered case of per- 
fectly conducting boundaries, such an equation can be 
expressed in terms of the vectorial density j of the cur- 
rents induced on Y. One of these equations, derived by 
Maue, reads: 


k Bue ethhe 1 
—ulP) X f f IO G+ 


etkPQ 


X grad ii dO, div J(Q) PO. 
= iu(P) X epr(P), 


and can even be applied to the boundaries of the in- 
finitely thin screens considered in the simplest diffrac- 
tion problems. Unfortunately, a straightforward attack 
through integral equations of this type or its equivalent 
has been successful in exceptional cases only, for in- 
stance, in Bouwkamp’s treatment of the diffraction 
through a small circular aperture, in which the addi- 
tional edge conditions along the rim play an important 
role.” But these equations may be the starting point for 
successive approximations (see Sections G and H). 


c. Solutions Usable for Waveguide Problems 


An example of a solution of this type is the so-called 
“configuration lens,” which consists of two cylindrical 
curved parallel metallic plates between which TE waves 
can be propagated. The principle can be described with 
the aid of curvilinear orthogonal co-ordinates. The cor- 
responding solution of Maxwell’s equations works out 
to a two-dimensional wave propagation along the geo- 
desics of the mean surfaces parallel to the metallic faces. 
The design of microwave configuration lenses thus de- 
pends on differential geometry. The theoretical and ex- 
perimental aspects have been reviewed by Myers. 


7 Philips Res. Reps., vol. 5, pp. 401; December, 1950. 
8 Myers, Jour. Appl. Phys., vol. 8, pp. 220; 1947. 
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D. APPROXIMATIVE METHODS DEPENDING 
ON KIRCHHOFF INTEGRALS 


When the boundary values of one single parameter 
over the diffracting or scattering surfaces are unknown, 
it is necessary to resort to approximative scalar meth- 
ods, these being particularly suitable for discussing ex- 
perimental results. The material boundaries of the prob- 
lem in question are then replaced by a surface S that 
consists partly of these boundaries and partly of im- 
aginary sheets filling up the open space between them. 

Unfortunately, the boundary conditions used for the 
material surfaces of S must obviously be different from 
those used for the open parts of S. In the case of aper- 
tures in a plane, for instance, conditions of symmetry, 
originally applied by Bethe, make it apparent that the 
normal derivatives of the tangential electric field and of 
the magnetic field normal to it must be identical with 
those of the primary field on the open parts of S, whereas 
the field components themselves should of course vanish 
on the perfectly conducting metal parts of S. For the 
rest, these rigorous boundary conditions lead to the vec- 
torial form of Babinet’s principle, which connects the 
two complementary diffraction problems in which the 
apertures and metal parts have been interchanged.? To 
obtain a simple approximative boundary condition for 
one and the same quantity throughout S (assumed to be 
a plane), the exact boundary condition for a quantity 
nonvanishing in the diffracting apertures or disks might 
be combined with a reasonable assumption for the same 
quantity in the remaining part of S (see Severin, SM 
54). The quantity in question could then be evaluated 
for any field from a Kirchhoff integral over the aperture 
plane, provided it satisfies the ordinary wave equation. 

Among others who have discussed the possible im- 
perfections of this method is Neugebauer (SM 51), who 
has considered such limiting cases as vanishing or in- 
finite apertures. The results derived by means of these 
methods nevertheless prove to agree with many experi- 
mental data, as has been emphasized by Andrews (SM 
49) with respect to measurements up to the aperture 
plane. By a related approximation method, Severin 
(SM 66) has determined diffraction field from a Kirch- 
hoff integral over conic limit of shadow cast by a point 
source behind an aperture, using geometric optics ap- 
proximation for the field on the shadow limit. 

Another approximative method interprets the second- 
ary field as originating from wavelets radiated from the 
points of the aperture rim. This concept is not incon- 
sistent with representation by a Kirchhoff integral over 
the aperture since, as shown by Rubinowicz, such an 
integral can be split into a geometric optic part and a 
contour integral along the rim (see also SM 64). An 
empirical formula based on this idea has been suc- 


ee “Diffraction Theory,” Reps. Prog. Phys., vol. 17, 
p. 46. 
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cessfully applied by Andrews (SM 49). According to the 
saddlepoint method, contributions from the rim arrive 
mainly from special points on latter (see Section E). 

All the boundary conditions considered so far have 
reference to a single field parameter. Further complica- 
tions arise when the finite conductivity of diffracting or 
reflecting screens has to be taken into account. Approxi- 
mations valid for high conductivities lead to mixed 
boundary conditions, depending on both the parameter 
in question and its normal derivative. Little work has 
been done so far in this field (see Felsen, SM 6). 


E. SADDLEPOINT APPROXIMATIONS 


Exponential factor of diffraction integrals enables ap- 
plication of saddlepoint approximations. These are use- 
ful in working with wavelengths that are short compared 
with dimensions of diffracting or scattering body, and 
in evaluating radiation field (see Section F). 

This particular exponent represents phase of con- 
tribution from a special point on surface of integration, 
whereas extreme values at saddlepoints correspond to 
optical path lengths of ray trajectories that connect 
source of field with point of observation (Fermat’s 
principle). Geometric optics constitutes the limit of 
Wave optics inasmuch as the smaller the wavelength 
the more the contributions of the surface integral are 
restricted to the immediate vicinity of the saddlepoints. 

The problem presented by diffraction integrals over 
smoothly curved surfaces is not the same as that pre- 
sented by diffraction integrals over surfaces bounded by 
sharp edges such as apertures. On the one hand, with 
smoothly curved surfaces, the contributions due to the 
two-dimensional saddlepoints can be approximated by 
cutting off the Taylor expansion of the exponent, taken 
around the saddlepoint, after the second-order terms 
while extending the integration at both sides up to in- 
finite values of the two co-ordinates. The second-order 
approximation thus obtained proves to be identical with 
the geometric optics approximation. 

On the other hand, diffraction integrals that extend to 
an edge involve one-dimensional saddlepoints connected 
with positions on the rim. A corresponding ray tra- 
jectory with a turning point at such a position is char- 
acterized by an extreme path length. The trajectories 
through these so-called “critical points of the second 
kind”!° can be included in geometric optics considera- 
tions by extending Fermat’s principle. This extension 
worked out on a formal basis by Keller (SM 20). 

The second-order saddlepoint approximations of ordi- 
nary diffraction integrals over unbounded surfaces 
break down if the Gaussian curvature of the wavefront 
of the corresponding rays becomes infinite, that is, on 
caustic surfaces. The simplest approximation then takes 
into account the smallest number of high-order terms of 


10 Van Kampen, Physica, vol. 14, pp. 575; 1948-1949, 
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the mentioned Taylor expansion, which leads to a finite 
expression even on the caustic. In the vicinity of ordi- 
nary points on caustics, this procedure results in an 
Airy integral, whose value involves a maximum of the 
field near the caustic, and an exponential decrease at 
one side (shadow region impenetrable to rays). The field 
on the caustic increases as k!/6 as \—0, with a phase 
shift of 7/4 relative to the geometric-optics value. At 
foci, that is at cusps of the caustic surface, the field in- 
creases as k!/4, the phase shift being 7/2. Other singulari- 
ties may lead to different phase shifts and powers of k. 
Further, since the diffraction fields associated with the 
critical points of the second kind are proportional to 
k/?, they become negligible with respect to the ordinary 
geometric optics fields as k> ~. 

The field distributions near caustic surfaces have, in 
particular, been investigated for rotationally symmetric 
systems (in connection with spherical aberration) and 
for systems with cylindrical symmetry (see Kay, SM 42, 
for rays reflected by a cylindrical wall). Derivation of 
geometric properties, and of position of caustic surfaces 
is a matter of differential geometry (Parke, SM 38). 


F. RADIATION PATTERN 


The limiting field distribution existing at distances 
that are large by comparison with both the wavelength 
and dimensions of diffracting and scattering bodies can 
be obtained, for example, by substituting asymptotic 
expansions for the integrals or series representing the 
field throughout space. Distance beyond which such an 
expression may be used for checking experimental re- 
sults can be derived from the condition that remainder 
of expansion should represent contributions beyond 
accuracy of measurements (see Siegel, SJ 8). 

Although the radiation field can be expressed in terms 
of the current distributions on the conducting surfaces, 
for microwave applications it is more important to ex- 
press it in terms of its dependence on the field distribu- 
tion in an aperture. For a plane aperture, the far field 
proves to be proportional to the amplitude of a single 
Fourier component of the two-dimensional aperture dis- 
tribution. This Fourier component concerns periodic 
repetitions of the field along intersections of the aper- 
ture with planes normal to the direction of observation, 
occurring one wavelength apart. The radiation pattern 
of the aperture distribution is not affected by details 
that are finer than one wavelength; the part of the field 
that corresponds to the Fourier spectrum of these finer 
details does not radiate (evanescent waves). These re- 
sults can be derived from a saddlepoint approximation 
for the rigorous integral representing any quantity that 
satisfies the wave equation. 

These illustrations of the relationship between the dis- 
tant field and the near field show that from any known 
expansion of the near field it is possible to derive corre- 
sponding terms for the far field (Stevenson, SM 53). 
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The relationship between these terms is simplest if the 
terms of a rigorous expansion have well-defined limiting 
values at great distances. This is true, for example, of 
the expansion in Mathieu functions of the field for an 
infinite slit in a plane aperture (Leonard, SM 19). 
Knowledge of the far field is not only important for 
antenna design, but also for study of scattered fields 
that are far fields produced by reflection. Investigation 
of these fields by Crispin (SM 10) shows that radar cross 
sections can in special circumstances be increased by 
moving receiver away from the transmitter. The reci- 
procity theorem makes it possible to replace determina- 
tion of radiation patterns by that of currents induced by 
a plane incident wave on a receiving system; this may be 
useful in study of slot antennas (see Sinclair, SM 58). 


G. INTEGRAL EQUATIONS AS A BASIS 
FOR FURTHER COMPUTATIONS 


Integral equations such as Maue’s (see Section C) are 
suitable not only for numerical approach to diffraction 
problems in general but for a proper physical under- 
standing of scattering phenomena. Successive terms of 
the Liouville-Neumann expansion of such an equation, 
whose right-hand side depends on primary field, can be 
interpreted as contributions that are due to repeated 
scatterings. This is particularly clear if current distribu- 
tions to be derived from the equation belong to two or 
more spatially separated bodies, for instance, an an- 
tenna feed and a reflector opposite it. Iterative kernels 
then contain contributions that depend on integrations 
over different bodies; interaction through scattering is 
expressed in a mathematical form. Scattering produced 
by a set of identical cylindrical reflectors has been in- 
vestigated by Twersky and by Karp (SM 63). 

The wave theory of microwave reflectors is connected 
with the first two terms of the Liouville expansion. For 
reflectors with radii of curvature that are large as com- 
pared with the wavelength, the second term proves to be 
approximately equal to the first term; and the first term 
itself corresponds to the primary current distribution, 
defined here as the distribution induced directly by the 
primary field in the absence of any interaction with the 
other parts of the scatterer. The far field can thus be 
derived approximately by assuming a current distribu- 
tion equal to twice the primary distribution. The field 
derived by geometric optics from the reflection against 
the scatterer proves to be identical with this far field. 
In other words, the corrections in the geometric optics 
theory correspond to the amounts that the induced cur- 
rents deviate from twice the value of the primary dis- 
tribution. An example of an integral equation for an- 
other parameter is the equation starting from a scalar 
Green’s formula, derived by Franz (SM 61) for the field 
on an infinite circular diffracting cylinder. The evalua- 
tion of the integral term by the saddlepoint method 
enables the investigation of a creeping wave, that is, a 
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wave propagated along the surface of the cylinder in the 
part not illuminated by the incident wave (see also 
Imai’s paper! mentioned in Section Gd). 


G. Various Mretuops APPLIED IN MICROWAVE 
DIFFRACTION PROBLEMS 


In addition to the general principles discussed before, 
the following special techniques have been applied and 
may also be useful for future work. 


a. Expansion in Modes 


Mode solutions of the proper differential equations 
satisfy the necessary boundary conditions, including the 
radiation condition at infinity, but their amplitudes 
have to be derived from data concerning the primary 
field. For a point source this amounts to the construc- 
tion of a linear combination having the required singu- 
larity at the source. The method has been applied to the 
theory of wave propagation around the earth, and has 
been used by Jasik in the wave theory of the inhomo- 
geneous Luneburg lens (SM 32). 

For extended sources, such as the aperture of a slot 
antenna, the amplitudes can be derived through the aid 
of the Lorentz relation 


ff (Ey X H2)ndO = ff (E. X H1),dO 


for the integral over a closed surface of the indicated 
normal components which refer to two solutions of 
Maxwell’s equations. This identity replaces Green’s 
theorem for scalar wave problems. If we identify the 
two fields with the actual field and with that of a single 
mode, respectively, and integrate over the perfectly 
conducting screen and over an aperture in it, this rela- 
tion yields the amplitude of the mode in question in the 
expansion of the actual field. A possible complication is 
that the analytic form of the radial factor of the mode 
often has to be taken differently in the near field and in 
the far field (see Bailin and Silver, SM 1; and Tilston, 
SM 11). 


b. Watson Transformations 


These transformations concern reduction of mode 
and other expansions into an integral in the complex 
plane. This integral may be the basis for a suitable 
saddlepoint approximation or for derivation of a new 
mode series converging more rapidly than the original 
one. Originally introduced by Watson for the wave 
propagation around the earth, the method has been suc- 
cessfully applied by Row (SM 57), by Felsen (SM 56) to 
the diffraction problem of a wedge, and by Riblet (SM 
62) and by Imai," to that of a cylinder. The method can 
often be interpreted as a transformation from an expan- 
sion in azimuthal modes into one in radial modes. 


1 Imai, Z. Phys., vol. 137, pp. 31; 1954. 


1956 


c. Variational Principles 


Numerical solution of integral equations, such as for 
current distributions on conducting parts of the system, 
can be facilitated by deducting a quantity that, relative 
to actual distribution, remains stationary with respect 
to fictitious changes in it. According to Rayleigh-Ritz’s 
principle, approximations can then be derived by start- 
ing from a reasonably assumed analytic expression for 
the unknown distribution and choosing its parameters 
so as to minimize the quantity that has to remain sta- 
tionary (for an example, see Levine, SM 66). General 
stationary quantities derived by Levine and Schwinger 
are complicated but stationary expression deducted by 
Kouyoumjian (SM 14) for radar cross sections in terms 
of current distribution on a scatterer hit by a plane- 
polarized wave is rather simple. The variational method 
has enabled Huang (SM 13) to verify a satisfactory 
agreement between theory and experimental results for 
transmission coefficient of plane waves diffracted 
through apertures of various shapes. 


d. Wiener-Hopf Technique 


The two-dimensional integral equations of diffraction 
problems generally depend on, among others, the dis- 
tance of an integration point to a fixed point, that is, on 
the differences of three-dimensional co-ordinates. Owing 
to special conditions of symmetry the general equations 
may occasionally be reduced to one or more one-dimen- 
sional equations depending on a difference kernel only, 
at least for a limited range of the main variable. Explicit 
analytic solutions can then be derived by means of the 
well-known Wiener-Hopf technique. Three examples 
are: (1) the determination by Heins” of the current dis- 
tribution in an infinite set of parallel plates hit by a 
plane-incident wave; (2) the investigation by Kay and 
- Kotik (SM 16) of the field of a surface wave guided by a 
thin dielectric panel of limited width in one direction, 
but semi-infinite in the other; (3) the reduction by Blass 
(SM 22) of the scalar wave propagation through a semi- 
infinite cylinder with elliptic cross section to an infinite 
number of Wiener-Hopf equations. 


e. Applications of the Rules of the Operational Calculus 
to the Fourier Transforms Involved in the Connection be- 
tween the Radiation Pattern and Field Distributions 
across an Aperture 


(See Spencer, SM 43.) For instance, coefficients of 
expansions of the radiation pattern can be expressed in 
terms of higher-order moments of the aperture distribu- 

tion by applying such rules. Further, the influence of 
phase errors represented by a factor exp lip(x, y)| 
amounts to the multiplication of the undisturbed trans- 
form f(w:, #2) by a symbolic factor 


2 Heins, Quart. Appl. Math., vol. 8, pp. 281; 1950. 
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exp [i6(— 16/51, —i5/6we) J. 


f. Procedures for Improving the Convergence of Special 
Series 


As an example we refer to Schensted’s (SM 9) trans- 
formation of a known expansion for the backscattering 
by a half-infinite cone into a more rapidly converging 
series by means of the formula 

vay). 


co 1 0 
Oye eis a 
Soon | at 3 Cont 
g. Procedures for the Solution of Integral Equations De- 
pending on a Small Parameter (Usually the Product of the 
Wave Number and Some Characteristic Distance) 


Yen (SM 7) pointed out that the method worked out 
by Bouwkamp for the rigorous vectorial solution of the 
diffraction by a circular aperture can also be applied to 
the other problems. If the kernel of the equation 


f(s) f 166, s’)ds’ 


admits an expansion 
KG yes) = 2K ts, 5), 
then the coefficients of 
I(s) = Ze"I,(s) 


are to be derived from the equation: 


¥ f 160K, S )aseStae its), 


with 6n°=0 for n+0, and 60°=1. All these relations re- 
duce in succession to inhomogeneous integral equations, 
with different known right-hand members but with the 
same kernel Ko(s, s’). Expansions of this type have been 
derived by Stevenson" in connection with (3), 2 being 
a scattering boundary. 


H. THE PARAMETER OBSERVABLE IN 
FIELD-STRENGTH MEASUREMENTS 


The well-known difficulty in such measurements is 
that the quantity being measured is disturbed by the 
measuring equipment. An application of the reciprocity 
theorem, as given by Burgess, shows that the voltage 
measured across the input terminals of a linear receiving 
antenna becomes proportional to the integral 


ih RRO A OLE 


where ds is a line element of the receiver, Eine is the 
local value of the component parallel to it of the undis- 
turbed incident field we want to measure and J;, is the 


18 Stevenson, Jour. Appl. Phys., vol. 24, pp. 1134, 1143; 1953. 
4 Burgess, Wireless Engr., vol. 21, pp. 154; 1944. 
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current that would exist if the receiver were used as a 
transmitter energized at its input terminals. 

The evaluation by Woonton (SM 56) of this integral, 
in the case of a received field depending on a Kirchhoff 
integral, shows that the contributions arriving from the 
different parts of an aperture are modified by a weight- 
ing factor depending directly on the amplitude and 
phase of the radiation pattern in the direction of arrival. 
This plausible result is particularly important in meas- 
urements near and in the aperture, such as those possible 
only for microwaves. The directions of arrival from the 
different parts of the aperture vary greatly; this in- 
fluences a significant effect on the voltage induced in the 
receiver, which Borts and Woonton were able to verify 
by empirical experiments. 
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The preceding integral has to be replaced by the cor- 
responding three-dimensional integral 


fff eine tonar 


in the general case of nonlinear receiving antennas. On 
the basis of the reciprocity theorem, this expression 
equals the equivalent one over the conducting parts of 
the source if the roles of transmitter and receiver have 
been interchanged. The integral has been termed “reac- 
tion” by Rumsey (SM 2), who was able to derive a 
variational principle for its numerical evaluation. The 
previously mentioned determination of radar cross sec- 
tions is an example of such a variational principle. 


Part III—Recent Researches on the Foundations of Geometric 
Optics and Related Investigations in Electro-magnetic Theory 
E. WOLFt 


INTRODUCTION 


HIS report surveys the researches carried out in 

recent years on the foundations of geometric 

optics and discusses some of their implications in 
rigorous electromagnetic theory. 

Geometric optics has as a rule been presented inde- 
pendently of Maxwell’s theory, as a system based on 
Fermat's principle of a stationary path. According to this 
principle, light is propagated along certain curves (the 
light rays), which are the solution of the variational 


problem 
sf nds = 0, (1) 


n being the refractive index of the medium. The surfaces 
orthogonal to these curves are the geometric wavefronts 
(say S(r)=constant) and satisfy the eikonal equation 


(VS)? = n?. (2) 


Although an approach based on (1) or (2) may be 
completely satisfactory for some problems, for example, 
in estimating the approximate performance of a light- 
optical or microwave-optical system, it is well known to 
be unsatisfactory in others, and the knowledge of the 
range of validity of the geometric optics approximation 
is therefore of utmost importance. 

t Formerly, Dept. of Mathematical Physics, Edinburgh Uni- 


versity, Edinburgh, U. K. Now at the Phvsical Labs., University of 
Manchester, Manchester, England. ; ‘sila 


It is clear that to determine this range, the connection 
between geometric optics and electromagnetic theory 
must be investigated. 

Although the basic laws concerning rays and wave- 
fronts were known for a very long time, it was not until 
1911 that Sommerfeld and Runge,! following a sugges- 
tion of P. Debye, derived the eikonal equation from the 
wave equation, in the limit of infinitely high frequency 
(negligible wavelength). More precisely, what they 
showed was this, that if 


f(r) — a(r)e**S@) (3) 


is a solution of the time-independent scalar wave equa- 
tion 
V2f + n?k*f = 0, (4) 


then as k = 217/X=w/c— ~, (4) goes over into the eikonal 
equation, (2). This confirmed Kirchhoff’s well-known 
theorem that geometric optics is the limit of physical 
optics for negligibly small wavelengths. 

The Sommerfeld-Runge paper was an important step 
toward the understanding of the connection between 
geometrical optics and Maxwell’s theory. But it was 
only a first step, for their analysis was based on the 
scalar wave equation as opposed to the vector equations 
that govern the behavior of an electromagnetic field. It 
was not until the last ten years or so that the analysis of 

‘A. Sommerfeld and J. Runge, “Anwendung der Vectorrechnung 


auf die Grundlagen der geometrischen Optik,” Ann. Phys., vol. 35 
pp. 277-298; 1911, P ys., vol. 35, 
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Sommerfeld and Runge was extended to take the vec- 
torial aspect into account. These investigations, which 


led to a number of interesting new results, will now be 
surveyed. 


RESEARCHES ON THE FOUNDATIONS 
OF GEOMETRIC Optics 


Friedlander, in 1947,? suggested that when the wave- 
length is sufficiently small, a time-harmonic electro- 
magnetic field may be approximated in the form 


E = e(r)etl*S (r)—wt] | H = h(r)etltS @)—et] (5) 


where S is the eikonal, and he showed that the intensity 
law of geometric optics then follows from Maxwell’s 
equations in the limit as k>~«. Moreover, Friedlander 
found that the amplitude vectors e and h obey a simple 
propagation law along each ray. In a somewhat different 
notation, Friedlander’s equation may be written as 


de 1 1 

Sao Sa nY (—vs)e + (e-Vlog n) VS = 0; 

Or 2 be 

Oh 1 1 

any See (—vs)s + (A-vlogn) VS =0. (6) 
ai € 


Here 7 denotes a certain parameter that specifies the 
position along each ray: 


na = VS-Y, (7) 

Or 
€ is the dielectric constant and yu is the magnetic per- 
meability. Thus, if eand A are known at one point of 
a ray, their values at any other point on the ray may be 
determined by integrating (6). Eq. (6) is thus a trans- 
port equation for the vector amplitudes e and A and 
gives information about the change in the state of 
polarization along each ray. 

Friedlander also suggested that (5) should perhaps 
be regarded as first terms in the asymptotic expansion 
of the field vectors, in the form (omitting as usual the 
time harmonic factor) 


i) 1 oo 1 
u ~ etks e,, y ~ eiks “> Asynih8) 
~ (ik)” yao (1k)” 


a suggestion confirmed for a wide class of fields by other 
investigations about which more will be said later. 

The relations in (6) have an interesting implication in 
terms of non-Euclidian geometry. With a medium speci- 
fied by the refractive index function u(r) =~V/eu one 
may associate a non-Euclidian space whose line element 
ds’ is given by 


ds’ = nds = nv/dx? 4+ dy? + dz’. (9) 


From (1) and (9) it follows that the geometric rays in 
ordinary space correspond to the geodesics in the non- 


2, G. Friedlander, “Geometrical optics and Maxwell’s equa- 
tions,” Proc. Camb. Phys. Soc., vol. 43, pp. 284-286; 1947. 
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Euclidian space, and it may be shown that (6) implies 
that each of the directions of polarization remains 
parallel to itself along each geodesic in the sense of Levi- 
Civita. This implication of (6) was pointed out by 
Luneburg,’ although the result was stated much earlier 
in a somewhat different form by Bortolotti.4 

The following authors studied the connection between 
Maxwell’s theory and geometric optics in anisotropic 
and, in some cases, nonconducting media: Arley,® 
Suchy,® Mandl,’ and Cap and Réver.’ 

The paper by Arley is concerned with the derivation 
of Huygens’ principle and Fermat’s principle. This 
paper, as well as that of Suchy and Mandl, considers 
also the limits of validity of geometric optics. Cap and 
Réver were concerned mainly with the formulation of 
Fermat’s principle for an absorbing medium. 

Some of the results obtained in these papers are not 
entirely satisfactory. For example, Arley and Mandl 
overlooked the fact that in a time-harmonic vector field 
it is necessary to introduce, not one, but three phase 
functions, one for each Cartesian component. Also, the 
conditions obtained by Arley and Suchy for the validity 
of geometric optics do not appear to be all independent 
of each other, and are probably too strong in some cases. 
According to their results, geometric optics can be a 
good approximation only if the fractional changes of the 
refractive index over a linear domain of the order of a 
wavelength is small compared to unity. If this result is 
true, it is difficult to understand how the law of refrac- 
tion can be valid, because across a refracting surface the 
refractive index changes discontinuously by an appre- 
ciable amount. Actually, no satisfactory derivation of 
the geometric optics law of refraction from Maxwell’s 
theory appears to have been given so far. 

Luneburg*? and Born and Wolf® derived explicit ex- 
pressions, in terms of the eikonal, for the variation of the 
intensity along each ray. If J, and J: are the intensities 
at two points of the ray, and m, and mz are the corre- 
sponding values of the refractive index, then 


Tp n> S225 

7 =e (-f — ds) = 

Tq, Ny Sy n? 

(10) 


3R. K. Luneburg, “Propagation of Electromagnetic Waves,” 
Notes of lectures delivered at New York University, 1947-1948. 

4E. Bortolotti, “Sulle rappresentazioni conformi, e su di una in- 
terpretazione fisica del parallelismo di Levi-Civita,” Rend. Acc. Linc., 
vol. 4, pp. 552-556; 1926. : ; 

5N. Arley, “A note on the foundations of geometrical optics,” 
Det. Kgl. Danske Vidensk. Selsk., vol. 22; pp. 1-21; 1945. : 

6K. Suchy, “Schrittweiser iitbergang von der wellenoptik zur 
strahlenoptik in inhomogenen anisotropen absorbierenden medien,” 
Ann. Phys., vol. 11, pp. 113-130; 1952; vol. 12, p. 423; 1953; vol. 13, 
pp. 178-197; 1953. ‘ 

7G. Mandl, “Zur begriindung der strahlenoptik aus der Max- 
wellschen feldtheorie,” Acta Phys. Aus., vol. 7, pp. 365-389; 1953. 

8F, Cap and W. Rover, “Lichtwege in inhomogenen, absorbier- 
enden, isotropen medien,” Acta Phys. Aus., vol. 8, pp. 346-355; 
1954, 

9M. Born and E. Wolf, “Principles of Optics,” Pergamon Press, 
London, England. In press. 
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the integral being taken along the ray. Similar relations 
also hold for the variation of the squared amplitudes of 
electric and magnetic vectors. In the special case when 
wavefronts are concentric spheres, in a homogeneous 
medium, (10), of course, reduces to inverse-square law. 

A fuller account of results of these researches will ap- 
pear in a forthcoming book by Born and Wolf.® 


PuLSE SOLUTIONS OF MAXWELL’S EQUATIONS 
AND ASYMPTOTIC EXPANSIONS OF 
STEADY-STATE FIELDS 


It is well known that the eikonal equation, (2) may 
also be regarded as the equation of the characteristics of 
the wave equation and describes rigorously the propa- 
gation of discontinuities in the solutions of these equa- 
tions. Physically, the discontinuous solutions represent 
the following situation: 

A source begins to radiate at a particular instant 
t=to. The field then spreads into the surrounding space, 
and at any later time t>?o, the disturbance reaches a 
well-defined region. Let 


PAT eal (11) 


be the equation of the boundary surface of this region. 
Across this surface the field changes discontinuously 
from a finite value in the interior to zero value outside it. 
The statement that the eikonal equation is the charac- 
teristic equation of the wave equation is physically 
equivalent to the statement that the surfaces F(p) =con- 
stant obey the eikonal equation 


(VE) = 97. (12) 


Hence, there exists a formal analogy between the 
steady-state field at small wavelengths and the field of 
discontinuities (representing a so-called pulse solution), 
latter being a rigorous solution of Maxwell’s equations. 

In important investigations carried out by Lune- 
burg*° and extended by Kline,! it was found that in 
some cases it is possible to derive formally the whole 
asymptotic expansion for a steady-state field from the 
knowledge of the associated pulse solution. We shall 
briefly indicate the main results of these investigations. 

Let the sources of the field be specified by a function 
Q(r, t) in terms of which the charge density p and the 
enforced current density j are given by 


1 owe 
=-—vQ, j=—@. 13 
p es Ome ine 2 Gish 


If the charges remain at a fixed location in space but 
vary in strength with the time, as is usually the case, @ 
will be of the form 


 R. K. Luneburg, “Mathematical Theory of Optics,” Notes of 
lectures delivered at Brown University; 1944. Also “Asymptotic De- 
velopments of Steady-State Electromagnetic Fields,” New York 
University Math. Res. Group. Rep. No. EN. 14; 1949. 

"M. Kline, “An asymptotic solution of Maxwell’s equations, ” 
Comm. Pure and Appl. Math., vol. 4, 225-262; 1951. Also “Sym- 
posium on Theory of Electromagnetic Waves,” Interscience Pub- 
lishers, New York, N. Y., pp. 225-262; 1951. 
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Q(r, t) = K(r)T(2). (14) 
In particular, for a pulse solution, @= Qo, where 
Q.(r, t) = 0 when t < ‘4 
= Kir) whenites.0; (15) 


In terms of the pulse solution Eo, Ho, the field EZ, H 
that arises from a @ distribution with any prescribed 
time dependence 7(t), may be derived from the Duhamel 
principle: 


art 
Ea if E(i—v) Tat: 
Aah 


H = < f aw — 1) T(t)dt'. | (16) 
dt 0 


Using this result, Luneburg and Kline showed that 
when there is a time-harmonic dependence T(t) =e’, 
with increasing time E and H approach the form 


B= ure <. Hi vine e (17) 


where u and v may be developed into the following 
asymptotic series” 


1 1 
ua Kile + oF fey — 1 +7 ee a 
iw (iw)? 
; 1 1 
y= 8 Ng — — hs + : shy by (18) 
iw (iw)? 
where 


0’ Ey 0’ Ho 
e, = and A, (19) 
ol” or” 
are the changes in the successive derivatives of Eo and 
Hy across the discontinuity surface and 


(20) 


Hence, each coefficient in the asymptotic expansion of a 
steady-state field is given by the discontinuous change in a 
derivative of a field that represents a certain pulse solution 
of Maxwell's equations. 

For sufficiently large values of w, in the region where 
kK=0, 

pS, [Ey leit*F—-et]; Tis [Hp JettF—-et, (21) 

This confirms the assumption of Friedlander about the 
form of E and H. See (5). 


Luneburg and Kline also found that each term in the 


expansion of (18) satisfies the following transport equa- 
tions: 


“It is assumed here that only one discontinuity surface passes 
through each point. In some cases, such as when reflection takes place 
at bodies embedded in the medium, several such surfaces may pass 
through some points. Eq. (18) must then be modified, each term re- 


placed by the sum of such terms, each of which arises from one of the 
partial fields. 
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Be) 1 1 
3 +5uv-(—vr)e,+ le, vlog n]vF = —C, 


70) 
Bye nat 1 
Bee, ¢Y (G vr) h, + [h,- Vlog n|F = — D,, (22) 
where 
Co = Dy) = 0 
and for y>0, 
1 
C, = —c <uv X (—v xe) 
Me 
1 
— +[— v-(ce.) |! ; 
€ 
1 1 
Di=——c Jew x (—vx hs) 
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1 
- v|— V- (uh,_) } (23) 
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T being as before a parameter that specifies the position 
on each “ray”: 


= = VF -Y. (24) 
OT 
In the special case v=0, (22) is identical with (6). 

In addition to the work of Luneburg and Kline, in 
which the transport equations for the pulse solutions 
were derived by using the full Maxwell’s equations, two 
related investigations must be mentioned. Copson' 
showed that the transport equations for e9 and Ay can 
be derived from the second-order wave equations for the 
electric and magnetic vectors. Bremmer,' in a very 
elegant analysis, further demonstrated the mathemati- 
cal equivalence of the geometric optics solution and the 
pulse solution by using the delta function techniques. 

The connection between the discontinuous solutions 
and the geometric optics field was also studied by 
Mandl’ in the paper already referred to. 

Since there is a tendency in the literature to identify 
the pulse solution with the geometric optics solution, it 
should be stressed that this identity is a purely formal 
one, based on the mathematical equivalence of the cor- 
responding equations. From the point of view of physzcs, 
the two solutions are of course quite distinct. The 
steady-state field, where geometric optics is a high- 
frequency limit, represents a wave that exists for an 
infinite length of time. On the other hand, the disconti- 


_ nuity surface of the pulse solution passes through each 


point in an infinitesimal instant of time—and is in a 
sense somewhat fictitious, for one can never measure the 
field on such a moving surface. Nevertheless, the mathe- 
matical analogy demonstrated by the work of Luneburg 


13 —. T. Copson, “The transport of discontinuities in an electro- 
magnetic field,” Comm. Pure and Appl. Math., vol. 4; pp. 427-433; 
1951. ; 

“4 YH, Bremmer, “The jumps of discontinuous solutions of the 
wave equation,” Comm. Pure and Appl. Math., vol. 4; pp. 419-426; 
1951. 
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and Kline gives us a much clearer insight than was 
hitherto possible into the nature of the geometric optics 
approximation. It is also likely to prove of value in 
other problems concerned with high-frequency electro- 
magnetic fields. 


SOME ANALOGIES BETWEEN GEOMETRIC OPpTIcs 
AND RIGOROUS ELECTROMAGNETIC THEORY 


Unlike the geometric optics approximation, the pulse 
solution represents a rigorous solution of Maxwell’s 
equations. The question naturally arises whether some 
of the geometric optics laws do not also hold, with ap- 
propriate modifications, in time-harmonic fields of an 
arbitrary frequency. 

An essential property of the geometrical optics ap- 
proximation is the orthogonality between the energy 
flow and a family of cophasal surfaces, and one must 
therefore consider whether such a relationship of orthog- 
onality is found in any of the known rigorous solutions 
of Maxwell’s equations. 

It was first pointed out by Braunbec® and illustrated 
in beautiful diagrams by him and Laukien" that orthog- 
onality of this type exists in the classic half-plane 
problem of Sommerfeld (see Fig. 1). If we bear in mind 
the complexity of the field near the diffracting edge the 
relationship of orthogonality seems to be a rather sur- 
prising feature. 

Braunbek" investigated whether such a relationship 
exists in more general fields and came to the conclusion 
that only in a limited class of fields (two-dimensional, or 
those with rotational symmetry) are the lines of the 
energy flow orthogonal to a family of cophasal surfaces. 

Somewhat more general results in this direction were 
recently obtained by Nisbet and Wolf!” in studies of 
electromagnetic fields where the waves are of arbitrary 
shape but where at least one of the field vectors is 
linearly polarized.1* They found that such fields show a 
close resemblance to the field in geometric optics. For 
example, if the linearly polarized vector is the electric 
vector E and we write 


E = e et(kS-#t), 


(25) 


where e is a real vector function of position, the energy 
flow is at each point found to be along the normal to the 
cophasal surfaces S(r) =constant, and the “electromag- 
netic rays” (lines of energy flow) so defined satisfy the 
generalized Fermat’s principle 


5 [ vas = 0, 


1 W. Braunbek, “Zur darstellung von wellenfeldern,” 
Naturf., vol. 6a, pp. 12-15; 1951. 

16 W. Braunbek and G. Laukien, “Einzelheiten zur halbenen- 
Beugung,” Optik, vol. 9, pp. 174-179; 1952. 

17 A, Nisbet and E. Wolf, “On linearly polarized electromagnetic 
waves of arbitrary form,” Proc. Camb, Phil. Soc., vol. 50, pp. 614-622; 
1954. 

18 This is to say that with varying time, the field vector at any 
particular point remains always in the same direction; unlike in the 
case of a homogeneous plane wave, however, this direction may be 
different at different points of the field. It has been shown by Nisbet 
and Wolf that many fields of practical interest are of this type. 


(26) 
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where 


yin nt —e [yx (—vx e) |. (27) 
ek? mM 

As k> ©, N—n, and (26) then correctly reduces to the 
ordinary Fermat’s principle. Nisbet and Wolf also 
showed that along the “electromagnetic rays” the vector 
amplitude e rigorously obeys a transport equation of 
the form of (6). 


Fig. 1—The classic Sommerfeld half-plane diffraction problem (H- 
, polarization). (a) Curves of constant phase of the H-vector; 
(6) Lines of the average energy flow (Poynting vector). (After 
Braunbek and Laukien."*) 


Another approach to the question of finding a simple 
model for the transport of energy arises from the fact 
that the field vectors do not uniquely define the flow. 
The usual definition is in terms of the Poynting vector 


C 
4a 
which obeys the conservation law 
dw 
By Thy Oe (29) 
where 
1 
= —(E-D+ H-B) (30) 
8a 
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is the energy density and 


Q=cEF’ (31) 
is the Joule heat, « being the conductivity. From (29) it 
is evident that the flow could also be defined by 


S’=S+4+V XA, (32) 
where A is any arbitrary vector, since on account of the 
identity V- (VX A) =0, the A vector contributes nothing 
to the total outward flux from any closed domain. 

__ According to modern field theories, freedom in defin- 
ing the flow density and the energy density is still 
greater; any definition consistent with the addition of a 
four-divergence to the Lagrangian density L of the field 


b= —(ED-H-B) (33) 


is permissible. As long ago as 1902, definitions of w and 
S which differ from the usual ones in this way were pro- 
posed by McDonald.!* Recently, Hines?® asserted that 
McDonald’s definition of the flow shows definite ad- 
vantages over that of Poynting’s. 

Green and Wolf”! recently showed that in vacuum, in 
regions free of currents and charges, the field may be de- 
rived from a single complex scalar wave function V(r, f), 
and that in terms of this function the energy density and 
the flux density may be defined by means of expressions 
analogous to the quantum-mechanical expressions for 
the probability density and the probability current. 

The class of transformations introduced by Green and 
Wolf was studied in detail by Bremmer,” who also in- 
vestigated the relationship between the V function and 
other potential functions previously considered by 
E. T. Whittaker and R. K. Luneburg. Nisbet recently 
showed that the real and imaginary parts of the V func- 
tion are closely related to certain Hertz vectors. A 
paper by Nisbet”? fully clarifies the old question relating 
to the number of components of a Hertz vector needed 
in the solution of different types of electromagnetic 
problems. 

From the very brief account of the researches dis- 
cussed in this report, it will be apparent that the subject 
has attracted a good deal of attention in recent years 
and that these investigations have considerably clarified 
our understanding of some general properties of electro- 
magnetic fields. 


* H. M. MacDonald, “Electric Waves,” Cambridge University 
Press, Cambridge, p. 72; 1902. See also “Electromagnetism,” Bell & 
Sons, Ltd p. a7 1 9oa, 

°C, O. Hines, “Electromagnetic energy density and flux,” 
Canad. Jour. Phys., vol. 30, pp. 123-129; 1952. 

*1 H. S. Green and E. Wolf, “A scalar representation of electro-— 
magnetic fields,” Proc. Phys. Soc., vol. 66, pp. 1129-1137; 1953. 

2H. Bremmer, “Remarks on the complex scalar function in- 
troduced by H. S. Green and E. Wolf in the theory of electromagnetic 
fields,” Philips Res. Reps. In Press. 

*8 A. Nisbet, “Hertzian electromagnetic potentials and associated 
gauge transformations,” Proc. Roy. Soc. A, vol. 231, pp. 250-263; 
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